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ABSTRACT 


At  16  Me / s ,  the  backscatter  coefficient  for  sea  is 
found  to  be  10  dB  higher  than  for  land  for  angles  of  ele¬ 
vation  between  25°  and  15°,  At  lower  angles  there  is  a 
knee  effect,  and  the  backscatter  coefficient  decreases  very 
rapidly.  The  knee  angle  is  lower  for  sea  than  for  land. 

For  a  given  surface,  and  at  angles  above  the  knee,  the  var¬ 
iation  of  backscatter  coefficient  with  elevation  angle  is 
consistent  with  the  type  of  scatter  expected  from  upright 
objects  such  as  trees  or  wave  crests,  and  the  fading  of  the 
echoes  is  ascribed  to  the  Doppler  movement  of  these  objects. 

The  frequency  of  observation  of  Sporadic-E  by  back¬ 
scatter  sounding  is  strongly  influenced  by  whether  the  back¬ 
scatter  occurs  on  the  land  or  on  the  sea.  Sporadic-E  itself 
appears  to  be  uniformly  distributed. 

In  observations  of  F  region  propagation  by  the  back¬ 
scatter  technique,  the  nature  of  the  terrain  should  be  taken 
into  account.  This  is  particularly  important  in  the  case  of 
trans-equatonal  one -hop  propagation. 

The  experimunta]  procedures  developed  here  include  a 
new  method  of  measuring  the  vertical  radiation  patterns  of 
large  high  frequency  antennas,  a  computer  program  to  calcul¬ 
ate  the  power  of  backscatter  echoes,  and  the  application  of 


vertical  incidence  data  to  oblique  propagation.  New  insights 
are  gained  into  the  interpretation  of  backscatter  records. 
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BACKSCATTER  OF  RADIO  WAVES  FROM  THE  GROUND 

1.  INTRODUCTION 

Ground  backscatter  is  the  process  by  which  high  fre¬ 
quency  waves  are  scattered  at  the  ground,  so  that  some  of 
the  energy  goes  back  along  its  original  path  to  the  trans¬ 
mitter.  This  process  is  used  in  oblique  incidence  sounding 
of  the  ionosphere.  A  pulse  of  radio  energy  may  undergo 
backscatter  at  the  ground  after  one  or  more  hops  via  the 
ionosphere,  and  the  time  the  echo  takes  to  return  is  a  meas¬ 
ure  of  the  range  at  which  backscatter  occurs.  The  various 
uses  of  backscatter  sounding  are  outlined  in  PETERSON  (1954) . 

Calculations  of  the  strength  of  backscatter  echoes  have 
been  hampered  by  the  lack  of  information  about  the  back¬ 
scatter  coefficient.  This  coefficient  is  a  parameter  sim¬ 
ilar  to  the  radar  cross-section  per  unit  area.  It6  relation 
to  other  parameters  will  be  considered  in  Chapter  4. 

The  composition  and  roughness  of  the  ground  where 
backscatter  takes  place  may  be  expec1 °d  to  influence  the 
backscatter  coefficient.  An  understanding  of  the  variation 
of  the  backscatter  coefficient  with  elevation  angle  and  type 
of  ground  is  essential  to  the  correct  interpretation  of  back¬ 


scatter  records,  and  is  the  main  concern  of  this  proiect. 
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(a)  Variation  with  angle  of  elevation  For  a  given 

backscatter  sounder,  a  particular  type  of  echo  is  not  norm¬ 
ally  detectable  beyond  a  certain  maximum  range,  and  this 
range  may  be  determined  by  the  angular  variation  in  the  back¬ 
scatter  coefficient.  A  knowledge  of  this  variation  facil¬ 
itates  the  identification  of  echo  types  and  the  recognition 
of  anomalous  behaviour  of  such  types.  This  will  be  discussed 
with  special  reference  to  Sporadic-E  echoes. 

(b)  Variation  with  type  of  ground  :-  Neglect  of  any 

variation  with  ground  type  may  result  in  misleading  conclus¬ 
ions  on  the  nature  of  the  ionosphere  in  certain  directions. 
For  example,  as  a  result  of  backscatter  observations,  it  has 
been  thought  that  occurs  more  frequently  over  the  sea 

than  over  the  land  (EGAN  and  PETERSON,  1961,  1962) ,  and  that 

transequatorial  one-hop  propagation  occurs  preferentially  to 
the  north-east  of  Brisbane  rather  than  to  the  north-west 
(THOMAS,  1962).  Knowledge  of  the  backscatter  coefficient 
may  cause  these  conclusions  to  be  modified. 

In  the  present  investigation,  the  16  Me / s  backscatter 
sounder  at  Brisbane  was  used  (THOMAS  and  McNICOL,  1960a, 
1960b,  1962) .  It  consists  of  a  rotating  array  of  four 
horizontally-polarized  Yagi  antennas  mounted  half  a  wave¬ 
length  above  the  ground  (Fig.  1).  The  antennas  were  directed 
to  the  east  and  to  the  west,  to  obtain  echoes  from  the  sea 
and  the  land  in  turn.  Rectangular  pulses  of  600  microseconds 
duration  were  transmitted,  at  a  peak  power  of  5  kW .  Back- 
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PiC.  1 


16  He/ s  rotating  array 
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scatter  echoes  were  recorded  on  a  range  — amplitude  display, 
by  the  use  of  a  swept-gain  unit.  As  the  recording  film 
moved  slowly  past  the  brightness-modulated  oscilloscope  trace 
the  echo  intensity  in  dB  was  plotted  against  the  oblique 
range  of  the  echo. 

For  comparison  with  these  range — amplitude  records,  the 

range - amplitude  relationship  was  calculated,  assuming  the 

backscatter  coefficient  to  be  independent  of  the  angle  of 
elevation  of  the  rays.  The  calculation  involved  ionospheric 
parameters,  estimated  from  vertical  incidence  soundings  at 
Brisbane,  and  the  vertical  radiation  pattern  of  the  antenna 
array,  measured  by  means  of  a  balloon-borne  transmitter. 

The  results  of  the  calculations  were  compared  with  the 
corresponding  observed  backscatter  records  to  deduce  the 
relative  backscatter  coefficient  at  the  calculated  angle  of 
elevation . 
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2.  PREVIOUS  HORK 


2.1  Variation  with  angle  of  elevation 

DIEMINGER  (1951)  ,  experimenting  at  Lindau  between  1  and 
20  Mc/s,  found  that  the  echo  amplitude  decreases  rapidly  with 
range,  and  reasoned  that  this  could  only  occur  if  the  back- 
scatter  coefficient  depends  mainly  on  the  number  and  areas  of 
surfaces  larger  than  a  wavelength.  As  such  areas  usually 
have  a  very  small  angle  of  tilt,  the  backscatter  coefficient 
should  decrease  very  rapidly  as  the  angle  of  elevation 
decreases . 

SHEARMAN  (1961)  at  Slough  recorded  15  and  21  Mc/s  back- 
scatter  echoes  on  a  range — amplitude  display.  For  comparison, 
he  calculated  the  echo  patteri.  ,  neglecting  any  angular  depend¬ 
ence  of  backscatter  coefficient,  and  making  other  simplifying 
assumptions,  and  found  fairly  good  agreement  with  observat¬ 
ions.  In  reply  to  discussion  he  stated  that  the  major  unknown 
was  the  variation  of  backscatter  with  the  angle  of  incidence, 
and  that  this  could  only  be  measured  satisfactorily  from  an 
aircraft.  Meanwhile,  NIELSON  et  al .  (1960)  published  some 

results  of  measurements  made  with  an  airborne  radar.  They 
pointed  out  that  when  ground-based  backscatter  sounders  are 
used  to  find  ground  backscatter  parameters,  at  least  two  ion¬ 
ospheric  reflections  are  involved,  and  the  precise  nature  of 


these  reflections  cannot  be  determined.  To  avoid  these  ion¬ 
ospheric  reflections,  they  installed  a  4  kW  backscatter  radar 
operating  at  32.8  Me / s  in  an  aircraft,  which  was  flown  over 
California,  the  Pacific  Ocean,  Alaska  and  the  Arctic  Ocean. 
Echoes  were  recorded  on  an  A-scope  display,  and  compared  with 
echo  strengths  calculated  on  the  basis  of  a  constant  back- 
scatter  coefficient  of  one.  The  ratio  of  the  observed  to  the 
calculated  values  was  taken  as  a  measure  of  the  backscatter 
coefficient.  They  constructed  curves  of  backscatter  coeffic¬ 
ient  against  angle  of  elevation  from  60°  to  4°,  showing  a 
slope  of  about  0.6  dB  per  degree  for  horizontally  polarized 
waves.  At  angles  below  20°,  the  slope  was  greater  than  0.6 
dB  per  degree  for  ice,  and  less  than  0.6  for  desert.  The 
antennas,  which  were  attached  to  the  nose  and  tail  of  the 
aircraft,  were  difficult  to  calibrate  accurately,  as  the  cal¬ 
ibration  had  to  be  done  while  the  aircraft  was  in  flight, 
using  a  receiver  on  the  ground.  It  was  difficult  to  know 
the  aircraft's  position  and  orientation  during  the  calibrat¬ 
ion,  and  the  overall  errors  amounted  to  perhaps  6  dB  (NIELSON 
ft  al .  page  34,  errata)  ,  and  therefore  12  dB  for  transmission 
and  reception. 

HAGN  (1962)  analysed  t ne  work  of  NIELSON  et  al .  (1960). 

Ife  published  a  i evised  set  of  curves  lor  the  backscatter 
coefficient.  He  rejected  results  below  angles  of  10°  for 
sea  and  20°  for  land,  because  at  lower  angles  the  curves 
show  a  "knee"  effect  (HAGN,  1962,  p.8  and  p.85)  where  the 
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backscatter  coefficient  decreases  very  rapidly.  This  knee 
did  not  always  occur  at  the  same  angle,  even  for  the  same 
terrain.  Further,  the  presence  of  mountains  at  the  horizon 
often  made  it  impossible  to  assume  that  the  ground  was  uni¬ 
formly  rough  for  rays  near  the  tangent  ray.  Towards  higher 
angles  of  elevation,  the  backscatter  coefficient  increased 
rapidly  for  horizontal  polarization  (about  1  dB  per  degree) , 
but  decreased  very  slowly  for  vertical  polarization  (about 
-0.1  dB  per  degree) . 
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2.2  Variation  with  type  of  ground 

VILLARD  and  PETERSON  (1952a,  1952b),  at  3-30  Me / s , 
found  no  noticeable  change  in  echo  amplitude  from  land  to 
sea.  They  suggested  that,  unlike  u . h . f .( 300-3000  Mc/s) 
echoes  which  do  show  some  change,  h.f,(3-30  Mc/s)  back- 
scatter  echoes  involve  comparatively  long  wavelengths  and 
comparatively  enormous  echoing  areas.  Both  these  factors 
would  tend  to  prevent  any  discrimination  between  land  and 
sea  in  backscatter  records. 

DIEMINGER  (1951) ,  while  observing  backscatter  echoes 
at  1-20  Me / s  swept  frequency,  identified  an  echo  group  which 
remained  at  almost  constant  range,  in  contrast  to  the  group 
associated  with  the  skip  focusing  area,  which  was  varying  in 
range.  The  range  of  the  fixed  group  corresponded  to  back¬ 
scatter  from  the  Northern  slopes  of  the  Alps.  SILBERSTEIN 
(1954) ,  using  3-25  Mc/s  swept  frequency  backscatter,  found  a 
similar  group  of  echoes  constant  in  range,  thought  to  be  due 
to  the  Rocky  Mountains.  Me CUE  (1956)  found  that  at  angles 

of  elevation  near  30°  land  is  a  more  prominent  source  of 
backscatter  than  sea,  for  frequencies  near  6  Mc/s. 

SHEARMAN  (1956b)  reported  that  observations  at  Slough 
between  10  and  27  ac / s  confirmed  the  lack  of  dependence  on 
terrain  found  by  VILEARD  and  PETERSON.  He  considered  that 
this  was  surprising,  in  view  of  the  large  difference  between 
sea  and  land  echoes  noticed  with  airborne  centimetric  radar 
equipment.  He  suggested  that  any  differences  are  masked  by 
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the  large  variations  in  echo  amplitude  introduced  by  the  ion¬ 
osphere.  In  the  discussion  following  this  paper  it  was 
pointed  out  (by  W.  R.  Piggott)  that  the  ordinary  countryside 
is  quite  rough  from  the  radio  point  of  view,  and  that  a  moun¬ 
tain  is  not  necessarily  a  better  echoing  object  than  the 
ground  as  a  whole.  The  only  cases  which  have  been  reported 
of  abnormal  backscatter  from  distant  mountain  masses  are 
those  for  which  the  phase  of  the  reflected  wave  is  held  con¬ 
stant  over  a  wide  area,  so  that  we  have  a  coherent  mirror- 
type  reflection.  For  example,  the  ionospheric  station  at 
Lindau  is  almost  at  the  centre  of  curvature  of  the  Alps,  and 
frequently  receives  an  abnormal  echo  coming  from  that  arc. 

The  beginning  of  the  Alps  forms  a  sort  of  curved  mirror,  and 
while  the  mountains  themselves  do  not  contribute  much,  the 
gradual  rise  at  their  beginning  does.  In  contrast,  no  echo 
from  the  Alps  was  received  at  Slough. 

SHEARMAN  considered  vertically  polarized  waves,  and 
developed  a  model  for  the  scattering  source,  namely  a  system 
of  hemispherical  bosses  on  a  perfectly  conducting  ground,  and 
found  that  a  density  of  about  1000  bosses  /km2  could  account 
for  the  observed  strength  of  a  backscatter  echo  from  a  ground 
range  of  670  km  West  of  Slough,  that  is,  from  the  Atlantic 
Ocean.  In  the  discussion  it  was  noted  that  these  scattering 
sources  would  be  common  for  average  sea  conditions,  and  a 
serious  decrease  in  scattering  would  only  be  expected  for 
exceptionally  calm  conditions.  SHEARMAN  concluded  the  di s- 


c  us s ion  by  stating  that  there  was  a  clear  need  for  a  thorough 
investigation  of  the  relative  amplitude  of  scatter  from  land 
and  sea,  and  that  the  direct  measurement  from  an  aircraft  of 
the  scattering  properties  of  various  types  of  terrain  at 
these  wavelengths  would  provide  much  useful  evidence. 

WILKINS  and  SHEARMAN  (1957)  reasserted  the  view  that, 
in  general,  terrain  effects  tend  to  be  masked  by  much  larger 
variations  of  echo  strength  due  to  the  variable  efficiency 
of  the  ionosphere  as  a  propagation  medium.  They  added  that 
with  the  wide  beams  in  use  with  rotating  aerials,  azimuthal 
resolution  of  geographical  features  is  unlikely  to  be  pract¬ 
icable.  SHEARMAN  (1961)  said  that  the  variation  with  the 
type  of  ground  appears  surprisingly  small,  and  that  the  only 
ob-jects  which  have  so  far  been  identified  are  large  mountain 
ranges . 

RANZI  and  DOMINICI  (  1959)  ,  although  using  a  wide  beam 
observed  that  22,3  Mc/s  backscatter  echoes  from  continents, 
and  especially  from  the  Sahara  Desert,  have  a  considerably 
lower  intensity  than  those  from  the  sea.  They  suggested  a 
difference  of  about  10  dB  between  land  and  sea.  DOMINICI 
(1962,  1963)  qualified  this  by  remarking  that  the  intensity 

distribution  of  echoes  was  influenced  by  the  variation  in  the 
vertical  radiation  pattern  of  the  antenna,  owing  to  the  vary¬ 
ing  topography  around  the  antenna. 

NIELSON  et  al.  (I960)  found  that  for  horizontal  polar¬ 
ization,  and  at  32.8  Mc/s,  the  backscatter  coefficient  for 


desert  was  about  b  dB  less  than  that  for  smooth  ice  or  sea 
at  angles  of  elevation  above  30°,  but  that  at  about  10° 
desert  gives  stronger  backscatter  chan  sea.  Vertical  pol¬ 
arization  gives  stronger  return  than  horizontal  from  sea, 
but  land  is  indifferent  to  polarization.  As  the  roughness 
of  the  sea  increases,  polarization  effects  become  less  not¬ 
iceable  . 

HAGN  (1962)  ,  in  revising  the  results  of  NIKLSON  et 
al . ,  found  that  for  both  horizontal  and  vertical  polariz¬ 
ations,  the  backscatter  coefficient  for  land  is  about  20 
dB  below  that  for  sea .  Owing  to  the  uncertainty  concern¬ 
ing  the  knee  in  the  curves,  no  comparison  could  be  made  at 
angles  lower  than  20°. 

2.3  Summary 

After  many  years  during  which  ground-based  sounders 
found  no  difference  between  sea  and  land  scatter,  one  ex¬ 
ception  emerged  in  which  the  sea  appeared  to  give  stronger 
bac  k  sc  a  t  t  <  r  , 

The  widely  publicized  need  for  aircraft  measurements 
of  the  back  sea  tter  coefficient  has  been  fulfilled,  and  the 
results  show  two  major  effects,  the  existence  of  a  knee  j.n 
the  curve  of  backscatter  coefficient  against  elevation 
angle,  and  the  predominance  of  sea  scatter  strength  over 
that  of  land  scatter.  However,  only  comparatively  smooth 


surfaces  were  considered. 


Also,  the  calibration  of  air- 
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borne  antennas  proved  to  be  quite  difficult  and  the  results 
may  contain  large  errors, 
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3.  THE  RELATIONSHIP  OF  THE  PRESENT 

INVESTIGATION  TO  PREVIOUS  WORK 

The  present  work  concerns  elevation  angles  from  about 
7°  to  30°,  and  the  polarization  at  the  echoing  area  is  ellip¬ 
tical  and  variable,  owing  to  Faraday  rotation  in  the  ionos¬ 
phere  . 

Concerning  the  variation  of  the  backscatter  coefficient 
with  angle  of  elevation,  the  knee  effect  has  been  discovered 
independently  of  HAGN  (1962) ,  and  the  angle  at  which  this 
occurs  has  been  investigated  for  the  terrain  types  available. 
For  angles  above  the  knee,  the  coefficient  decreases,  and 
this  trend  agrees  with  IlAGN's  results  for  vertical  polarizat¬ 
ion. 

The  variation  with  the  type  of  ground  is  quite  marked, 
and  the  present  work  indicates  that  for  land,  the  backscatter 
coefficient  is  about  10  dB  less  than  for  sea.  This  supports 
the  jbservations  of  RANZI  and  DOMINIC  I  (1959).  For  low 
angles,  HAGN's  reappraisal  of  NIELSON  et  aJ  .  has  removed  the 
most  serious  disagreement  with  the  present  work. 

Trie  method  used  here  has  certain  advantages  over  prev¬ 
ious  met  hods . 

(a)  Narrow  Azimuthal  Beam . 

(i)  As  tlie  antenna  beam  is  narrower  than  in  previous  in- 
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investigations,  there  is  less  azimuthal  spread  of  the  rays  in 
the  ionosphere,  and  therefore  the  ionospheric  parameters  in 
the  path  of  the  rays  can  be  estimated  more  precisely. 

(ii)  The  minimum  group  range  of  a  backscatter  echo  serves 
as  e  guide  to  the  ionospheric  configuration  at  the  control 
point,  and  this  range  is  well  defined  when  the  antenna  beam 
is  narrow. 

(lii)  Owing  to  the  narrower  beam,  it  is  easier  to  discrimin¬ 
ate  between  different  areas  of  the  ground  from  which  back¬ 
scatter  is  received. 

( b)  Vertical  Radiation 
P Pat  tern  of  Antenna . 

(  1)  A  ground-based  antenna  is  easier  to  calibrate  than  an 
airborne  antenna. 

(ii)  Changes  in  the  topography  near  the  antenna  can  be  taken 
into  account,  so  that  the  variations  of  the  vertical  pattern 
with  azimuth  are  well  understood. 

(c)  Location  of  the  Sounder . 

(i)  The  sounder  is  well  placed  for  comparing  land  ar d  sea 
backscatter,  as  the  land  and  sea  extend  without  interruption 
for  several  thousand  kilometers  to  the  west  and  east 
respect ive ly . 

(ii)  It  is  particularly  easy  to  observe  the  boundary  between 
land  and  sea  scatter,  as  the  coastlines  recede  almost  rad¬ 
ially  from  Brisbane  at  ranges  beyond  500  km, 

(m)  The  absence  of  noise  such  as  that  due  to  aircraft 
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\ 


engine  ignition  facilitates  good  recording  of  echoes. 

(d)  Applicability  of  Results . 

(  1)  The  polarization  of  the  waves  at  the  echoing  area  is 
typical  for  operating  conditions  of  backscatter  sounders, 
therefore  the  backscatter  coefficient  derived  here  should  be 
more  applicable  to  normal  backscatter  sounding  than  that  de¬ 
rived  from  airborne  experiments. 

(li)  The  frequency  used  here  (16  Me / s )  is  typical  for  back¬ 
scatter  sounding,  whereas  the  results  obtained  at  higher 
frequencies  cannot  confidently  be  extrapolated  to  this  reg¬ 
ion  . 

(iii)  The  echoes  rr  ty  include  returns  from  many  large  topog¬ 
raphical  features,  as  the  echoing  area  is  large.  In  con¬ 
trast,  aircraft  measurements  involve  a  much  smaller  sample 
of  terrain,  and  therefore  comparatively  smooth  surfaces  only 
were  selected  ( HAGN ,  1962)  . 
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4.  DISCUSSION  ON  THE  BACKSCATTER  COEFFICIENT 


4.1  Def init ions 


Backscatter  coefficient 

The  radar  equation  is  of  the  form 
0  G  2  P  X  2 

P  =  - =~9— (CLAPP,  1946;  SHEARMAN,  1956b; 

(4Tt)3  R4 


COSGRIFF 


where 


et  al  .  ,  1960;  NIELSON  et  al .  ,  1960), 
P^  =  power  transmitted 

P„  =  power  received 

R  =  range  of  target 
G  -  antenna  gain 

radar  cross-section 


4  T,  t  Power  scattered  per  unit  solid  angle 
Power  incident  per  unit  area 

(  N  IELSON  et  al .  ,  1960)  . 

A  parameter  commonly  used  in  radar  is  C*  ,  the  radar 

cross-section  per  unit  area  of  the  surface. 

The  backscatter  coefficient  Y  (COSGRIFF  et  al . ,  1960) 

is  related  to  0  ,  and  is  the  radar  cross-section  per  unit 
o  v 

area  normal  to  the  direction  of  propagation  at  the  scatter¬ 


ing  area  (  HAGN  ,  196  2)  . 

For  an  element  of  scattering  area,  the  incident  rays 
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are  nearly  parallel,  and  Fig.  2  illustrates  the  relationship 

between  Y  and  0  .  Let 
1  o 

A  =  angle  of  elevation  of  the  rays 

ip  =  -  A  =  angle  of  incidence 

A  =  area  of  scattering  surface 
s 

A  =  A  sinA  -  area  projected  normal  to 
ns 

the  rays  =  "incidence  area" 


The  n 


and  Y 


C 

A 

s 

a 

A 


a  a 

o  s 


a 


or 


A  sin  A  cos  9 

n  T 


P  = 


NIFLSON  et  al  .  (  1960)  and  HAGN  (1962)  use  the  symbol 

■  .  SHEARMAN  (  19  56b)  uses  the  symbol  A  =  0  .  The 
l  o  o 


parameters  Y,  P,  and  A^  are  independent  of  the  radar  par¬ 

ameters  such  as  transmitted  power,  pulse  width,  antenna 
pattern  and  receiver  gain.  They  depend  only  on  the  terrain 
and  the  angle  of  elevation. 


Isotropic  scatter 

If  energy  is  scattered  uniformly  into  the  space  avail¬ 
able  the  scattering  is  called  isotropic.  For  isotropic 
scattering  into  the  hemisphere  above  the  scattering  surface, 
y  is  constant  for  all  angles  of  elevation.  If  all  the  in¬ 
cident  energy  is  scattered  in  this  way,  y  =  2  and  p  =  1. 

To  show  t  his  ,  let 

P  =  power  scattered 

s 

power  incident 

dO  .  element  of  azimuthal  angle 


Rays  falling  on  element  of  scattering  area 


Quasi snecular  scatter  from  suitably 
inclined  surface. 
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d  SI  =  element  of  solid  angle, 
-  d  Ad  ©  . 


From  the  definition  of C  , 

d_Fa 

dTJ 


a  =  4  Tt 


P.  /A 

i  s 


There  fore , 


dP 


_Pi 

4TC 


—~z  -r-6  0. 


4Tt 


y  sin  A  d  Ad  9 


If  energy  is  scattered  uniformly  into  the  volume  above  the 
surface , 

Tt/  2  I  2  Tt 

I  sin  A  d  A  d  0 

O  'o 

“  *  Pi  Y  • 

I  f  P  =  p .  ,  Y  -  2  . 
s  l  1 

Some  workers,  for  lack  of  information  about  y  ,  have 
assumed  isotropic  scatter  (SHEARMAN,  1961;  THOMAS  and 
Me INNES  ,  1962)  ,  and  others,  when  investigating  y,  have 

assumed  it  constant  for  comparison  with  experimental  results 
(CLAPP,  1946;  C05GRIFI'  et  al  .  1960;  NIELSON  et  al  .  1960). 

An  exception  is  SHEARMAN  (1956b),  who  assumed  A^  constant. 


total  P 


_P 
4  Tt 


Y 


Relative  backscatter  coefficient 

The  present  work  concerns  the  relative  backscatter 
coefficient.  To  find  the  absolute  backscatter  coefficient, 
it  is  necessary  to  find  absolute  values  of  the  radar  para¬ 
meters,  but  that  is  to  be  dealt  with  in  another  project. 
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The  relative  backscatter  coefficient  does,  however,  make 
possible  a  comparison  of  backscatter  from  different  terrains 
and  at  different  angles  of  elevation. 

Angular  power  spectrum 

For  radar  echoes,  some  authors  (EVANS  and  PETTENGILL, 

1963;  LYNN  et  al .  ,  1964)  use  the  parameter  P(ip),  the  angular 

power  spectrum,  which  is  defined  as  the  angular  distribution 

of  echo  power  per  unit  area  of  surface.  It  is  the  angular 

part  of  0  . 

o 
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4.2  Laws  of  illumination  an  \  omission 

In  the  formulation  of  echo  power,  use  is  made  of  the 
following  laws,  which  apply  to  a  unit  element  of  surface 
area,  that  is,  for  the  case  when  A^=l,  and  hence 

A  =  sin  A  =  cos  (J 
n  T 

Let  F  =  incident  flux. 


Cosine  law  of  illumination 

Power  incident  on  unit  element  of  surface  =  F  A 

n 

*  F  cos  ip  . 


Lambert's  cosine  law 

The  flux  emitted  from  a  unit  element  of  surface  is  the 
same  in  all  directions. 

Hence,  power  emitted  in  any  direction  from  unit  area 
of  surface  =  (emitted  flux)  . 

=  (constant)  ,  cos  (p 

For  the  echo  power,  the  above  laws  are  combined,  so 
2 

that  P(ip)  =  tor,"  ip  , 


Lomme 1-Se e 1 l nger  law 


For  diffuse  re  1’ lection, 

where  0  =  angle  of  emission, 

hence  P(tp)  =  cos  rp  . 


P(«p,  0  ) 


cos  cp  cos  9 
cos  ip  +  cos  0 


For  backscatter,  0  (p, 


Specular  scatter 


I  f  tp  =  O,  the  above  laws  reduce  to  P(tp)  -  1. 
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For  some  surfaces,  this  may  hold  even  when  ip  is  not 

zero.  <p  is  the  angle  between  the  ray  and  a  line  normal  to 

the  mean  level  of  the  terrain,  and  Fig. 3  shows  that  even  if 

(p  is  not  zero,  the  local  angle  of  incidence  <p  '  on  any 

given  scatterer  may  be  zero.  For  example,  if  a  surface  is 

covered  with  hemispheres,  P((p)  does  not  vary  with  ip  ,  since 

provided  the  radius  a  is  greater  than  0.2 X  ,  0^  is  approx- 

2 

imately  equal  to  n  a  for  each  scatterer.  This  type  of 
specular  scatter  from  suitably  orientated  surfaces  of  suff¬ 
icient  extent  is  referred  to  os  "semispecular "  or  "quasi- 
specular"  scatter. 


The  relationship  of  y  and  Cn 
to  the  scattering  laws 

P  ( tp )  is  the  angular  part  of  0  ,  so 

oc  y  cos  if  , 

If  is  independent  of  ip  ,  P  ( q? )  is  constant,  that  is, 

specular  type. 

If  v  is  independent  of  ip  (i.e.  isotropic  scattering)  , 
P(ip)  oc  cos  ip  ,  that  is,  Lomme  1  -See  1  inge r  scattering. 

Table  1  gives  a  summary  of  the  laws: 
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TABLE  1 

SUMMARY  OF  SCATTERING  LAWS 


Scattering  Law 

P(  tp  )  or  CT 

0 

Oo 

1  COS  (p 

Lambe r  t 

2 

X  COS  <p 

oc  cos  tp 

Lommel-Seelinger 

x.  cos  tp 

OC  constant 

Specular 

constant 

oc  1/cos  tp 

The  use  of  the  scattering 
laws  in  the  present  work 

The  calculations  follow  the  method  of  SHEARMAN  (1956b) , 

with  one  difference.  He  assumed  A  constant:  here  v  is 

o  ' 

assumed  constant. 

In  SHEARMAN,  we  find  the  equivalent  of  dP  C-  A^  dS , 

where  dP  =  power  returned  from  an  element  of  ground  area  dS  , 

and  A  =  0  ,  A  was  taken  as  the  radar  cross-soction  of  n 
o  o  o 

hemispherical  bosses  per  unit  area  of  the  ground,  which  is 

o 

n  Un"  (a  =  radius  of  boss),  independent  of  p 

In  the  present  work,  A^  has  been  replaced  by  y  sin  A  . 
As  a  first  assumption,  y  is  taken  to  be  independent  of  A  . 
This  brings  the  presentation  into  line  with  NIELSON  et  al , 

(  1960)  ,  ilAGN  (1962)  and  COSGR  IFF  et  al.  (1960),  and  there¬ 
fore  the  results  will  be  easily  comparable  wath  theirs. 

The  difference  between  observations  and  calculations  will 
be  a  measure  of  the  deviation  of  the  backscatter  coefficient 
y  from  isot"opic  scattering. 


4.3  Attempts  to  find  Y  theoretically 

The  usual  approach  is  t  '  assume  a  statistical  model 
of  the  surface  and  derive  an  expression  for  y  that  depends 
on  <p  and  the  wavelength  X  ,  and  on  the  statistical  para¬ 
meters  assumed,  such  as  the  distribution  of  heights  of  sur¬ 
face  irregularities  and  the  auto-correlation  coefficient  of 
the  irregularities  across  the  surface.  The  classical  work 
in  this  field  is  by  DAVIES  (1954,  1955) ,  and  most  work  sub¬ 
sequent  to  his  merely  alters  the  statistics.  DAVIES  con¬ 
sidered  a  statistical  model  of  a  surface  with  height  z 
following  a  Gaussian  distribution  with  standard  deviation  C 7  . 
For  radiation  scattered  from  a  "slightly  rough"  surface 
(z  and  0  X)  ,  he  deduced  that  the  radiation  can  be  div¬ 

ided  into  a  coherent  and  an  incoherent  component.  For  the 

2 

latter,  P  ( tp )  cos  tp  ,  which  corresponds  to  Lambert's  Law. 
For  a  "very  rough"  surface  (  CT  Y>  X)  ,  DAVIES  had  some  succ¬ 
ess  in  predicting  the  sea  clutter  of  centimetric  radar  from 
(p  =  0°  to  30°,  His  treatment  failed  to  predict  P((p)  near 
<p  -  T  /  2  ,  which  would  be  of  greater  interest. 

KATZ  and  SPETNER  (1958)  ,  quoted  by  COSGRIFF  et  al . 
(1960)  ,  produced  a  theory  for  angles  near  grazing.  The 
theory  takes  account  of  the  fact  that  there  is  a  tendency 
for  scatter  to  be  specular  in  type,  from  upright  obiects 
such  as  tree  trunks .  Their  result  is  of  the  form 

y  cJO  - i - 

2 

cos  (p  (  A  +  cos  <p  ) 


which,  depending  on  the  value  of  A,  will  be  between  1/  cos  tp 


3 

and  1/  cos  p  .  An  experimental  result  given  by  KATZ  and 
SPETNER  (1960)  is  of  this  type.  A  similar  theoretical  re¬ 
sult  was  derived  by  PEAKE  (1957a,  1957b,  1958) ;  it  applies 

to  a  model  consisting  of  vertical  cylinders,  and  for  p  near 
2 

It/ 2,  A  »  cos  :p  ,  so  Y  os  1/  cos  p  . 

The  usefulness  of  a  statistical  approach  is  doubtful. 
HAGN  (1962)  held  the  view  that  no  simple  mathematical  model 
would  suffice  to  describe  the  actual  physical  process. 
COSGRIFF  et  al  .  (1960)  pointed  out  that  even  with  satisfact¬ 

ory  theoretical  models,  there  would  be  no  unique  way  of  ass¬ 
ociating  a  given  type  of  terrain  with  one  of  the  models. 

The  theories  do,  however,  give  some  idea  of  the  ways  in 
which  y  might  be  expected  to  vary  with  p  .  It  is  to  be 
expected  that  different  types  of  surfaces  and  different 
radio  wavelengths  will  give  rise  to  laws  ranging  from 
Lambert  to  various  forms  of  specular,  and  y  will  corres¬ 
pondingly  depend  on  powers  of  cos  p  ranging  from  1  to  -3. 

4.4  Summary 

The  backscatter  coefficient  is  a  parameter  similar  to 
the  radar  cross-section  per  unit  area,  and  is  related  to  it 
by  the  expression  y  =  0^/cos  P  .  For  isotropic  scattering, 
y  is  constant,  and  in  the  present  work  this  is  assumed  in 
the  calculations,  so  that  the  difference  between  observat¬ 
ions  and  cn’culations  will  be  a  measure  of  the  deviation  of 
y  from  isotropic  scattering.  Theoretical  considerations 

suggest  angular  dependencies  ranging  from  y  Xcos  p  to 
3 

y  c  1/  cos  p. 
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5,  SPECULAR  AND  DIFFUSE  ECHOES 

Specular  and  diffuse  echoes  correspond  to  the  coher¬ 
ent  and  incoherent  components  derived  theoretically  by 
DAVIES  (1954,  1955). 

5,1  Evidence  from  lunar  echoes 

Various  authors  (BROWN,  I960,  HAG FOR S ,  1961,  WINTER, 
1962)  consider  that  the  total  lunar  return  at  any  instant 
of  time  is  composed  of  two  components,  which  are  the  result 
of  two  general  types  of  terrain  features.  The  specular  com¬ 
ponent  arises  from  an  area  where  ip  =0,  or  from  a  large-scale 
irregularity  giving  a  coherent  echo  such  that  cp  is  effect¬ 
ively  zero.  The  diffuse  component  is  predominant  as  ip 

o  2 

approaches  90  ,  and  follows  a  law  such  as  P(  <p)  oc  cos  ip 

( FETTENG ILL ,  I960)  or  cos  t p  (EVANS  and  PETTENGILL,  1963, 

LYNN  et  al . ,  1964) ,  or  some  intermediate  law  such  as  cos 

3  /  2  vp  (BROWN,  i960,  EVANS  and  PETTENGILL,  1963).  The  diff¬ 
use  component  : s  presumably  due  to  a  small-scale  roughness 
superimposed  on  the  large-scale  structure  ( HAGFOKS ,  1961)  . 

The  relative  importance  of  specular  and  diffuse  com¬ 
ponents  depends  on  the  wavelength.  At  optical  wavelengths, 
the  moon  is  a  completely  diffuse  scatterer  ,  Moon  radar  ex¬ 
periments  (LYNN  et  al.,  1964)  show  that  as  the  wavelength 
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is  increased  from  8  mm  to  68  cm,  the  diffuse  component  be¬ 
comes  weaker  relative  to  the  specular  component.  As  might 
be  expected,  the  longer  waves  are  insensitive  to  the  small- 
scale  irregularities.  The  specular  component,  which  at  8  mm 
is  predominant  only  for  tp  close  to  zero,  becomes  predominant 
at  longer  wavelengths  for  values  of  <p  up  to  50°.  While  it 
is  risky  to  extrapolate  this  trend  to  higher  values  of  (p  at 
yet  longer  wavelengths,  it  is  tempting  to  predict  that  at 
16  Mc/s  (A  =  18.75  m)  the  specular  component  of  lunar  echoes 
would  predominate  at  all  angles.  Recent  work  by  DAVIS  and 
ROHLFS  (1964)  at  wavelengths  between  11m  and  22m  confirms 
the  virtual  absence  of  any  diffuse  component, 

5.2  Terrestrial  echoes  at  u.h.f . 

COSGRIFF  et  al ,  (1960) ,  using  wavelengths  near  1  cm, 

found  that  for  surfaces  such  as  grass,  that  appeared  rough 
in  terms  of  \  ,  y  was  fairly  independent  of  incidence  angle 
between  10°  and  80°.  This  corresponds  to  a  predominance  of 
the  diffuse  component.  For  a  smooth  surface  such  as  a  con¬ 
crete  road,  y  decreased  by  as  much  as  25  dB  (COSGRIFF  et  al . 
page  43)  when  A  varied  from  80°  to  10°,  indicating  a  de¬ 
crease  in  the  diffuse  component.  A  longer  wavelength  would 
have  the  same  effect  as  a  smoother  surface  (neglecting 
changes  in  dielectric  constant).  It  might  therefore  be  ex¬ 
pected  that  at  longer  wavelengths,  the  diffuse  component  of 
ground  scatter  will  bo  small,  and  therefore  the  echo  wilJ 
depend  mainly  on  irregularities  contributing  specular  type 
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returns  . 

5.3  Terrestrial  echoes  at  h.f  . 

NIELSON  et  al .  (1960,  page  2)  assumed  specular  echoes 

to  be  negligible  except  near  vertical  incidence,  although  it 
was  noted  that  W.  S.  Ament  believed  tnat  there  is  a  trans¬ 
itional  range  between  specular  (at  vertical  incidence)  and 
diffuse  (at  grazing  incidence)  (BLAKE,  1950) .  The  subse¬ 
quent  report  ( HAGN ,  1962)  admitted  that  a  large  return  could 
come  from  a  combination  of  specular  and  diffuse  components; 
no  attempt  was  made  to  distinguish  between  them,  apart  from 
cases  where  mountains  and  coastlines  gave  obvious  specular 
enhancements.  Such  cases  were  not  used  to  determine  p  , 
and  only  smooth  surfaces  such  as  flat  deserts  and  icefields 
were  studied.  If  we  combine  HAGN '  s  results  for  horizontal 
and  vertical  polarizations,  we  find  that  for  mixed  polariz¬ 
ations  there  is  a  marked  decline  in  p  as  A  varies  from  50° 
to  20°,  indicating  a  weak  diffuse  component.  In  addition, 
a  knee  effect  cuts  off  the  echo  at  still  lower  values  of  A . 

In  the  present  work,  the  knee  effect  is  also  observed. 

Above  the  knee,  however,  y  increases  as  A  decreases  from  30° 

to  10°.  This  increase  cannot  be  due  to  a  diffuse  component 

alone,  since  for  diffuse  scattering  y  should  be  constant 

(  Lomme  1  -See  1  i  nger )  or  vary  as  cos  ip  (Lambert),  but  here  it 

2 

varies  approximately  as  1/  cos"  tp  .  The  conclusion  may  be 
that  the  echoes  are  predominantly  quas i-spec ul or  .  The  diff¬ 
erence  from  HAGN ' s  result  probably  lies  in  the  fact  that  he 


selected  smooth  surfaces  only,  for  which  the  specular  com¬ 
ponent  is  concentrated  near  vertical  incidence;  the  present 
work  involves  very  large  areas  with  a  mixture  of  many  types 
of  surfaces,  the  more  rugged  of  which  could  cause  specular 
type  scattering  even  at  angles  near  grazing  incidence. 

5,4  The  nature  of  the  irregularities 

If  the  irregularities  are  hemispheres  (  SHEARMAIi , 

1956b)  ,  we  should  observe  specular  type  echoes  such  that 

2 

y  or  1/  cos  ip  ,  but  in  fact  we  get  1/  cos  ip  .  It  is  there¬ 
fore  considered  that  the  irregularities  responsible  for  the 
strong  quasi-spccular  component  are  more  rugged  in  shape 
than  hemispheres,  and  cause  comparatively  stronger  back- 
scattering  at  low  angles  of  elevation.  It  has  been  observed 
that  at  low  angles,  vertically  polarized  signals  are  scatt¬ 
ered  back  more  strongly  than  horizontally  polarized  signals 
( HAGN ,  1962) .  In  the  present  work,  with  polarization  ellip¬ 
tical  after  reflection  in  the  ionosphere,  it  seems  likely 
that  the  echo  consists  of  a  large  quasi-spccular  component 
due  to  the  vertical  component  of  the  radiation  at  the  scatt¬ 
ering  sources,  superimposed  on  a  much  weaker,  and  perhaps 
negligible,  diffuse  component. 

The  theory  of  KATZ  and  SPETNER  (Id  58)  and  also  of 
PEAKE  (1967a,  1957b,  lnsg)  considers  scatterers  which  are 

primarily  vertical,  such  as  vertical  cylinders.  Trees  would 
vary  in  shape  from  approximate  hemispheres  to  vertical  cyl¬ 
inders  analogous  to  top-leaded  vertical  antennas.  The  form- 
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er  would  be  expected  to  give  y  k  1/  cos  <p  ,  and  the  latter 
3 

probably  1/  cos  <p  ,  as  it  resembles  the  type  of  antenna  most 

o 

effective  at  low  angles.  An  average  of  1/  cos“cp  is  consist¬ 
ent  with  this  picture.  For  the  sea,  it  is  likely  that  back- 
scatter  arises  primarily  at  the  wave  crests  (KERR,  1951) . 

The  scatterers  would  then  be  conical  in  shape,  that  is,  in¬ 
termediate  between  hemispheres  and  vertical  cylinders,  and 

2 

again  probably  consistent  with  y  1/  cos  <p  . 

In  terms  of  irregularities  such  as  trees  and  wave 
crests,  certain  phenomena  mentioned  elsewhere  in  this  work 
may  be  interpreted: 

(i)  Sea  echoes  are  stronger  than  land  echoes.  In  general, 
the  sea  is  rougher  than  the  land,  and  there  are  more  wave 
crests  effective  in  producing  backseat  ter  than  effective 
trees,  per  unit  area. 

(li)  Sea  echoes  fade  regularly,  but  land  echoes  fade  in  a 
random  fashion,  giving  a  patchy  appearance  on  the  records. 
This  may  be  because  the  wave  crests  tend  to  move  with  a  con¬ 
stant  velocity,  and  therefore  give  rise  to  a  constant  Dopp¬ 
ler  fading  of  the  echoes.  On  the  other  hand,  the  trees 
move  to  and  fro  in  the  wind  in  a  random  fashion,  so  that  the 
Doppler  frequency  is  never  constant.  Fading  ir.  discussed  in 
Section  7.2. 

(lii)  The  knee  in  the  y  (A)  curve  occurs  at  about  13°  for 
land  and  at  lower  angles  for  sea.  Although  re-radiation 
from  a  short  vertical  dipole  is  most  effective  at  low  angles 
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of  elevation,  ground  losses  seriously  diminish  the  radiat¬ 
ion  near  the  grazing  angle.  Ground  losses  are  greater  for 
land  than  for  sea,  hence  the  higher  knee  angle  for  land. 

The  variation  in  the  knee  angle  for  sea  may  be  due  to  var¬ 
iations  in  the  height  of  the  wave  crests  above  the  surround¬ 
ing  surface.  This  is  discussed  in  Section  9.4. 

5 . 5  Summar  y 

Table  2  compares  the  theories  and  experiments  ment¬ 
ioned  in  this  chapter. 

Both  lunar  and  terrestrial  echoes  show  a  trend  away 
from  diffuse  scattering  as  the  wavelength  is  increased  from 
1  cm  to  about  30  m.  The  scattering  irregularities  at  16 
Mc/s  are  therefore  considered  to  be  quasi-specular  scatter- 
ers.  They  are  probably  most  sensitive  to  vertically- 
polarized  waves,  and  are  probably  to  be  identified  with 
trees  and  the  crests  of  sea  waves.  Such  a  picture  is  con¬ 


sistent  with  other  observed  features  of  backscatter  echoes. 
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TABLE  2 

SUMMARY  OF  THEORIES  AND  EXPERIMENTS 


Type 

Y 

Theory 

Experiment 

D  i  f  fuse- 

cos  <p 

L  a  mbe  r  t  ; 
DAVIES , 
"slightly 
rough" , 

«  .h.f  . 

Optical  frequencies. 

1 

Lommel- 
Seelinger . 
"isot  ropic  " 

LYNN  et  al . ,  diff¬ 
use  component  of 
lunar  echoes  at 
u . h  .  f  .  ;  COSGRIFF 
et  al .  ,  rough  terr¬ 
estrial  surfaces  at 
u  .  h  .  f  . 

1/  cos  y 

SHEARMAN , 
hemi spheres. 

HAGN  ,  smooth  sea 
and  land  at  32  Mc/s, 

vert .  poln ,  low 
angles  . 


o 

1/  cos"  q? 


PEAKE  ,  and 
KATZ  and 
SPETM-R  , 
vertical 
cylinders  , 
low  angles. 


Present  work,  random 
sea  and  land  at  16 
Mc/s,  probably  vert, 
poln. ,  low  angles. 


1/ 


3 

cos 


9 


Specular 


» 
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6.  ANTENNA  CALIBRATION 


6.1  Introduction 

In  the  calculation  of  the  power  returned  in  backscatt- 
er  echoes  it  is  necessary  to  know  the  vertical  radiation 
pattern  of  the  antenna  system.  The  antenna  comprised  a  ro¬ 
tating  array  of  four  three-element  Yagi  antennas,  horizont¬ 
ally  polarized.  The  Yagis  were  spaced  1.1  \  apart  and  were 
designed  to  produce  a  narrow  azimuthal  beam.  The  height  of 
the  Yagis  was  0.47  \  above  the  ground,  with  the  intention 
of  achieving  a  broad  lobe  in  the  vertical  plane  with  its 
maximum  at  about  25°  (THOMAS  and  McNICOL,  1960a) . 

The  vertical  radiation  pattern  could  not  be  calcul¬ 
ated  simply,  as  the  ground  was  level  to  only  about  100  feet 
from  the  centre  of  the  array  (about  1,5  \  ) ,  and  also  be¬ 
cause  the  electrical  properties  of  the  ground  were  not 
known.  The  radiation  pattern  in  a  southerly  direction  was 
measured  using  a  transmitter  in  a  helicopter  hovering  at 
selected  heights  at  a  distance  of  2  km  (THOMAS  and  McNICOL, 
1960b).  The  pattern  in  a  northerly  direction  was  measured 
using  a  transmitter  in  an  aeroplane  traversing  a  line  to 
the  north  at  a  series  of  fixed  heights  and  fixed  distances 
(THOMAS  and  McINNES,  1962) .  The  maximum  radiation  was  found 
to  occur  at  an  elevation  of  3  5*  to  the  south,  and  20°  to  the 
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north.  The  discrepancy  between  these  results,  as  well  as 
t ne  need  for  measurements  to  the  east  and  west,  necessitated 
a  more  accurate  method. 

With  a  helicopter  or  aeroplane  the  main  errors  arise 
in  plotting  the  position  of  the  aircraft  during  the  experi¬ 
ment.  The  helicopter  was  not  capable  of  hovering  at  a  con¬ 
stant  altitude.  The  aeroplane  did  not  fly  exactly  on  the 
desired  course.  In  both  cases  theodolite  readings  of  the 
position  of  the  aircraft  could  not  be  taken  at  the  antenna 
site  for  elevation  angles  below  about  20°,  as  the  view  was 
obscured  by  trees  and  haze.  The  aeroplane  was  particularly 
difficult  to  track  with  a  theodolite  owing  to  its  speed. 
These  problems  were  overcome  by  tethering  a  balloon  at  a 
range  of  up  to  5000  feet  (80  \  ) ,  and  operating  a  balloon 

theodolite  at  the  tethering  point.  The  height  of  the  ball¬ 
oon  was  calculated  from  the  amount  of  tethering  line  releas¬ 
ed,  and  when  possible  this  was  checked  by  another  theodolite 
near  the  antenna.  In  calm  weather,  the  balloon  remained  in 
any  desired  position  for  long  periods. 

b  ,  2  Gqu lpme  n  t 

(a)  Balloon .  Meteorological  balloons  of  300,  500  and  800 

g  were  used,  the  most  suitable  being  a  500  g  neoprene  ball¬ 
oon  specially  designed  for  use  as  a  captive  (tethered)  ball¬ 
oon.  The  balloon  was  filled  with  Hydrogen  to  a  diameter  of 
6  or  7  feet,  that  is,  to  a  volume  of  about  350  cubic  feet. 
This  provided  a  free  lift  of  over  10  pounds;  as  the  balloon, 
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transmitter  and  tethering  line  weighed  abc  t  4  pounds,  the 
tension  on  the  tethering  line  was  about  6  pounds.  It  was 
found  from  experience  that  a  nylon  line  of  50  pounds  break¬ 
ing  strain  was  needed,  to  allow  for  wind  gu6ts  and  braking 
during  ascent,  as  well  as  damage  to  the  line  due  to  snags. 

The  line  was  wound  manually  by  means  of  a  heavily-constructed 
reel  one  foot  in  diameter.  Figure  4  illustrates  the  balloon 
and  reel.  In  good  weather  the  balloon  could  be  recovered 
after  its  flight  and  taken  by  truck  to  another  site  for  fur¬ 
ther  use  . 

( b)  Transmit  ter .  The  transmitter  is  illustrated  in  Fig. 5. 
The  box  contained  a  16  Mc/s  crystal-controlled  oscillator, 
modulated  by  a  130  c/s  signal.  Fig.  6  shows  the  miniatur¬ 
ized  construction  of  the  circuit,  and  Fig,  7  is  the  circuit 
diagram.  The  9  volt  battery  shown  was  capable  of  operating 
the  transmitter  for  several  hours.  It  was  found  that  after 
the  first  few  minutes  of  operation,  the  transmitter  output 
decreased  at  the  rate  of  about  8%  per  hour,  and  this  corr¬ 
ection  was  applied  to  the  data  taken  during  the  experiment. 

The  transmitter  was  suspended  25  feet  below  the  ball¬ 
oon  on  a  separate  line,  to  allow  it  to  hang  free  of  the 
tethering  line,  and  to  allow  it  to  rotate  about  a  vertical 
axis  with  the  aid  of  a  swivel.  The  signal  was  radiated  by 
a  short  horizontal  dipole  of  length  1.25  m.  A6  this  is 
much  shorter  than  a  wavelength,  its  radiation  pattern  was 
assumed  to  be  the  same  as  for  an  elementary  dipole  in  free 


1 


Pig.  4.  Balloon  and  reel  used  in  antenna  measurement 


Walloon-,  ornf 


circiu*  cons  true  l  ion. 


Fig.  7.  Circuit  diagram  for  16  V.c/s  balloon-borne 
transmitter  with  130  c/s  modulation. 
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space,  that  is,  radiating  isotropically  in  the  vertical 
plane,  and  having  maximum  azimuthal  radiation  in  the  broad¬ 


side  directions.  Two  celluloid  cones  were  used  as  wind 
vanes  to  keep  the  dipole  rotating,  and  the  result  was  a  kind 
of  turnstile  antenna.  BRUECKMANN  (1963)  had  a  similar  idea 
when  he  proposed  that  an  artificial  earth  satellite  be 
equipped  with  a  turnstile  antenna  for  use  in  calibrating 
large  antenna  systems. 

6 . 3  Recording 

The  signal  was  received  at  the  array  on  the  ground, 
and  the  130  c/s  waveform  was  displayed  on  an  oscilloscope. 
Earphones  were  used  to  aid  identification  and  tuning.  The 
array  was  rotated  continuously,  to  ensure  that  every  read¬ 
ing  was  taken  with  the  maximum  of  the  azimuthal  pattern  of 
the  array  directed  towards  the  balloon.  This  rotation  took 
nearly  4  minutes,  whereas  the  balloon  transmitter  rotated 
about  once  per  second.  Twice  in  each  rotation  of  the  turn¬ 
stile  antenna,  its  maximum  was  directed  towards  the  array, 
and  this  maximum  was  independent  of  the  angle  of  elevation 
of  the  transmitter  with  respect  to  the  receiver,  because  of 
the  uniform  vertical  pattern  of  the  turnstile  antenna.  The 
effect  of  rotating  both  antennas  was  that  the  fast  (2  c/s) 
fading  of  the  transmitted  signal  was  modulated  by  the  slow 
fading  due  to  the  rotation  of  the1  array.  When  the  received 
signal  reached  a  maximum  during  any  one  rotation  of  the 
array,  the  peak  voltage  was  read  from  the  oscilloscope 
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screen. 

To  allow  time  for  rotation  of  the  array  and  for  the 
winding  of  the  reel ,  the  balloon  was  raised  or  lowered  at  5 
minute  intervals.  Each  new  position  corresponded  to  some 
new  angle  of  elevation,  so  new  theodolite  readings  were  tak¬ 
en.  Theodolite  tracking  was  made  possible  at  night  by  att¬ 
aching  a  small  6  volt  bulb  to  the  transmitter,  with  a  separ¬ 
ate  battery  supply.  A  red  globe  was  used  so  that  the  light 
could  be  distinguished  from  the  stars,  and  the  regular 
twinkling  due  to  rotation  of  the  transmitter  also  aided  id¬ 
entification.  Most  of  the  work  had  to  be  done  at  night  to 
satisfy  aviation  regulations,  but  it  was  found  that  at  night 
the  calm  atmospheric  conditions  were  much  superior  to  day¬ 
time  conditions,  and  so  night  flying  was  preferred. 

Communications  between  the  balloon  site  and  the  rot¬ 
ating  array  were  effected  by  means  of  a  radio  link  from  a 
truck  at  the  balloon  site  to  a  hut  near  the  rotating  array, 
and  thence  by  underground  telephone  to  the  control  cabin  of 
the  array.  The  schematic  diagram  of  Fig,  8  shows  this 
ar  rangement . 

6 , 4  Geometry 

The  range  and  elevation  angle  of  the  pilot  transmitt¬ 
er  were  obtained  by  graphical  construction  based  on  theod¬ 
olite  readings  and  the  length  of  the  tethering  line.  The 
length  of  the  line  was  measured  by  counting  the  number  of 
revolutions  of  the  winding  reel.  Although  the  line  was  sub- 
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Pig,  8.  Schematic  diagram  of  method  of  measuring 
vertical  pattern  of  array. 
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ject  to  stretching,  it  was  found  that  after  the  first  use 
of  the  line  no  correction  for  stretching  was  needed,  as  the 
line  was  then  permanently  in  tension  on  the  reel  and  further 
stretching  was  negligible  under  normal  conditions. 

Figure  9  is  a  scale  drawing  illustrating  the  graphical 
construction.  In  the  elevation  at  the  right  of  the  figure, 
the  line  length  and  theodolite  elevation  were  used  to  give 
the  height  of  the  transmitter  and  its  ground  range  from  the 
reel.  This  range  was  then  used  in  conjunction  with  the 
theodolite  azimuth  to  construct  the  plan  and  thereby  find 
the  ground  range  of  the  balloon  from  the  array.  Finally, 
this  ground  range  and  the  balloon  height  were  used  to  con¬ 
struct  the  left-hand  elevation,  from  which  the  range  r  and 
elevation  A  of  the  transmitter  with  respect  to  the  array 
were  read.  The  elevation  with  respect  to  a  second  theodol¬ 
ite  near  the  array  was  also  read  from  the  figure,  and  found 
to  agree  with  readings  from  that  theodolite  to  within  1°. 

In  the  figure,  the  crosses  represent  positions  of  the  trans¬ 
mitter;  the  line  and  balloon  are  sketched  in,  and  A  is  in¬ 
dicated,  for  the  highest  position  only. 

Figure  9  is  taken  from  an  actual  balloon  flight,  with 
the  balloon  tethered  to  the  west  of  the  array  at  a  point 
277°  from  true  north. 

Figure  10  shows  another  flight,  with  the  tethering 
point  44°  from  true  north.  The  dots  represent  positions  of 
the  transmitter,  and  the  numbers  represent  minutes  after 
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the  hour,  the  times  at  which  readings  were  taken  at  the  arr¬ 
ay.  Throe  points  represent  readings  taken  during  descent  of 
the  balloon,  the  remainder  during  ascent.  The  elevations 
show  heights  above  ground  level  at  the  array,  which  was  55 
feet  above  sea  level.  The  azimuth  of  the  balloon  with  res¬ 
pect  to  the  array  varied  by  more  than  40°,  but  this  did  not 
invalidate  the  results,  as  the  ground  in  a  north-easterly 
direction  near  the  array  was  of  fairly  constant  slope,  and 
also  because  the  greatest  deviation  from  44°  was  for  high 
angles  of  elevation.  Provided  A  was  less  than  20°,  the  az¬ 
imuthal  deviation  was  less  than  10°.  For  higher  angles  the 
ground  ray  was  reflected  well  within  the  level  circle  of 
radius  100  feet,  so  the  vertical  pattern  was  considered  to 
be  azimuthally  invariant  above  20°. 


Table  3  summarizes  the  main  balloon  flights. 

TABLE  3 

BALLOON  FLIGHTS 


Dcsignat ion 

Date 

Time  of 
Start 
( hours) 

Azimut  h 

( degree  s ) 

Range 

(ft) 

Max . 
Ht  . 
(ft) 

Max . 

E le v  . 

N  1 

28 .11.62 

0044 

1 

4900 

1200 

22°30 ' 

N  2 

4. 12.62 

2252 

355 

19  40 

7  50 

51° 30 ' 

NE 

5.12.62 

0  209 

44 

1730 

1450 

44°00 ' 

SE 

6  .  2.63 

2248 

143 

1890 

1600 

25°45 ' 

W 

7  .  2.63 

0212 

277 

2800 

1500 

23°00 ' 
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The  aerial  photograph  of  Fig.  11  shows  the  five  sites.  The 
rotating  array  appears  as  a  circle  in  the  centre  of  the  pic¬ 
ture.  The  main  topographical  feature  in  the  vicinity  of  the 
array  is  a  small  creek  that  flows  from  the  south-west  in  a 
clockwise  direction  to  the  north-east,  as  far  as  the  river 
which  is  itself  comparatively  unimportant  in  the  formation 
of  the  radiation  patterns  owing  to  its  greater  range. 

The  ground  profiles  between  the  array  and  the  balloon 
sites  are  plotted  in  Fig.  12.  The  \ertical  scale  is  greatly 
magnified,  and  heights  are  measured  above  sea  level..  The 
proximity  of  the  creek  on  the  western  side  of  the  array  is 
particularly  apparent. 

The  Fresnel  zone  condition  requires  the  range  of  the 

2 

pilot  transmitter  to  be  greater  than  2D  /  \  ,  where  D  is  the 

lateral  dimension  of  the  antenna  (BRUECKMANN,  1955;  KEAY  and 

GRAY,  1964) .  The  horizontal  length  of  the  array  =  3.8  \  , 

2 

so  2D  /  \  =  1728  feet.  For  the  five  balloon  flights,  the 

balloon  was  tethered  at  ranges  greater  than  this,  although 
on  one  occasion,  N2 ,  the  balloon  was  blown  to  half  this 
range  by  the  wind.  The  measured  vertical  pattern  may  there¬ 
fore  be  in  error  for  the  higher  angles  of  elevation  in  this 
instance  . 

6.5  Results 

G(A)  ,  the  antenna  power  gain,  was  obtained  as  follows. 
P(  A )  oc  G  (  A)  /  r2, 
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f'ig.  11.  Aerial  Photograph  of  Radio  Research  Station 
Roggill,  and  environs.  Arrows  indicate  sites  at  which 
the  balloon  was  ;ethered. 
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where  P(A)  =  power  received  at  the  array 

G(A)  =  antenna  power  gain 

r  =  oblique  range  of  pilot  transmitter. 

But  P  0C  S2, 

where  S  =  signal  strength  received  (in  volts) . 

Hence,  G(A)  x  S2  r2. 

For  the  five  flights,  G(A)  was  plotted  against  A  ,  and 
for  comparison  the  five  curves  were  normalized  at  the  maxi¬ 
mum  of  the  main  lobe,  which  occurred  at  about  25°.  The  re¬ 
sults  are  presented  in  Fig.  13,  giving  G  in  dB  relative  to 
the  gain  at  the  maximum  of  the  lobe.  Also  shown  is  a  curve 
obtained  by  measuring  the  pattern  of  a  model  Yagi  at  v.h.f. 

(A.R.P.L.  ,  1956,  Fig.  4 - 55)  .  The  model  was  of  similar 

proportions  to  the  Yugis  on  the  array,  and  was  mounted  0.5\ 
above  a  flat,  perfectly  conducting  surface,  and  its  maximum 
radiation  was  at  28°. 

Only  N2  and  NK  extend  the  measurements  to  high  angles 
above  the  lobe  maximum,  and  their  lack  of  agreement  can  be 
attributed  to  the  declining  accuracy  of  the  measurements  as 
the  balloon  ascends.  As  most  backscatter  echoes  at  16  Me / s 
arrive  at  angles  below  30°,  the  results  are  adequate  for 
the  purpose  of  the  present  project. 

Interest  centres  on  the  disagreements  between  the 
curves  below  the  maximum  of  the  main  lobe.  NE  and  SE  foll¬ 
ow  the  v.h.f.  result  fairly  close’y,  as  the  ground  is  fair¬ 
ly  level  in  those  directions  in  t  e  neighbourhood  of  the 
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array.  N1  agrees  well  with  the  v.h.f,  result  down  to  5°  , 
and  the  diminishing  radiation  below  3°  is  probably  duo  to 
shadowing  by  the  intervening  hill  (See  Fig.  12). 

In  contrast,  N2  and  W  show  enhancement  of  the  pattern 
between  10°  and  4°.  This  is  thought  to  be  due  to  focusing 
of  the  ground-reflected  rays  in  the  concave  mirror  formed 
by  the  creek.  The  western  pattern  will  be  discussed  in 
greater  detail  in  a  later  section. 

6.6  Correction  for  the  Finite 
Range  of  the  Balloon 

Owing  to  the  finite  range  of  the  balloon  from  the  arr¬ 
ay,  it  is  not  always  a  good  approximation  to  assume  that  the 
rays  from  the  pilot  transmitter  arriving  at  the  array  are 
parallel.  If  the  ground  is  flat,  it  is  fairly  easy  to  cal¬ 
culate  the  error  in  the  measured  pattern  due  to  this  ass¬ 
umption,  for  various  balloon  ranges.  This  was  done  for  the 
case  of  a  half-wave  dipole  at  a  height  of  0,5  \  above  the 
ground,  and  assuming  a  ground  reflection  coefficient  of  -1. 
The  correction  terms  obtained  should  be  sufficiently  applic¬ 
able  to  the  vertical  pattern  of  the  array. 

The  geometry  is  shown  in  big.  14.  It  is  seen  that 
tan  t  =  -y  —  tan  A  +  \  /  L, 

Path  difference  between  t  lie  direct  and  ground-reflected  rays 
t  s 

L  L 

d  - —  -  — — 

cos  &  cos  A 


1  (  so  c  ‘ 


s  e  c  A  ) 


54 


=  L  |  s  e  c  tan  (tan  A  +  A  / L )  -  sec  A  J 

Phase  difference  (p  =  2  TC  d  /  A  . 

Let  the  field  at  the  array  due  to  the  direct  ray  =  E^  sin  U)  t 
Then,  for  a  reflection  coefficient  of  -1, 

field  at  the  array  due  to  the  reflected  ray  =  -E^  sin(Wt+<p  )  , 
since  the  assumed  dipoie  is  isotropic  in  the  vertical  plane. 


Total  field  E  =  E  )  sin  tJ  t  -  sin  (  a)  t  + 


<P  )J 


=  2  E^  cos  (  u  t  +  <p  /  2 )  sin  (-  (p  j  2)  , 
i.e.,  the  amplitude  of  the  field  varies  as  -  sin  (p  /  2 , 


E  -C  -  sin  T  d  /  A  , 


sm 


sec  tan  (tan  A  +  A/L)  -  sec  A 


This  amplitude  was  calculated  on  a  computer  for  diff¬ 
erent  values  of  L/\  from  15  to  1000,  and  for  A  varying 
from  1°  to  34°.  The  results  were  plotted  as  a  set  of  curves 
of  E  (in  dB)  against  A  (Fig.  15) .  The  L/ A  =  1000  curve 
was  taken  to  be  approximately  the  situation  for  parallel 
rays,  and  it  was  found  that  the  necessary  correction  to  be 
applied  to  the  experimental  results  was  greatest  at  low  an¬ 
gles  of  elevation,  and  negligible  at  the  maximum  of  the 


lobe .  For  L/ \ 


30,  which  is  the  case  for  measurements 


N2,  NE  and  SE ,  the  correction  was  0.7  dB  at  10  ,  1.5  d3  at 
5°,  and  3  dB  at  2°.  This  correction  was  applied  to  the 
measured  antenna  patterns.  In  cases  where  the  ground  slopes 
away  from  the  array  at  an  angle  of,  say,  2°,  the  correction 
normally  applicable  at  5°,  say,  would  apply  at  about  3°. 


As  angles  below  5  were  not  used  in  calculating  t  lie 


heoretical  antenna  patterns  for  various  balloon  ranges 


strengths  of  backscatter  echoes,  the  correction  is  important 
chiefly  between  5°  and  10°. 

To  compare  the  magnitude  of  this  correction  with  the 
experimental  accuracy,  Fig.  1G  shows  the  experimental  points 
obtained  from  the  balloon  flight  SO  (  L  =  30.8  X  ) ,  as  well 
as  the  v.h.f.  result  referred  to  above.  Also  shown  are  two 
curves  from  Fig.  15  ,  those  calculated  for  L  =  1000  X  and  L 

30  X  .  Above  15°,  the  scatter  of  points  about  a  mean  is 
-  0.5  dB,  which  is  greater  than  the  correction  calculated 
for  L  =  30  X  .  Relow  15°  the  experimental  points,  though 
fewer,  are  probably  more  accurate,  and  the  correction  becom¬ 
es  significant  compared  to  the  experimental  error. 

The  v.h.f.  result  illustrates  the  performance  of  a 
typical  Yagi  under  ideal  conditions.  Even  when  the  experi¬ 
mental  points  are  corrected,  they  differ  from  the  v.h.f. 
result,  because  of  the  uneven  ground  and  the  slightly  diff¬ 
erent  nominal  height  of  the  antennas. 

G  .  7 _ Influence  of  the  Western 

Ground  Profile 

For  the  western  pattern,  the  ground  near  the  array  is 
so  irregular  that  it  was  doubted  whether  a  correction  for 
the  finite  range  of  thf  balloon  as  described  above  would 
suffice.  A  calculation  similar  to  the  above,  but  taking  in¬ 
to  account  the  western  ground  profile,  was  carried  out.  The 
geometry  is  shown  in  Fig.  17.  The  ground  profile  was  spec¬ 
if  ted  in  terms  of  x,  y  and  a  ,  and  the  figure  shows  a  ground- 


Relative  importance  of  range  correction  and 
experimental  accuracy. 
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Ray  geometry  for  v/estern  ground  profile. 
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reflected  ray,  reflected  at  the  point  (  x,  y  )  where  the 
slope  of  the  ground  is  CL  .  Double  reflections  of  the 
ground  ray  were  neglected,  and  the  antenna  pattern  was  cal¬ 
culated  for  a  horizontal  dipole,  \  /2  above  the  ground,  with 

7 

the  ground  range  L  equal  to  2800  feet  and  10  feet,  repre¬ 
senting  the  balloon  flight  W  and  the  case  for  parallel  rays 
respect ively  . 

The  results  are  shown  in  Fig.  18.  Curve  (a)  is  the 
pattern  measured  at  a  range  of  2800  feet,  and  the  dots  den¬ 
ote  actual  points  measured.  Curve  ( b)  was  calculated  for  a 
flat  earth  and  a  range  of  2800  feet,  and  also  applies  for 
the  western  pattern  down  to  15°,  as  the  area  of  the  ground 
where  reflection  occurs  is  flat  for  the  higher  angles. 

Curves  (c)  were  calculated  for  the  western  ground  profile 
and  indicate  the  pattern  below  15°,  for  a  balloon  range  of 
2800  feet.  Curves  (d)  were  calculated  for  the  western  pro- 

7 

file  and  a  range  of  10  feet.  Discontinuities  in  curves  (c) 
and  (d)  were  prominent,  and  where  gaps  were  left  in  the 
patterns  they  may  be  taken  to  be  6  dB,  being  the  contribut¬ 
ion  of  the  direct  ray  only. 

It  appears  from  Fig.  18  that  the  correction  for  the 
finite  range  of  the  balloon  may  be  best  achieved  by  raising 
the  pattern  at  3°  by  about  3°,  similar  to  the  way  in  which 
curves  (c)  are  raised  to  (d) .  The  correction  would  be 
smaller  at  higher  angles. 

The  method  did  not,  however,  give  a  very  accurate 
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18.  Antenna  patterns  for  various  ranges,  west. 

(a^  Pattern  measured  with  "balloon  at  2300  ft. 
lb)  Pattern  calculated  for  2800  ft,  flat  earth. 

(c)  Pattern  calculated  for  2800  ft,  western  profile 
(dj  Pattern  calculated  for  10'  ft,  western  profile. 
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picture  of  the  vertical  radiation  pattern,  even  though  it 
sufficed  to  give  an  idea  of  the  correction  to  be  applied  to 
the  measured  pattern.  It  was  possible  to  improve  the  cal¬ 
culation  of  the  pattern,  provided  parallel  rays  only  were 
considered.  The  improved  calculation  had  four  advantages 
over  the  earlier  method: 

(i)  The  gain  of  the  Yagi  antenna  at  different  angles  in 
free  space  (WILLIAMS,  1950,  page  115)  was  tabulated  in  the 
computer  and  used  to  modify  the  simple  dipole  result. 

(ii)  The  ground  profile  was  taken  at  smaller  intervals  of 
range,  and  the  slope  at  each  point  was  calculated  from  the 
adjacent  heights. 

(iii)  The  focusing  of  the  ground  ray  due  to  the  curvature  of 
the  ground  was  taken  into  account. 

(iv)  Where  more  than  one  reflected  ray  contributed  to  the 
pattern  at  a  given  angle  of  elevation,  all  were  added  vect- 
orially  with  the  direct  ray. 

Fig.  19  shows  the  calculated  pattern  (  Curve  (a)  )  , 

with  the  points  measured  a t  2800  feet  plotted  for  compari¬ 
son.  Above  16°  the  curve  is  due  to  the  sum  of  the  direct 
and  reflected  rays.  Between  16°  and  10°,  only  the  direct 
ray  is  involved.  At  10°  the  reflected  ray  reappears,  in 
such  a  phase  that  it  strongly  enhances  the  pattern,  and  at 
8°  a  second  reflected  ray  gives  further  enhancement.  The 
discontinuities  involve  jumps  as  large  as  6  dB, 

The  discontinuities  do  not,  however,  appear  in  the 
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Fig.  19.  Antenna  patterns,  west. 

Dots  represent  pattern  measured,  with  balloon  at  280C  ft. 
fa)  Pattern  calculated  for  2800  ft,  western  profile. 

(b )  Most  probable  corrected  pattern. 
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measured  pattern.  The  measured  curve  is  fairly  smooth,  and 
even  if  corrected  for  the  finite  range  of  the  balloon,  would 
still  be  a  smooth  curve.  The  fact  that  the  discontinuities 
do  not  appear  in  practice  is  probably  due  to  a  combination 
of  causes: 

( i)  The  ground  profile  is  not  constant  across  the  lateral 
dimension  of  the  antenna.  The  ground  100  feet  on  either 
side  of  the  assumed  profile  line  will  have  a  different  pro¬ 
file,  and  this  is  particularly  important  close  to  the  array 
at  the  edge  of  the  levelled  circle,  where  the  rays  contrib¬ 
uting  to  the  pattern  at  about  16°  are  being  reflected.  The 
actual  pattern  is  therefore  the  average  of  a  number  of  diff¬ 
erent  patterns,  hence  the  observed  smoothness  of  the  measur¬ 
ed  pattern  . 

( li)  Fresnel  diffraction  at  an  edge  causes  appreciable  ill¬ 
umination  in  the  shadow  region.  KNIGHT  et  al ,  (1964)  showed 

that  calculated  antenna  patterns  which  exhibited  sharp  dis¬ 
continuities  could  be  smoothed  by  taking  the  diffracted  en¬ 
ergy  into  consideration  in  the  calculations.  They  were  con¬ 
cerned  with  relatively  simple  profiles,  and  it  is  doubtful 
whether  their  treatment  could  be  applied  to  the  present 
problem,  but  it  does  at  least  justify  the  smoothing  of  a 
calculated  pattern. 

The  curve  ( b)  was  taken  as  the  most  probable  corrected 
pattern  . 

The  corrected  patterns  for  both  east  and  west,  as 


adopted  for  use  in  calculating  relative  backscatter  coeff¬ 
icients,  are  given  in  Fig.  20. 

6.8  Summary 

The  vertical  radiation  patterns  of  the  16  Mc/s  an¬ 
tenna  array  at  Moggill  have  been  measured,  using  a  new 
method  involving  a  balloon-borne  transmitter  with  a  rotat¬ 
ing  antenna.  The  patterns  differ  according  to  the  azimuth¬ 
al  orientation  of  the  array,  and  the  differences  have  been 
interpreted  as  due  to  the  ground  profiles  between  the  array 
and  the  pilot  transmitter. 

For  fairly  smooth  ground,  the  measured  pattern  can  be 
quite  simply  corrected  for  the  finite  range  of  the  balloon. 

For  the  western  ground  profile,  the  theoretical  patt¬ 
ern  of  the  array  is  in  general  agreement  with  the  measured 
pattern,  and  contains  discontinuities  which  can  be  explain¬ 
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Fig.  20.  Corrected  patterns  adopted  for  use  in  calculating 
relative  "backscatter  coefficients. 


ELEVATION  ANGLE  (degrees) 


BACKSCATTER  RECORDS 


7_;_ 

7.1  Production  of  records 

The  various  facilities  of  the  16  Me / s  backscatter 
sounder  have  been  described  by  THOMAS  and  McNICOL  (1962) . 

In  the  present  work,  the  sounder  was  operated  with  rectang¬ 
ular  pulses  of  600  microseconds  duration,  a  pulse  repetition 
frequency  of  25  pulses  per  second,  ani  peak  pulse  power  of 
about  5  kW .  The  antennas  were  aimed  in  a  fixed  direction, 
usually  to  the  east  or  to  the  west . 

The  type  of  data  presentation  was  a  rectangular  range 
—  time  display,  with  an  amplitude  sweep  superimposed.  The 

basic  range - time  display  was  formed  in  a  similar  way  to  a 

type  B  presentation  in  radar  (range - azimuth) ,  but  without 

rotation  of  the  antennas.  Backscatter  echoes  to  a  range  of 
4500  km  were  displayed  on  an  oscilloscope  screen  by  bright¬ 
ness-  or  bl ac kout -modul a t ion  of  the  time-base  trace,  and 
were  recorded  on  film  moving  continuously  past  the  trace  at 
1/2  inch  per  minute.  Such  a  record  enlarged  three  times  is 
illustrated  in  1  ig .  21  (20/12/63,  1530  hours,  192°  mag,  az, , 
no  gain  sweep)  .  The  st:  •  ng  horizontal  lines  are  at  inter¬ 
vals  of  1000  km,  and  the  weaker  ones  at  intervals  of  250  km. 
The  echo  has  its  leading  edge  at  about  1900  km  (p'm^n»  'vhere 
p'  designates  equivalent  fret* -  space  range). 
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Fig.  21.  Backscatter  range — time  record,  20/12/63,  1530 
hours,  192°  mag.  az.,  no  gain  sweep.  Horizontal  scale  is 
about  l-j  inches  per  minute  in  this  enlargement.  Leading 
edge  of  echo  is  at  1900  km. 


Fig.  22.  As  for  Fig.  21,  but  with  sv/ept-gain  unit  operating 
continuously.  Each  3weep  traverses  40  d3  in  1  minute, 
producing  a  range — amplitude  frame. 
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This  range - time  display  was  modified  by  the  use  of 

a  swept-gain  unit,  by  means  of  which  the  receiver  gain  was 
regularly  diminished  by  about  40  dB,  and  then  gradually  re¬ 
stored  to  its  maximum  value  in  a  series  of  11  steps.  Fig.  22 
shows  the  same  film  as  Fig.  21,  a  few  minutes  later,  after 
the  swept-gain  unit  had  been  switched  on.  The  minimum  range 
of  the  echo  is  increasing  at  the  rate  of  about  10  km  per 
minute,  owing  to  slow  changes  in  the  ionosphere,  and  super¬ 
imposed  on  this  variation  is  the  gain  sweep,  traversing  40 
dB  per  minute  every  minute.  Each  stage  of  the  gain  sweep 
lasts  5  seconds,  as  this  was  found  to  be  the  minimum  time 
needed  to  average  out  fading  effects.  The  1  minute  sweep 
is  therefore  a  compromise  between  two  requirements  ,  the  need 
to  record  at  maximum  gain  at  least  once  per  minute,  and  the 
need  to  record  each  gain  stage  for  at  least  5  seconds  durat¬ 
ion  . 

Neglecting  any  variation  in  range  during  one  gain 
sweep,  we  can  regard  each  frame  of  the  swept-gain  records 

as  a  range - amplitude  display,  similar  in  format  to  a  type 

A  presentation  of  the  echo  pulse  with  a  logarithmic  amplit¬ 
ude  scale.  The  equivalence  of  the  swept-gain  and  A-scan 
records  is  illustrated  in  Fig,  23  for  an  idealized  echo 
pulse.  The  advantage  of  this  type  of  record  over  the  A-scan 
lies  in  the  fact  that  in  practice,  blackout-modulation  as 
used  here  appears  to  be  a  more  satisfactory  way  of  averaging 
out  the  echo  pulse  shapes  than  is  possible  with  A-scan. 
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Illustrating  the  equivalence  of  sv/ept-gain 
and  A-scan  records. 


Pig.  23 
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Successive  pulses  are  returned  with  different  shapes  as  a 
result  of  fading,  and  some  averaging  process  is  required. 
NIELSON  et  al .  (1960)  relied  on  film  exposures  long  enough 

to  superimpose  a  large  number  of  fading  shapes  on  the  A-scan 
frame.  This  has  the  disadvantage  that  the  larger  echo  pul¬ 
ses  obscure  the  smaller  ones,  and  what  remains  is  an  envel¬ 
ope  of  the  larger  pulses.  The  present  method  is  thought  to 
provide  a  more  useful  type  of  averaging,  as  the  limits  of 
the  strong  and  weak  echoes  can  be  discerned,  even  if  the 
echo  fades  to  zero,  and  some  estimate  of  the  fading  rate 
can  also  be  obtained. 

In  the  present  work,  the  maximum  amplitude  visible  on 
the  film  was  chosen,  for  the  desired  swept-gain  frame.  Oth¬ 
er  frames  nearby  were  checked  to  see  if  the  reading  was  con¬ 
stant  over  a  few  minutes.  Although  the  echo  appeared  very 
faint  at  the  edges  of  the  swept-gain  films,  different  people 
scaling  the  films  were  able  to  obtain  consistent  results. 

It  may  be  noted  that  the  films  would  be  easier  to  scale  if 
there  were  a  greater  contrast  between  the  recorded  signal 
and  that  which  falls  below  the  recording  threshold  determ¬ 
ined  by  the  film  sensitivity.  A  steeper  threshold  could  be 
achieved  by  making  the  receiver  reject  all  signals  below  a 
given  level  for  each  stage  of  the  swept-gain  unit. 

The  error  in  reading  the  films  was  probably  ^  2  d  B , 

so  the  amplitude - time  curves  were  smoothed,  as  in  Fig. 49, 

and  the  smoothed  values  were  used  in  the  calculations. 
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The  readings  obtained  depend  on  the  depth  of  fading 
(ADAM  and  WHITEHEAD,  1960) .  The  present  work  assumes  that 
the  depth  of  fading  is  fairly  constant  (See  Section  7.2) . 

7 , 2  Fading  Types 

The  example  given  in  Figures21  and  22  show  a  typical 
F-region-propagated  one-hop  backscatter  echo  (designated  IF) 
from  a  direction  of  192°  magnetic.  At  the  range  shown,  the 
backscattering  occurs  on  the  Southern  Ocean  near  Tasmania. 

The  fairly  regular  fading  (about  1  c/s)  noticeable  at  the 
leading  edge  and  the  trailing  edge  of  the  echo  is  fairly 
typical  of  sea  scatter  (SHEARMAN,  1962) .  The  fading  is 
fairly  constant  in  phase  across  the  range  spread  of  the  echo. 

By  contrast,  Fig.  24  (9/4/63,  1800  hours,  260°  mag. 
az.)  is  a  IF  echo  from  the  land  near  Birdsville .  The  fad¬ 
ing  in  general  has  no  preferred  frequency,  and  he  phase  is 
irregular,  giving  rise  to  a  patchy  appearance  of  the  echo. 

It  is  considered  that  this  fading  is  due  to  t ne  motion  of 
the  trees,  which  are  thought  to  be  the  main  scattering  ob¬ 
jects.  CLAPP  (1946,  p.4)  refers  to  the  "wind  in  the  trees" 
observed  by  radar  at  u.h.f.  In  contrast,  SHEARMAN  (1962) 
believed  that  any  fading  of  echoes  from  land  must  be  due  to 
the  ionosphere  alone,  while  fading  of  sea  scatter  represents 
a  combination  of  ionospheric  and  sea  wave  effects.  SHEAR¬ 
MAN'S  opinion  is  probably  correct  for  land  surfaces  bare  of 
vegetation  or  other  moving  objects. 

Fig.  25  (20/12/63,  1200  hours,  192°  mag.  az.)  gives 
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Pig.  24.  Backscatter  grom  land.  Range — time  record, 
9/4/63,  1800  hours,  260°  mag.  az.,  no  gain  sweep. 
Horizontal  scale  is  about  inches  per  minute.  Leading 
edge  of  IP  echo  is  at  1300  km. 


Pig.  25.  Backscatter  from  land  and  sea.  Range — time 
record,  20/12/63,  1200  hours,  192°  mag.  as.,  no  gain  sv/eeo. 
Echo  at  700  km  is  1E3  from  land.  Echo  at  2000  km  is 
IP  from  sea. 
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an  interesting  comparison  of  land  and  sea  backscatter.  The 
echo  from  600  to  1000  km  is  E^-propagated  ground  scatter 
(designated  1E^)  from  the  eastern  highlands  of  New  South 
Wales,  and  shows  the  patchiness  characteristic  of  land 
scatter.  The  echo  at  2000  to  2600  km  is  a  IF  echo  from  the 
Southern  Ocean. 

The  possibility  that  the  regular  fading  attributed  to 
sea  scatter  originates  in  the  ionosphere  is  ruled  out  by  the 
fact  that  when  lE^  and  IF  echoes  are  both  present  in  the 
east,  the  same  kind  of  fading  appears  in  both  echoes  (Fig. 
26,  5/4/63,  1620  hours,  80°  mag,  az.).  Even  if  the  phenom¬ 
enon  causing  the  fading  could  occur  simultaneously  in  both 

the  E  and  F  regions,  the  control  points  for  E  and  F  propa- 

s 

gat  ion  are  at  different  ranges  from  the  sounder.  The  fading 

is  probably  due  to  the  motion  of  the  sea  waves. 

A  similar  fading  record,  obtained  by  direct  back- 

scatter  from  the  sea  at  low  angles,  was  reported  by  DOWDEN 

(1957),  but  his  explanation,  based  on  the  beat  between 

Doppler-shifted  echoes  from  opposite  directions,  will  not 

suffice  here,  as  the  sea  is  viewed  in  only  one  direction. 

Nevertheless,  the  fading  frequency  is  similar  to  the  Doppler 

shift  expected  for  16  Me  /  s  ,  if  sea  waves  of  length  \/  2  are 

cor  si  derod  (CROMBIE,  195  5). 

1  he  velocity  of  a  sea  wave  of  length  L  is  given  by 
F  L 

v  t|—  ’  g  being  the  acceleration  due  to  gravity.  The 

Doppler  shift  A  f  produced  by  sea  waves  of  length  L  =  \  /2 


Fig.  26.  Backscatter  from  aea.  Sweot-gain  record.,  5/4/63, 
1620  hours,  80°  mag.  az.  Duration  of  each  frame  is  1 
minute,  and  fading  rate  is  approximately  48  cycles  per 
minute.  Similar  fading  occurs  in  both  IF  and  1ES  echoes. 


Fig.  27.  Direct  backscattt r  from  field-aligned  E_. 
Swept-gain  record,  20/12/63,  1700  hours,  192  mag.  az. 
Range  of  echo  is  about  400  km,  and  a  IF  sea  echo  appears 
at  2400  km,  showing  possible  magneto-ionic  splitting. 
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travelling  radially  to  the  observer  will  be 

-  f  -  UTT  -0- 

At  16  Me /s,  A  f  =  0.41  c/s. 

In  a  Fourier  analysis  of  the  sea  waves  there  will  be 
two  components  with  wavelength  L  =  \  / 2.  These  will  have 

velocities  +v  and  -v ,  and  so  there  will  arise  two  Doppler- 
shifted  echoes  with  frequencies  f  +  0.41  and  f  -  0.41  c/s. 
The  fading  can  be  attributed  to  the  beat  between  these  two 
signals,  as  follows: 

(a)  The  fading  rate  for  sea  echoes  in  Figures  22,  25  and 
26  is  about  48  per  minute,  or  0.8  c/s.  This  compares  well 
wi t  h  2  A  f  . 


( b)  The  depth  of  fading  may  be  as  great  as  20  dB  between 
maximum  and  minimum.  This  implies  that  there  are  two  sig¬ 
nals  of  approximately  equal  intensity  interfering  with  each 
other  . 


INGALLS  and  STONE  (1957)  measured  the  Doppler  power 

spectrum  of  sea  echoes  at  18  and  24  Mc/s,  and  found  three 

spectrum  peaks.  One  peak  was  at  the  carrier  frequency,  and 

the  others  were  narrow  sidebands  shifted  by  A  f.  The  band- 

_  2 

width  of  the  Doppler-shifted  sidebands  was  about  10  c/s, 
and  the  distribution  of  the  echo  power  among  the  three  peaks 
was  a  function  of  sea  state. 

In  the  present  work,  it  would  seem  that  when  deep  fad¬ 
ing  occurs,  the  sidebands  are  of  equal  power  and  the  carrier 
frequency  component  is  negligible.  This  appears  to  corre- 


late  with  fairly  calm  conditions.  When  the  fading  is  not 


pronounced,  one  of  the  sidebands  is  probably  suppressed  in 
favour  of  the  other.  The  carrier  frequency  component  is 
again  negligible,  as  fading  is  not  observed  at  the  frequency 
A  f  but  only  at  2  A  f. 

It  is  important  to  note  that  the  fading  can  be  explain¬ 
ed  almost  entirely  in  terms  of  the  terrain.  Any  fading  due 
to  ionospheric  effects  is  normally  slow  and  weak,  indicating 
that  the  ionosphere  is  a  comparatively  stable  medium. 

Direct  backscatter  from  field-aligned  irregularities 
in  is  illustrated  in  Fig.  27  (20/12/63,  1700  hours,  192° 
mag.  az.) .  This  occurs  at  a  range  of  300  to  600  km,  and  is 
often  detected  in  a  southerly  direction  when  1E^  echoes  are 
present.  It  usually  shows  fairly  regular  fading,  but  faster 
than  that  of  sea  echoes. 

7.3  Skip  Distance  Focusing 

In  F  region  propagation,  strong  focusing  of  energy 
occurs  at  the  skip  distance,  and  in  IF  backscatter  echoes 
this  occurs  in  both  ionospheric  reflections  (out  and  back) , 
giving  considerable  enhancement  of  echo  power.  The  skip 
ray  returns  to  the  sounder  just  after  the  minimum  time  de¬ 
lay  ray  and  persists  for  a  time  1/2  6  t ,  where  6  t  is  the 
pulse  length.  For  600  microsecond  pulses,  this  corresponds 
to  a  range  interval  of  90  km.  The  echo  power  tends  to  fall 
off  rapidly  beyond  this  range.  The  effect  is  most  pronoun- 
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ced  in  short  range  echoes  (p'  about  1200  km)  and  is  an 

min 

aid  in  distinguishing  between  IE  and  IF  echoes,  as  lE 

s  s 

echoes  show  negligible  focusing.  Some  1E^  echoes,  however, 
are  limited  in  range  by  the  size  of  the  E^  cloud,  and  may 
appear  to  be  focused.  Fig.  28  illustrates  the  effect  of 
skip  distance  focusing,  for  a  IF  echo  (28/3/63,  1556  hours, 
80°  mag.  az.) .  Fig.  29  shows  a  multiple-hop  echo  (23/11/ 
62,  1500  hours,  80°  mag.  az.) ,  when  the  E^  cloud  was  fairly 
large;  focusing  is  not  apparent. 

7.4  Ionospheric  Irregularities 

(a)  Tilts:  The  iso-ionic  contours  in  the  F  region  are  of¬ 
ten  tilted  from  the  horizontal  by  1  or  2  degrees,  causing 
an  increase  or  decrease  in  the  range  of  backscatter  echoes. 
This  effect  is  important  tor  long  range  echoes,  correspond¬ 
ing  to  low  angle  rays  which  are  most  sensitive  to  tilts. 

For  a  quiet  ionosphere,  the  amount  of  tilting  can  be  gauged 
from  the  rate  of  change  of  h'F  measured  at  Brisbane  as  the 
ionosphere  passes  overhead.  The  tilts  can  be  taken  into 
account  in  the  calculations  (See  8.3  (a)  ) ,  but  normally  it 

was  found  preferable  to  select  data  from  days  on  which  tilt¬ 
ing  was  negligible. 

( b)  Ri ppl e s :  If  the  ionosphere  is  disturbed  by  large-scale 
ripples  such  as  travelling  ionospheric  disturbances,  the 
leading  edge  of  the  ripple  becomes  a  second  region  of  focus¬ 
ing  in  addition  to  the  focusing  of  the  skip  ray.  Such 


Pig.  28.  Skip  distance  focus 
ing  in  a  IP  echo.  28/3/63, 
1556  hours,  80°  mag.  az.  One 
minute  sweep. 


Pig.  29.  Multiple-hop  Es 
echoes.  23/11/62,  1500  hours, 
80°  mag.  az.  One  minute 
sweep. 


ripples  may  be  500  km  in  extent  and  may  move  at  about  600 
km  per  hour,  so  that  the  secondary  focused  echo  may  move  at 

twice  this  speed.  Fig.  30  is  a  range - time  record  taken 

at  slow  film  speed  (22/7/SO,  1100  hours,  180°  mag.  az.) , 

The  secondary  focused  echo  moved  towards  the  sounder  until 
it  coalesced  with  the  skip  region  and  temporarily  reduced 
its  range.  A  true  height  analysis  of  vertical  incidence  da¬ 
ta  at  Brisbane  and  Canberra  on  22/7/60  was  performed,  using 
the  method  of  SCHMERLING  and  VENTRICE  (1959) .  Fig.  30a  shows 
the  result.  The  two  disturbances  which  passed  over  Brisbane 
at  1050  and  1210  hours  are  probably  the  same  that  passed  over 
Canberra  at  0930  and  1050  hours  respectively.  Each  disturb¬ 
ance  moved  at  the  rate  of  about  600  km  per  hour,  and  can  pro¬ 
bably  be  identified  with  the  disturbances  on  Fig.  30  which 
reached  their  minimum  recorded  ranges  at  1040  and  1120  hours 
respectively.  In  routine  swept-gain  records  taken  at  i  inch 
per  minute,  the  presence  of  such  ripples  is  not  always  ob¬ 
vious,  and  care  must  be  taken  to  detect  them.  The  secondary 
focusing  appears  or.  swept-gain  records  as  in  Fig.  31  (22/3/ 
63,  1800  hours,  80°  mag.  az.) .  Most  ripples  move  from  south 
or  south-east  to  north  or  north-west,  so  for  records  of  ech¬ 
oes  from  east  or  west,  the  ripple.;  were  usually  observed  end- 
on,  and  their  influence  on  the  range  and  power  of  echoes  was 
not  serious  in  this  case. 

Small-scale  ripples  have  a  focusing  effect  similar  to 
travelling  disturbances,  but  as  a  large  number  are  illumin- 
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ated  at  one  time  by  the  sounder,  their  effect  on  the  back- 
scatter  echoes  is  averaged  out,  and  the  result  ia  "range 
spreading",  similar  to  that  observed  on  vertical  incidence 
ionograms , 

(c)  Sporadic- E :  Eg  clouds  are  usually  classed  as  irregular¬ 

ities.  IE  echoes  at  16  Mc/s  occur  between  500  and  1500  km, 
s  ’ 

the  lower  limit  depending  on  f qE  ,  and  the  upper  limit  being 
determined  by  the  backscatter  coefficient  and  the  amount  of 
range  spreading.  Either  of  these  limits  may  be  modified  by 
the  size  of  the  cloud.  Complete  E^  blanketing  is  rare;  a 
small  cloud  or  one  with  low  f^E^  allows  F-propagated  back¬ 
scatter  to  occur  as  well  as  E^,  provided  IF  propagation  con¬ 
ditions  are  suitable,  but  there  is  always  some  doubt  as  to 
how  much  the  IF  echoes  have  been  attenuated  and  such  cases 
are  best  avoided  in  deriving  the  backscatter  coefficient. 

7 , 5  Summary 

For  the  present  purpose,  swept-gain  records  are  pre¬ 
ferred  to  A-scan.  They  allow  the  echo  power  to  be  read  more 
easily.  They  enable  fading  rates  to  be  recorded,  and  facil¬ 
itate  the  identification  of  propagation  modes  such  as  IF, 

IE  or  direct  backscatter  from  field-aligned  E  .  The  dis- 
s  s 

tinctive  types  of  fading  allow  the  source  of  the  scatter  to 
be  identified,  whether  sea,  land  or  field-aligned  E^.  The 
fading  rates  give  support  to  the  idea  that  backscatter  arises 
from  trees  and  wave  crests. 
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Irregularities  such  as  ripples  in  the  F  region  can  be 
recognized  on  swept-gain  records. 


8.  CALCULATION  OF  BACKSCATTER  ECHO  POWER 


8,1  Introduction 

The  method  of  calculation  is  similar  to  that  given  in 
SHEARMAN  (1956b) .  There  are  three  main  stages. 

( i)  The  Appleton  and  Beynon  equations  (APPLETON  and  BEYNON, 
1940)  are  used  to  calculate  the  equivalent  free-space  path 

p'  and  the  ground  range  D  for  a  single-hop  F-propagated  ray, 
for  given  ionospheric  parameters  and  a  given  angle  of  elev¬ 
ation  A  of  the  ray  at  the  transmitter.  The  equations  ass¬ 
ume  that  the  ionospheric  profile  is  parabolic  and  that  it  is 
constant  with  range  throughout  the  region  traversed  by  the 
rays,  that  is,  tilts  and  ripples  are  absent.  The  earth's 
magnetic  field  is  neglected.  The  ionospheric  parameters  re¬ 
quired  are  the  base  height  h^ ,  the  critical  frequency  f^  and 
the  semi-thickness  y^  of  the  layer  (assumed  parabolic)  at 
the  ionospheric  control  point  of  the  ray.  For  the  skip  ray, 
the  control  point  is  at  range  ^  ,  where  is  the  skip 

distance . 

( ii)  The  focusing  factor  F  is  calculated,  involving  the 
first  and  second  differentials  of  D  with  respect  to  A 
(lii)  The  radar  equation  is  integrated  to  give  the  relative 
echo  power  received  at  the  transmitter. 

The  data  available  include  the  virtual  base  height 
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h  ' ,  the  ordinary  ray  critical  frequency  Cq  and  the  M(3000) 

factor  (scaled  from  vertical  incidence  soundings  at  Brisbane 

every  10  minutes) ,  and  the  16  Mc/s  minimum  equivalent  path 

o'  .  scaled  from  backscatter  records  taken  every  minute. 
r  min  ’ 

The  true  height  hQ  is  normally  about  20  km  less  than 

h  '  for  h  '  about  200  km.  This  was  ascertained  from  true 
o  o 

height  analysis  of  typical  cases.  This  empirical  adjustment 

of  h  1  to  obtain  h  may  be  in  error  by  as  much  as  10  km,  but 
o  o 

this  does  not  affect  the  result  by  more  than  about  0.6  dB 
(See  Section  8,3) ,  and  in  any  case  h^ '  is  only  accurate  to 
-  5  km . 

Two  methods  are  available  to  find  y  .  These  will  be 

m 

discussed  later  (Section  8.3)  . 

8,2  Matching  Vertical  Incidence 

and  Backscatter  Data 

It  is  assumed  that  the  ionosphere  is  quiet,  and  that 
the  ionospheric  configuration  at  any  place  in  the  ionosphere 
remains  constant  as  the  ionosphere  passes  over  Brisbane  from 
east  to  west.  The  ionospheric  parameters  at  Brisbane  at  a 
given  time  will  therefore  apply  at  the  control  point  of  the 
skip  ray  at  an  earlier  time  if  the  backscatter  sounder  is 
pointing  east,  or  at  a  later  time  if  the  sounder  is  pointing 
west.  It  is  required  to  find  the  appropriate  time  lag  be¬ 
tween  Brisbane  and  the  control  point,  and  hence  to  match  the 
given  ionospheric  parameters  at  Brisbane  with  the  appropri¬ 
ate  backscatter  recor  ,  taking  into  account  the  variation  of 
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p'  in  the  backscatter  records ;  and  finally  to  obtain  the 
r  mm 

desired  value  of  p'^  from  the  appropriate  record. 

A  graphical  method  was  devised  to  perform  this  match¬ 
ing  process.  A  set  of  curves  was  prepared,  plotting  calcul¬ 
ated  p'  against  time  lag  for  different  values  of  h  and 

v  min  o 

y  .  These  curves  are  shown  in  Fig.  32,  designated  (a) ,  ( b) 

and  (c) .  They  correspond  to  hQ  =  190,  220  and  250  km,  and 

y^  =  80,  90  and  100  km  respectively.  They  were  constructed 

by  taking  SHEARMAN'S  curves  of  p'  -  D  against  p' 

min  s  min 

(SHEARMAN,  1956b,  Fig  11) ,  and  replotting  them  as  p'  ^ 
against  ,  and  then  changing  the  axis  to  a  time  lag  ax¬ 
is  by  considering  that  at  Brisbane  (27^°  south  latitude)  the 
earth  rotates  at  about  1476  km  per  hour  relative  to  the  sub- 
solar  meridian,  so  that  the  ionosphere  moves  overhead  at 
this  rate  with  respect  to  Brisbane.  Therefore  if  D^=  1476 
km,  for  example,  the  time  lag  to  the  skip  zone  will  be  1 
hour  and  the  time  lag  to  the  control  point  of  the  skip  ray 
will  be  1/2  hour.  The  curves  fa) ,  ( b)  and  (c)  are  referred 

to  here  as  the  overlay. 

The  observed  P'm^n  on  anY  day  was  plotted  against  the 
time  of  day  on  the  same  scale  as  the  overlay,  and  the  inter¬ 
section  of  thus  curve  with  the  appropriate  curve  of  the 
overlay  gave  the  backscatter  frame  appropriate  to  the  given 
ionospheric  parameters,  and  also  the  value  of  P ' m^n  at  that 
time. 

For  example,  Fig.  32  shows  the  time  lag  curves  super- 
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Pig.  32.  Matching  vertical  incidence  and  backscatter  data 

(a)  h  =  190  km,  ym  =  80  km. 

v  '  o  m 

(b)  h  =  220  km,  ym  =  90  km. 

(c)  hQ  =  250  km,  y^  =100  km. 

(d)  p*  for  3/5/63,  80°  mag.  az. 


imposed  on  part  of  the  p'^  curve  for  3/5/63  (80°  mag,  az.) 
(Curve  (d) ) .  Suppose  it  is  desired  to  find  the  swept-gain 
frame  corresponding  to  a  control  point  having  vertical  in¬ 
cidence  parameters  measured  at  Brisbane  at  1700  hours.  The 
overlay  is  placed  with  its  zero  time  lag  over  1700  hours  on 
curve  (d) .  The  height  h  '  at  Brisbane  at  1700  hours  was  220 

km,  so  h  was  taken  to  be  200  km,  and  curve  (a)  was  chosen, 
o 

The  intersection  of  curves  (a)  and  (d)  occurred  at  a  time 

delay  of  39-J-  minutes,  when  o'  was  207  5  km.  The  time  of 

mm 

the  backscatter  frame  required  was  1700  -  39-£  or  1620^, 
which  could  be  read  directly  from  the  time  scale  of  curve 

(d)  .  Similarly,  the  vertical  incidence  data  at  1710  hours 
at  Brisbane  (say)  coaid  be  matched  with  another  swept-gain 
frame. 

It  is  to  be  noted  that  curve  (d),  representing  east¬ 
ern  backscatter,  is  plotted  with  time  progressing  from  right 
to  left.  For  western  backscatter  echoes,  the  time  is  plott¬ 
ed  from  left  to  right,  as  in  Fig.  33  for  14/4/63  (260°  mag. 
az.)  (Curve  (e)).  At  1700  hours,  h^  was  about  200  km,  so 
when  the  zero  of  the  overlay  was  placed  over  1700  hours, 
the  result  was  given  by  the  intersection  of  curves  (a)  and 

(e)  .  The  time  lag  was  32  minutes,  p'  was  1730  km,  and 

the  appropriate  backscatter  frame  was  1732  hours. 

If  the  estimated  value  of  h  is  in  error  by  10  km  (as 

o 

is  likely)  ,  the  value  of  p'  derived  will  be  in  error  by 

’  mm 

about  10  km,  and  the  backscatter  swept-gain  frame  chosen  may 
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Pig.  33.  As  for  Pig.  32,  except: 

(e)  P'nin  I",°1’  1 4/4/63,  260°  mag.  az. 
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be  1/2  minute  in  error,  which  is  less  than  the  probable  err¬ 
or  in  recording  and  scaling  the  backscatter  echoes. 

8.3  Estimation  of  y 

rn 

The  next  step  in  obtaining  data  for  calculation  was  to 

find  y  .  Two  methods  were  used,  and  the  second  was  found  to 

be  more  satisfactory  than  the  first. 

( a )  First  Method,  using  M(3QQQ) 

TAIEB  (19G2)  published  a  set  of  curves  parametric  in 

h^  ,  relating  y^  to  M{3000).  His  simple  theory  neglects  the 

earth's  magnetic  field,  and  his  more  developed  theory  has 

not  been  produced  in  graphs  except  for  a  few  examples  of 

particular  cases.  To  see  if  his  method  could  be  profitably 

employed  in  the  present  work,  values  of  h  '  and  M(3000) 

scaled  from  vertical  incidence  soundings  at  Brisbane  were 

used  in  conjunction  with  the  simple  theory,  neglecting  the 

Earth's  field  and  the  correction  for  hQ ,  and  corresponding 

values  of  y^  were  read  from  the  curves.  The  Appleton  and 

Deynon  equations  were  used  to  calculate  p'  for  a  range  of 

values  of  A  ,  and  p'  calculated  in  this  way  was  found 

to  disagree  with  the  observed  value  of  p'  recorded  on 

min 

the  corresponding  backscatter  record.  A  discrepancy  of  20% 
or  more  was  common,  and  the  average  calculated  p'  was 
systematically  about  10%  too  large.  Several  possible  ex¬ 
planations  were  considered. 

(i)  The  true  height  h^  should  have  been  used  instead  of 


h  *  .  But  when  this  was  done,  TAIEB's  curves  yielded  higher 
values  of  which  compensated  for  the  lower  values  of  h  , 

and  the  resulting  calculated  p'  was  still  in  error  by  the 

same  amount  , 

(ii)  The  earth's  field  should  have  been  taken  into  account. 

In  a  few  cases  when  this  was  considered,  y^  was  slightly  re¬ 
duced,  thereby  reducing  the  calculated  p'  and  the  system¬ 

atic  error,  but  the  random  discrepancies  were  still  quite 
large  . 

(iii)  The  ionospheric  parameters  were  assumed  to  be  constant 

over  a  large  area  of  the  ionosphere  from  range  v  to  in¬ 

finity,  or  at  least  to  the  range  of  the  control  point  of  the 
tangent  ray.  This  means  that  any  tilts  in  the  iso-ionic 
contours  were  neglected,  and  this  neglect  could  have  produc¬ 
ed  the  discrepancies,  A  computer  program  due  to  B.  A. 

Me  Innes  was  used,  taking  into  account  the  ionospheric  par¬ 
ameters  at  every  200  km  along  the  propagation  path.  The  re¬ 
sults  for  <14  cases  are  shown  in  Pig.  34,  where  the  calculat¬ 
ed  values  of  p'  ^  are  plotted  against  the  observed  values. 
The  full  line  represents  the  agreement  desired,  and  the  bro¬ 
ken  lines  indicate  -  10%  discrepancy.  The  scatter  of  points 
is  still  serious,  and  there  is  still  a  systematic  disagree¬ 
ment  of  10%. 

( iv)  The  M( 3000)  values  may  be  in  error.  SHEARMAN  (1956a) 
found  a  systematic  difference  of  10-15%  between  observed 

and  calculated  values  of  P'mjn  w'nch  agrees  with  the  present 
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findings.  Ho  was  inclined  to  attribute  this  discrepancy  to 
the  presence  of  echoes  from  azimuths  other  than  the  centre 
of  the  main  antenna  lobe,  but  such  an  explanation  will  not 
suffice  for  the  narrow  beam  in  use  in  the  present  work.  In 
the  discussion  following  SHEARMAN'S  paper  it  was  suggested 
by  W.  R.  Piggott  that  all  conventional  ionosondes  give  val¬ 
ues  of  h’  which  are  systematically  too  high  by  5-15  km,  and 
this  introduces  a  significant  consistent  error  in  the  M  fac¬ 
tor  deduced  . 

(v)  The  extraordinary  ray  should  be  considered.  The  ex¬ 
traordinary  ray  critical  frequency  for  vertical  incidence 
is  larger  than  f ^  by  an  amount  equal  to  J/2  fjj,  where  f^  is 
the  gyrofrequency  at  the  particular  locality  considered 
(provided  f ^  U  fQ).  At  Brisbane,  1/2  fH  is  approximately 
0.7  Mc/s,  and  when  fQ  +  1/2  f^  was  used  in  place  of  fQ  in 
the  calculations,  the  calculated  value  of  p'  was  reduced, 

and  the  average  systematic  error  was  virtually  eliminated. 

However,  this  approach  is  not  valid  at  oblique  incidence, 
and  in  fact  for  east-west  propagation  the  difference  between 
the  ordinary  and  extraordinary  ray  critical  frequencies  is 
negligible,  amounting  to  only  about  1/20  of  the  difference 
for  north-south  propagation. 

The  most  likely  sources  of  error  are  in  explanations 
(ii)  and  ( iv) ,  As  it  is  not  easy  to  apply  TAIEB's  method 
when  the  earth's  field  is  considered,  and  as  there  seems  to 
bo  some  doubt  aoout  the  accuracy  of  M(3000) ,  and  as  in  any 
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case  there  would  still  be  large  random  errors  even  if  the 

systematic  error  could  be  eliminated,  it  seemed  preferable 

to  seek  an  alternative  method  of  estimating  y  . 

m 

( b)  Second  Method,  using  observed  p' 

_ ' _ min 

As  explained  in  Section  8.2,  for  every  set  of  values 
of  h^  and  f  obtained  at  vertical  incidence,  the  correspond¬ 
ing  value  of  p1  can  be  obtained  from  backscatter  records, 
mm 

The  parameter  y^  can  be  obtained  by  calculating  back  from 
P  * 

K  min 

With  the  aid  of  a  computer,  p'  was  calculated  for 

K  K  min 

a  wide  range  of  values  of  h  ,  f  and  y  ,  and  the  results  are 
plotted  in  Figures  35  to  40.  Each  chart  applies  for  a  sing¬ 
le  value  of  h  ,  and  f  is  plotted  against  p'  for  a  range 

o’  o  v  min 

of  y^  values,  giving  a  set  of  curves  parametric  in  y^ .  y 
was  obtained  from  the  observed  p'  by  means  of  these 

charts,  and  then  used  in  subsequent  calculations. 

In  this  procedure  it  is  assumed  that  the  ionosphere 
is  not  tilted,  so  care  was  taken  to  reject  data  in  which  hQ 
was  varying  rapidly.  However,  a  slight  tilt  of  about  0.5° 
would  give  rise  to  an  error  of  1°  in  A  the  elevation 

angle  of  the  skip  ray,  and  an  error  of  this  magnitude  is 
acceptable.  Cn  the  other  hand  an  error  in  p’  such  as 

commonly  arises  when  the  first  method  above  is  used,  could 
have  a  large  effect  on  the  calculated  power  returned. 

The  second  method  also  tends  to  cancel  out  any  errors 


96 


PiC.  35.  Chart  for  deducing  y  . 
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Pig.  37.  Chart  for  deducing  y  . 
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Fig.  38.  Chart  for  deducing  y  . 
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Fig.  39.  Chart  for  deducing  y^. 
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Chart  for  deducing  y  . 
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Fig.  40 


in  h  .  If  h  is  10  km  too  small,  then  y  road  from  charts 

will  be  typically  7  km  too  large,  in  order  to  satisfy  the 

correct  P ' m ^ u •  In  consequence,  the  computed  value  of  A 

will  be  0.25°  too  small,  and  P  ,  the  calculated  peak  power, 

will  be  0.6  d3  too  large.  All  these  errors  are  tolerable. 

Fig.  41  illustrates  the  consequences  of  errors  in  h^  for  a 

typical  case  in  which  f  =7  Me / s  and  p*  .  =  1500  km.  P 

J  F  o  min  c 

is  given  in  dB  relative  to  its  value  at  h^=  200  km,  for  con¬ 
stant  G  and  y  ,  which  in  any  case  do  not  vary  much  over  the 
range  of  A  represented. 

In  practice,  the  accepted  variation  in  h^  was  5  km  in 
10  minutes,  that  is,  in  a  range  interval  of  245  km  on  the 
east-west  ionospheric  height  contours.  Even  10  km  was  acc¬ 
epted  if  it  was  an  isolated  variation.  A  tilt  of  this  mag¬ 
nitude  could  be  due  to  a  travelling  ionospheric  disturbance, 
ns  illustrated  in  Fig.  30a.  But  such  a  disturbance  advanc¬ 
ing  from  the  south  would  not  have  much  effect  on  east-west 
propagation,  so  the  10  km  variation  in  h^  can  be  regarded 
as  an  error  as  far  as  the  east-west  contours  are  concerned. 
This  error  would  be  partially  offset  by  the  method  of  cal¬ 
culating  y  and  would  be  virtually  eliminated  by  averaging 
adjacent  points  in  the  curve  of  y  against  A  .  An  estimate 
of  the  actual  magnitude  of  east-west  tilts  may  be  gained 
from  the  result  that  the  knee  angle  for  land  is  fairly  con¬ 
stant  (13°  -  1°) ,  which  suggests  that  tilts  were  effective¬ 


ly  about 


Pig.  41.  Effect  of  errors  in  h  on  y  ,  a  and  P  .  when 

o  m  s  c 

fQ=  7  Mc/s,  P'min=  1600  km,  and  0  end  ^  are  constant. 
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In  later  versions  of  the  main  computer  program,  the 
use  of  the  charts  to  find  y  was  obviated  by  computing 
from  p'  by  successive  approximations,  before  proceeding 

to  find  Pc  and  A  .  But  in  general  the  charts  were  pre¬ 
ferred,  as  they  enabled  one  to  recognize  and  eliminate  ab¬ 
surd  results,  such  as  negative  y  ,  more  effectively,  and 
also  because  the  use  of  the  charts  involved  a  large  saving 
in  computer  time. 

A  third  method  of  finding  y  would  be  by  calculation 

m 

from  p’f  curves  by  the  method  of  SCHMERLING  and  VENTRICE 
(1959).  But  the  true  height  of  the  maximum  ionization  ob¬ 
tained  in  this  way  is  subject  to  large  errors,  as  the  ruled 
overlay  must  be  placed  in  line  with  the  f  value,  and  any 
error  in  judging  this  is  greatly  magnified  by  the  concen¬ 
tration  of  rulings  on  the  overlay  near  this  frequency,  that 
is,  the  concentration  of  h’  values  to  be  read  from  the  curve 

where  h'  is  large  and  rapidly  varying.  Probably  y  could 

m 

bo  obtained  to  within  about  -  10  km,  but  this  would  lead  to 
large  errors  in  the  calculation  of  p' 

min 


8.4  Calculation  of  P  and  A 
_ c _ s 

The  complete  formulas  are  given  in  Appendix  II.  The 
method  is  summarized  here. 

The  input  data  were  h  ,  f  and  y  ,  and  the  antenna 
gain  G,  tabulated  for  every  degree  of  A  from  5°  to  40°. 

The  Appleton  and  Beynon  equations  were  stored  in  a  subrout- 
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ine  for  very  frequent  use;  also  in  this  subroutine  was  the 

formula  to  calculate  dD/d  A  . 

The  skip  distance  D  and  the  skip  angle  A  were  found 

s  s 

by  calculating  the  penetration  angle  and  then  decreasing  A 
in  steps  until  dD/d  A  was  zero,  when  the  corresponding  val¬ 
ue  of  D  was  D  ,  and  the  value  of  A  was  A  .  The  second 
s  ’  s 

2  2 

differential  d  D/ d  A  at  D^  was  then  calculated  numerically 
from  the  gradient  of  the  dD/d  A  vs  A  curve,  and  was  used  to 
calculate  the  focusing  factor  at  the  skip  distance.  The  el¬ 
evation  angle  A  was  decreased  in  further  steps,  and  at  each 
stage  the  first  differential  dD/d  A  was  calculated  and  used 
to  find  the  focusing  factor. 

Let  b  D  be  an  interval  of  ground  range  from  which  the 
echo  is  being  simultaneously  received. 

Then  the  actual  power  received  back, 


II)  . 


G  1  y 


sin  A  D 


dD  (See  Appendix 


The  interval  6  D  depends  on  the  pulse  width  6  t. 

If  focusing  is  neglected, 

6D  =  ^c-  ft  £  =  — km  for  6 1  =  600  microseconds. 
cosA  cosA 

For  a  thick  ionosphere,  b  D  is  considerably  increased  near 
D  ,  an  effect  known  as  time  focusing. 

The  data  was  chosen  for  the  control  point  of  the  skip 
ray,  because  it  is  desired  to  calculate  the  peak  value  of 
the  echo  power,  and  owing  to  skip  distance  focusing  this 
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power  travels  via  the  skip  ray  and  adjacent  ray9,  provided 

the  antenna  gain  is  reasonably  constant  over  the  range  of 

angles  involved.  To  illustrate  this  focusing  for  a  typical 

case  where  h  =  200  km,  f  =  7  Mc/s  and  y  =50  km,  the  var- 

iation  of  F  with  D  was  calculated  and  plotted  in  Fig.  42. 

The  elevation  angles  corresponding  to  various  points  on  the 

F - D  curve  are  indicated  on  the  graph.  At  D  ,  F  may  be  ten 

3 

times  its  value  at  D  +  6  D,  and  therefore  the  echo,  which 

s 

2 

depends  on  F  ,  may  be  enhanced  a  hundred  times  at  .  At 

long  ranges,  F  increases  again  (horizon  focusing) ,  but  at 

the  ranges  and  elevation  angles  concerned  this  effect  is  not 

important  in  the  present  work. 

The  value  of  F  at  D  has  to  be  calculated  by  a  more 

s 

2  2 

rigorous  method  than  at  longer  ranges,  and  involves  d  D/dA  , 
as  mentioned  before.  This  is  because  the  ordinary  formula 
is  an  approximation  which  results  in  F  going  to  infinity  as 
D  approaches  .  The  curve  of  Fig,  42  shows  that  there  is 
a  smooth  transition  between  the  rigorous  and  approximate 
formulas  as  applied  in  this  program. 

The  value  of  6  D  in  this  example  is  derived  in  the 
course  of  the  calculations,  and  turns  out  to  be  130  km  for 
a  600  microsecond  pulse  and  the  given  ionospheric  configur¬ 
ation  . 

The  calculated  relative  peak  power,  taking  y  as  con¬ 
stant  ,  i s  given  by  D  +  6  D 

P  =  ;  S  G2F2sin  A  D  dD. 

c 

D 
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The  integration  was  performed  by  calculating  the  area  of  the 
2  2 

strip  G  F  sin  A  D  dD  each  time  A  was  incremented  until 

such  time  as  p'=  p'  +  Ic  (that  is,  when  D  =  D  +  6  D)  . 

The  sum  of  the  strips  then  gave  the  total  integral  for 

subsequent  comparison  with  the  observed  P^  . 

The  optimum  size  of  the  steps  d  A  was  chosen  by  trial 

and  error  so  that  P^  varied  by  less  than  0.1  dB  from  its 

value  when  d A  was  very  small.  The  values  finally  chosen 

for  dA  were  0.2°  for  A  between  the  penetration  angle  and 

A  ,  and  also  until  p'  reached  p'  ,  and  then  0.5°  for  p* 
s’  min ’  K 

between  p'  and  p1  .  +  k  6  t. 

min  min 

It  might  be  thought  that  it  would  be  profitable  to 
calculate  the  complete  P^ - p'  echo  pulse  shape  for  compar¬ 
ison  with  range - amplitude  (swept-gain)  records,  to  give  a 

y - A  curve  for  a  wide  range  of  A  and  not  -just  for  A 

'  s 

This  was  attempted  for  1E^  echoes,  assuming  that  the  layer 
was  at  a  height  of  110  km  and  of  negligible  thickness,  but 
results  so  obtained  for  the  backscntter  coefficient  were 
very  inconsistent.  This  was  attributed  to  the  fact  that  an 
cloud  does  not  usually  cover  the  whole  area  of  the  E  reg¬ 
ion  illuminated  by  the  sounder,  and  the  echo  varies  in  shape 
as  the  cloud  passes  through  the  illuminated  area.  Ideally 
it  should  be  possible  to  select  IE  echoes  at  times  of  com¬ 
plete  E^  blanketing.  The  leading  edge  of  the  echo  then  de¬ 
pends  only  on  the  critical  frequency  (assumed  constant 

throughout  the  cloud),  and  the  t i ailing  edge  should  provide 
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data  useful  for  deriving  y  ,  but  unfortunately  the  tail  of 
the  echo  is  often  lengthened  due  to  range  spreading  due  to 
small  irregularities  in  the  cloud,  which  can  no  longer  be 
considered  to  act  as  a  smooth  mirror.  For  oblique  reflect¬ 
ion,  the  irregularities  have  their  greatest  influence  at  low 
elevation  angles  (long  ranges)  ,  although  it  might  be  poss¬ 
ible  to  overlook  them  at  the  leading  edge.  These  considerat¬ 
ions  make  it  very  hazardous  to  draw  conclusions  about  the 
backscatter  coefficient  from  the  shape  of  the  complete  range 
- amplitude  records. 

For  IF  records,  what  is  said  about  range  spreading  of 
1E^  echoes  applies  even  more,  although  the  F  layer  can  norm¬ 
ally  be  assumed  to  extend  across  the  entire  illuminated  F 
region.  The  F  region  irregularities  that  cause  spreading 
can  be  much  larger  than  E^  irregularities.  Not  only  may  the 
irregularities  cause  confusion  at  long  ranges,  but  also  the 
tail  of  the  IF  echo  is  reflected  in  the  ionosphere  at  a  much 

greater  range  than  IE  ,  and  the  Brisbane  vertical  incidence 

s 

data  appropriate  to  the  control  point  of  the  skip  ray  cannot 
confidently  be  assumed  at  the  longer  ranges.  Therefore  the 
ordinary  ionospheric  parameters  are  known  with  less  certain¬ 
ty,  and  there  is  no  way  of  checking  the  assumed  parameters 
such  as  is  afforded  by  the  observations  of  p'  at  the 

leading  edge. 

It  seems  preferable,  therefore,  to  concentrate  on  the 
peak  observed  value  of  ,  and  to  observe  its  .c.iations 
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with  the  time  of  day,  as  h  ,  f  and  y  vary,  and  hence  as 
A  vanes.  By  plotting  PR  -  against  As  it  should  be 
possible  to  obtain  the  relative  backscatter  coefficient  at 
the  angles  represented. 

Every  10  minutes  h  and  f  were  available  from  vert- 

o  o 

ical  incidence  soundings,  and  when  these  were  matched  with 
backscatter  records,  it  was  possible  to  draw  up  a  table  giv¬ 
ing  time  of  vertical  incidence  sounding,  h  ,  f  ,  correspond¬ 
ing  backscatter  record  time,  p'  ,  Pn  and  y  .  Calculat- 
ions  then  gave  P^  and  Ag,  hence  PR  -  P^  could  be  plotted 
against  A  Table  4  is  an  actual  example  (14/4/63,  260° 

mag.  az.) .  The  computer  time  involved  was  about  3  minutes. 

Pp ,  the  peak  power  scaled  from  the  backscatter  records, 
is  given  in  dB  relative  to  the  highest  gain  stage  of  the 
swept-gain  unit.  Thus  if  an  echo  is  just  discernable  on  the 
swept-gain  frame,  PR  =  0  dB,  and  if  it  extends  right  across 
the  frame,  PR  is  about  40  dB.  P^.  is  relative  to  an  equally 
arbitrary  zero,  but  PR  -  P^.  suffices  to  describe  the 
re 1 at ive  backscatter  coefficient. 

The  calibration  of  the  swept-gain  unit  appropriate  to 
the  period  22/10/52  to  16/5/63  is  shown  in  Fig.  43.  Gain 
stage  12  represents  the  lowest  receiver  gain,  and  therefore 
echoes  discernable  in  section  12  of  the  swept-gain  frame 
are  the  strongest.  The  errors  shown  on  the  calibration 
chart  .are  cumulative  errors  relative  to  stage  12.  Compar¬ 
ison  wit1-  the  dotted  line  shows  that  the  gain  steps  are 
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TABLE  4 


EXAMPLE  OF  CALCULATIONS 
14/4/63,  260°  mag,  az. 


TIME 
(  hr  ) 

h 

0 

(  km) 

f 

o 

(Me) 

TIME 

(  hr) 

P  ' 

min 
(  km) 

PR 

(dB) 

ym 
(  km) 

pc 

(dB) 

A  g 
(°) 

P  -P 

R  C 

(dB) 

1500 

2  20 

8.4 

1523 

1370 

31.0 

61 

-127.3 

24.6 

158.3 

1510 

210 

8 . 2 

1533 

1360 

31.0 

61 

-127  .7 

23  .7 

158  .7 

1520 

210 

8 . 2 

1543 

1390 

31.0 

66 

-128.0 

23.4 

159.0 

1530 

210 

8.2 

1553 

1370 

31.0 

63 

-127.4 

23.6 

158 .4 

1540 

210 

7.8 

1603 

1370 

30.6 

50 

-129.2 

22.4 

159.8 

1550 

210 

7  .7 

1615 

1450 

30.3 

59 

-129,8 

21.4 

160.1 

1600 

200 

7.4 

1626 

1500 

29.6 

49 

-130,2 

20 .7 

159.8 

1610 

200 

6.7 

1636 

1500 

29.3 

40 

-134.3 

17 .9 

163.6 

1620 

200 

6.7 

1647 

1550 

29  .0 

46 

-134,8 

17 .6 

163 .8 

1630 

200 

6.4 

1657 

1550 

28,6 

36 

-136.1 

16,7 

164.7 

1640 

200 

6.5 

1708 

1600 

28  .0 

44 

-135.7 

16 .7 

163.7 

1650 

200 

6.4 

1719 

1650 

27 .3 

46 

-136.6 

16 . 1 

163.9 

1700 

200 

6,3 

17  30 

1700 

26 .3 

47 

-137 . 4 

15 . 5 

163.7 

17  3  0 

200 

6,3 

17  44 

1900 

24.4 

67 

-137 .4 

14.5 

161.8 

17  20 

200 

6 , 3 

1755 

1950 

12  .0 

72 

-137 , 8 

14 . 2 

149.8 

17  30 

190 

5.9 

1810 

2200 

3  .0 

82 

-138.8 

11.8 

141 . 8 

17  40 

190 

5.4 

1822 

2300 

3  .0 

65 

-141.0 

10.0 

144.0 

approximately  3  dB  each,  except  Por  steps  1  to  2  and  2  to 
3,  which  in  any  case  could  not  be  calibrated  as  accurately 
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as  the  others,  owing  to  receiver  noise.  The  peak  echo  power 
is  normally  recorded  in  stages  3  to  12,  thereby  minimizing 
receiver  noise.  (Note  that  this  would  not  be  so  in  A-scan 
records  taken  with  maximum  receiver  gain)  . 

8,5  Selection  of  Data 

The  following  principles  governed  the  selection  of 

data. 

(a)  Both  backscatter  and  vertical  incidence  data  should  be 
simultaneously  available  and  clearly  recorded.  This  condit¬ 
ion  was  sometimes  not  fulfilled,  owing  to  breakdowns  in  eit¬ 
her  sounder  or  use  of  the  backscatter  sounder  for  other  work. 

( b)  Radio  interference  should  be  at  a  minimum.  Interfer¬ 
ence  sometimes  rendered  the  vertical  incidence  records  use¬ 
less,  even  when  good  backscatter  records  were  available. 

(c)  No  large  ripples  should  be  present  in  the  ionosphere. 
Large  ionospheric  disturbances  are  most  frequent  in  the  win¬ 
ter,  therefore  other  seasons  were  preferred. 

(d)  Sporadic-E  should  not  be  present.  is  most  preval¬ 

ent  in  summer.  In  view  of  (c)  and  (d) ,  the  most  useful  data 
was  obtained  in  the  autumn. 

(e)  D  region  absorption  should  be  negligible.  This  con¬ 
dition  was  satisfied  by  confining  attention  to  the  early 
morning  or  late  afternoon. 

(f)  The  F  region  should  approximate  to  a  simple  parabolic 
layer.  This  is  most  nearly  so  in  the  afternoon,  after  the 
FI  layer  has  merged  with  the  F;2  layer,  normally  after  1600 


114 


hours  . 

(g)  The  swept-gain  unit  calibration  should  be  known.  The 
most  reliable  calibration,  described  in  Section  8.4,  was 
performed  in  April,  1963. 

Attention  was  therefore  concentrated  on  afternoon  data 
taken  during  the  period  from  19/3/63  to  16/5/63,  when  back- 
scatter  data  to  the  east  or  west  was  taken  on  38  days.  Of 
these  days,  20  were  capable  of  yielding  results  in  conjunct¬ 
ion  with  vertical  incidence  data  and  calculations. 

8.6  Discussion  of  Calculations 

To  find  the  way  in  which  P„  and  A  might  be  expected 

S 

to  behave  as  f^  decreases  at  sunset,  a  sample  calculation 
was  performed  for  hQ=  200  km,  y^=  50  km  and  f  varying  from 
10  Me  /  s  to  5  i.lc/s;  G  and  y  were  assumed  constant.  The  re¬ 
sults  for  Pr  and  A  were  plotted  against  the  results  for 

O  G 

p'  ,  and  the  f  values  were  indicated  on  the  P„  curve, 
min  o  6 

Fig.  44  shows  these  results.  If  pr  increases  at  a  con- 

min 

stant  rate,  then  the  curve  is  a  measure  of  the  rate  cf 

change  of  P^,  with  time.  The  curve  suggests  that  P^.  varies 
slowly  and  is  fairly  constant  at  long  ranges.  This  is  for¬ 
tunate,  as  it  means  that  provided  A  is  greater  than  zero, 
any  sudden  changes  in  P^  will  be  due  al most  entirely  to 
changes  in  G  and  y 

G  and  y  are  both  functions  of  A  ,  so  it  is  interest¬ 
ing  to  note  the  range  of  values  of  A  that  may  be  expected 


(degrees)  P_  (dB) 
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Pig.  44.  Variation  in  and  A  with  p'  .  as 

c  s  rain 

assuming  G,  -y  »  hQ  and  ym  to  he  constant. 


f  varies, 
o 
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for  typical  data.  When  fQ  =  10  Mc/s ,  Ag  =  33° ,  and  as  f 
decreases,  A  approaches  zero  at  the  rate  of  about  5°/Mc 
(Fig.  45) .  On  many  days  when  records  were  suitable  for 
scaling,  f  was  below  8  Mc/s,  and  therefore  Ag  was  less 
than  about  23°.  At  about  4  Mc/s,  Ag  =  0,  although  even 
at  5  Mc/s  the  computer  program  may  fail,  as  after  Ag  has 
been  computed,  it  may  not  be  possible  to  compute  the  integ¬ 
ral  for  P^  as  A  becomes  smaller.  To  avoid  wastage  of  com¬ 
puter  time,  the  work  was  programmed  to  reject  data  whenever 
it  required  A  to  become  less  than  5°.  The  inability  of  the 
computer  to  cope  with  low  values  of  A  appears  to  manifest 
itself  in  the  Appleton  and  Beynon  equations.  Three  programs 
containing  these  equations  failed  to  cope  with  the  very  low 
angles,  including  programs  written  in  both  GAP  and  WTZ 
source  languages. 

The  program  used  to  find  the  curves  giving  y  from 

m 

h  ,  f  and  p'  also  yielded  the  value  of  A  for  the  p' 

o  o  min  min 

ray.  When  A  s  is  30°,  this  angle  is  1°  or  2°  less  than 

A  .  When  A  1S  small,  this  angle  is  very  slightly  less, 
s  s 

The  values  of  A  for  P'min  have  been  plotted  against  f  for 

h  =  180,  220  and  260  km,  and  various  values  of  y  ,  in  Fig- 

ures  46  to  48.  In  all  these  cases  it  can  be  seen  that  A 

becomes  zero  between  about  4  and  5  Mc/s.  The  rate  at  which 

A  approaches  zero  increases  at  low  angles,  so  that  a  very 

small  decrease  in  f  may  cause  the  p’  ray  to  fall  to  zero 

o  r  mm 

It  is  expected  therefore  that  the  disappearance 


elevat ion  . 


34  56789  10 

f0  (Mc/s) 


Pig.  46.  Variation  in  elevation  angle  of  the  minimum 
path-length  ray  as  f  varies. 
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Pig.  47.  Variation  in  elevation  angle  of  the  minimum 
path-length  ray  as  f  varies. 


3 


4  5  6  7 

f0  (Mc/s) 
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of  the  echo  as  A  approaches  zero  after  sunset  should  occur 
rapidly,  even  if  G  and  V  are  constant.  Any  sudden  decrease 
in  P^  when  A  is  small  is  probably  due  to  this,  as  an  error 
of  only  0.1  Mc/s  could  give  rise  to  an  error  of  3°  in  A  ; 
but  a  sudden  decrease  in  when  A  is,  say,  14°,  can  only 
be  due  to  G  or  y 

Fig.  49  (6/4/63,  80°  mag.  az . )  is  a  graph  of  echo  pow¬ 
er  against  time  of  day.  Each  gain  stage  represents  about  3 
dB.  The  echo  power  is  almost  constant  for  several  hours, 
then  drops  off  rapidly  near  sunset.  Within  half  an  hour, 
the  echo  has  disappeared  from  the  record.  Investigation 
shows  that  in  that  time  the  critical  frequency  dropped  from 
7  Mc/s  to  5.3  Mc/s,  p'  increased  from  1550  km  to  2250  km 

and  A^  dropped  from  17°  to  about  10°.  It  is  that  short 
space  of  time  (not  quite  as  short  for  western  echoes)  which 
can  reveal  most  about  the  variation  of  Y  with  A 

8,7  High  angle  rays 

The  present  calculations,  in  common  witli  SHEARMAN 
(1956b),  have  neglected  the  high  angle  rays,  that  is,  those 

with  angles  of  elevation  greater  than  A  .  Actually  at 

s 

ranges  beyond  the  high  angle  and  low  angle  rays  interfere, 

producing  a  pattern  of  fringes  about  1  km  per  fringe  in 
width  (BREMMER,  1948,  Fig.  72) .  At  ranges  considerably  be¬ 
yond  ,  the  fringes  virtually  disappear,  owing  to  geomet¬ 
rical  considerations  which  cause  the  intensity  of  the  high 
angle  rays  to  fall  away  (BREMMER,  1948,  Fig.  69) .  At  a 
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Pig.  49.  Variation  in  echo  power  with  time  of  day, 

6/4/63,  80°  mag.  az. 


range  100  km  beyond  ,  a  high  angle  ray  is  typically  3  dB 
weaker  than  the  low  angle  ray,  and  at  +  1000  km  it  is  20 

dB  weaker.  In  addition,  the  deviative  absorption  of  the 
high  angle  rays  is  comparatively  very  large  (BIXBY,  1953, 

Fig .  28 ) . 

As  the  present  work  is  concerned  with  echoes  from  near 
the  skip  distance,  it  may  be  assumed  that  high  and  low  angle 
rays  are  equally  strong,  and  that  well-defined  fringes  will 
be  formed,  although  the  width  of  the  transmitted  pulse  pre¬ 
vents  any  resolution  of  the  individual  fringes.  The  maxima 
are  about  twice  as  strong  as  the  low  angle  contribution  al¬ 
one,  that  is,  about  3  dB  stronger  than  the  calculated  inten¬ 
sity.  As  the  results  are  only  relative,  this  3  dB  error  is 
of  no  consequence.  Even  if  the  size  of  this  error  varies 
slightly  as  A  varies,  neglect  of  the  high  angle  rays  does 
not  give  rise  to  significant  errors  in  the  results. 

Therefore  the  high  angle  rays  have  been  neglected  but 
not  overlooked.  Further,  in  the  curves  of  y  against  A  , 
the  high  angle  rays  have  been  taken  into  account  as  follows: 

y  was  plotted  against  Ag,  but  in  any  echo  of  finite 
pulse-width,  the  echo  from  the  skip  distance  is  combined 
with  echoes  from  beyond  the  skip  distance  propagated  by  both 
high  and  low  angle  rays.  The  value  of  Ag  plotted  is 
therefore  not  only  the  angle  of  the  most  powerful  ray,  but 
also  the  approximate  mean  of  all  the  angles  represented. 
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8.8  Magne t o- ionic  splitting 

In  the  calculations,  only  the  ordinary  ray  has  been 
considered.  The  neglect  of  the  extraordinary  ray  is  justif¬ 
ied  in  the  case  of  east-west  propagation.  This  Section 
gives  a  summary  of  the  theory  leading  to  estimates  of  the 
errors  to  be  expected  if  magneto-ionic  splitting  is  neglect- 
e  d  . 

The  formulas  are  expressed  in  the  form  given  by 
RATCLIFFE  (1959) ,  and  the  relevant  symbols  are  as  follows: 

f  =  wave  frequency  =  16  Me / s 

f  =  critical  frequency,  neglecting  the  earth's 
field 

f  =  ordinary  ray  critical  frequency 

fg-  extraordinary  ray  critical  frequency 

f^=  gyrof requency  =  1.4  Mc/s  at  Brisbane  in  the 

F  region 

G  =  angle  between  ray  and  geomagnetic  field 
2  2 

X  =  f  ‘/f 

Y  =  f  /f  =  0.0875  here 

Y^. ,  Yl  =  components  of  Y  transverse  and  longitud¬ 
inal  to  the  ray 

The  treatment  here  is  consistent  with  BIXBY  (1953)  and 
C.R.P.L.  (1948).  Collisions  are  neglected. 

( a )  No  magnetic  field 

2  1  Y 

n  =  1-X 

2  2 

=  1-  fcVf 


(1) 
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The  lower  level  ray  is  the  E  ray,  for  which 


1  -  Y 


(1  +  Y)  f  since  Y  <T<  1. 
c 


(3) 


Also,  we  have 

,  2 


0 


hence 


(  1  +  Y) 


(1+Y)  f 


1.09  f 


(d)  Conditions  intermediate  bet 
between  ( b)  and  (c) 


0  ray:  f  =  f 

o  c 


E  ray:  Most  of  the  E  ray  reflection  takes  place 
at  the  lower  level,  where  X  =  1-Y,  hence  from  (3) , 

fE  -J 


1  +  Y  f 


1.043  f  since  f  =  f  , 
o  o  c 


( e )  Effect  of  neglecting  spiriting 

Suppose  f  -  7  Mc/s,  h  =  200  km,  y  =50  km. 


m 


(i)  Quasi-transverse:  f..=  1.004  f  =  7.028  Mc/s. 

t  o 

Fig . 37  shows  that  for  the  assumed  configuration, 

p'  (0)  =  1500  km,  and  p'  (E)  =  1490  km, 

mm  mm 


(n)  Quasi-longitudinal:  f  = 


f ,=  6.7  Mc/s  (from  (2)) 


°J  1  +  Y  ° 

fE  fc-  7,3  Mc/S  (from  (3)) 


Then  from  Fig. 37  ,  p'  (0)  1610  km  and  p*  (E)  =1430  km, 

'  min  mm 
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In  the  present  work,  the  pulse  width  contributes  to 
the  echo  an  interval  b  p'  of  90  km,  so  for  east-west  pro¬ 
pagation  ( quasi-transve rse )  the  two  parts  of  the  echo  can¬ 
not  be  resolved.  The  only  effect  of  neglecting  the  E  ray  is 

an  error  of  about  10  km  in  p’  ,  which  is  smaller  than  the 

min  ’ 

accuracy  of  reading  p'  from  the  films. 

For  north-south  propagation,  however,  the  rays  app¬ 
roach  the  quasi-longitudinal  condition.  The  E  ray  echo  may 
be  resolved  from  the  0  ray  echo,  especially  if  the  skip  dis¬ 
tance  focusing  is  strong  and  if  the  0  and  E  rays  are  of 
equal  intensity;  in  any  case,  p'  observed  is  that  for  the 

K  ray.  Magneto-ionic  splitting  should  not  be  neglected  in 
north-south  propagation.  Fig.  27  probably  illustrates 
splitting  in  a  IF  echo. 

8.9  Deviative  absorption 

The  calculations  neglect  any  correction  for  deviative 
absorption,  that  is,  the  absorption  that  takes  place  within 
the  F  layer  and  causes  the  F  layer  reflection  coefficient 
to  be  less  than  unity. 

The  amount  of  deviative  absorption  is  believed  to  be 
negligible  in  the  present  work.  An  estimate  of  its  magni¬ 
tude  can  be  gained  from  the  relation 
“  ln  P  =  ( P'-P) , 

where  p  =  reflection  coefficient 


V 


collision  frequency 


in3  -1 
10  sec 


approx,,  for  the 
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\ 


region . 

P' 

=  group 

path 

P 

=  phase 

path. 

This  formula  is  given  by  RATCLIFFE  (1959,  pp  116  and 
132)  and  applies  when  there  is  no  magnetic  field,  and  app¬ 
roximates  to  the  case  for  east-west  (quasi-transverse)  pro¬ 
pagation.  The  loss  of  power  (in  dB)  for  a  ray  from  the 
bottom  of  the  F  layer  to  the  top  of  its  path  is 


L  =  - 


10 


In  0 


2.3 

=  -  4 . 35  In  p 

For  propagation  up  and  down,  the  total  loss  will  be 

L  =  -  8.7  In  p 

8  .7  v  ,  . 

=  {p  “P) 

-  5 

=  1,45  .  10  (p'-p),  if  (p'-p)  is  in  metres. 

For  a  rough  estimate  of  p'-p,  let  us  assume  a  flat 
earth  and  flat  ionosphere.  In  the  usual  notation  (the  sub¬ 
script  1  denoting  the  ionospheric  portion  of  the  path) , 


P'-P  =  Pjl*  -  P  x  - 

p^'  is  simply  sec  A  .  To  estimate  p^,  we  note 

that 

J  DT°  -  <2V2  >  »i  >  Di 

Taking  p^  as  the  r.m.s.  value  of  these  limits, 


The  curves  in  THOMAS  and  MclNNES  (1962,  Fig. 3)  show 


D„  for  various  values  of  A  ,  given  h  =  200  km  and  y  = 

150  km. 

For  Ag  =  25°,  the  penetration  angle  is  about  29°, 
so  the  9  Mc/s  curve  is  appropriate.  For  this,  D^=  650  km. 

Then  p'-p  =  33  km  approximately,  for  which  L=0.5  dB, 

For  A  =  12°,  D^=  1000  km,  p'-p  =  0  approx.,  and 
L  =  0  approx. 

For  the  echo,  the  losses  would  be  txv.ice  as  large. 
Actually,  in  the  present  work  y^  is  nearer  to  50  km  than 
150  km,  therefore  the  absorption  is  less  than  estimated 
here. 

We  may  conclude  that  the  loss  of  echo  power  due  to 
deviative  absorption  is  no  greater  than  1  dB  in  the  present 
wo  rk  . 

8 . 10  Summary 

The  method  of  calculating  the  maximum  power  returned 
in  IF  echoes  takes  into  account  all  known  influences,  in¬ 
cluding  antenna,  ionospheric,  geomagnetic  and  ground  effects. 
Skip  distance  focusing  is  of  great  importance,  while  for  the 
data  chosen  it  was  possible  to  neglect  D  region  absorption, 

E  and  Fl  layer  effects,  deviative  absorption,  high  angle 
rays  and  magneto-ionic  splitting. 

A  graphical  method  of  matching  oblique  backscatter 


data  with  vertical  incidence  data  extrapolated  to  the  con¬ 
trol  point  has  been  developed,  as  well  as  a  satisfactory  way 
to  estimate  the  semi-thickness  of  the  assumed  parabolic 


layer.  The  heoretical  behaviour  of  the  echo  power  at  sun¬ 
rise  or  sunset  has  been  discussed  on  the  assumption  that 
backscattering  is  isotropic. 
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9.  RELATIVE  BACKSCATTER  COEFFICIENT - RESULTS 

Figures  50  to  61  show  P^-P^  (relat  .e  backscatter  co¬ 
efficient  y  )  plotted  against  A  for  most  of  the  days  for 
which  calculations  were  done.  To  save  space,  results  for 
sea  (80°  mag,  az. ,  crosses)  were  sometimes  plotted  with  oth¬ 
er  results  for  land  (260°  mag.  az.,  circles).  Three  notable 
features  emerge. 

9.1  Comparison  of  land  and  sea 

Backscatter  from  the  land  is  much  weaker  than  from  the 
sea.  The  y - A  curves  above  15°  are  approximately  para¬ 

llel,  but  there  is  a  10  dB  difference  between  land  and  sea. 
This  difference  is  much  greater  than  any  likely  errors  in 
the  antenna  patterns  or  in  the  choice  of  the  ionospheric 
parameters,  and  is  therefore  largely  due  to  a  difference  in 
y  .  The  asymmetry  of  east  and  west  propagation  at  sunset 
could  conceivably  have  affected  the  result,  so  a  sample  cal¬ 
culation  was  done  for  both  east  and  west  in  the  morning, 
when  the  asymmetry  was  reversed.  Fig.  59  shows  these  re¬ 
sults;  sea  scatter  is  still  stronger  than  land  by  about  10 
dB.  Morning  calculations  were  not  as  reliable  as  those  of 
the  evening,  owing  to  the  presence  of  the  morning  FI  layer, 
which  rendered  the  Appleton  and  Beynon  formula  not  strictly 
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Figures  50  and  51.  Relative  backscatter  coefficient. 
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Figures  54  and  55 


Relative  backscatter  coefficient 


?0  A  (degrees) 
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o 
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Figures  56  and  57 


Relative  "backscatter  coefficient 
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Figures  58  and  59. 


Relative  "backscotter  coefficient. 


Figures  60  and  01.  Relative  ‘backscatter  coefficient. 
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applicable.  Nevertheless  the  morning  results  were  suffic¬ 
iently  accurate  to  show  that  the  10  dB  difference  between 
sea  and  land  scatter  is  independent  of  ionospheric  asymmetry. 

This  difference  agrees  with  the  observations  of  RANZI 
and  DOMINICI  (1959) ,  although  their  estimate  was  a  very 
rough  one  based  only  on  a  visual  inspection  of  B-scan  re¬ 
cords,  and  was  subject  to  serious  doubts  as  the  behaviour  of 
the  vertical  radiation  pattern  of  their  antenna  varied  mark¬ 
edly  with  azimuth,  and  this  had  not  been  taken  into  account. 

The  present  result  also  agrees  qualitatively  with  1IAGN 
(1962) ,  who  found  land  backscatter  to  be  about  20  dB  weaker 
than  sea,  for  comparatively  smooth  land  surfaces.  In  the 
present  work,  it  is  possible  that  mountains  and  other  large 
features  contributed  to  the  higher  backscatter  coefficient 
than  that  given  by  HAGN . 

The  difference  between  sea  and  land  backscatter  has 
important  implications  which  will  be  dealt  with  in  Chapters 
10  and  11. 

9,2  The  knee  effect 

For  most  of  the  cases  calculated,  as  A  decreases,  y 
declines  suddenly  at  some  angle.  This  effect,  also  describ¬ 
ed  by  HAGN  (1962) ,  is  called  the  knee  effect.  According  to 
HAGN,  the  angle  at  which  the  knee  occurs  is  not  constant, 
even  for  similar  types  of  terrain.  For  example,  in  HAGN, 
the  knee  occurred  at  angles  between  0°  and  25°  for  dry  land 
of  apparently  the  same  type. 
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Here  we  are  concerned  with  only  two  terrains,  land 
and  sea,  and  it  is  probable  that  these  terrains  are  fairly 
constant  in  type,  as  the  location  of  the  sounder  does  not 
change,  and  as  in  the  directions  chosen  the  land  or  the  sea 
extends  to  long  ranges,  without  any  conspicuous  change  in 
type  . 

For  land,  the  angle  at  which  the  knee  occurs  is  usua¬ 
lly  between  14°  and  12°;  in  view  of  the  closeness  of  these 
limits,  the  variation  is  likely  to  be  due  to  experimental 
errors,  and  13°  may  be  taken  as  the  true  and  virtually  in¬ 
variant  knee  angle  for  the  given  terrain. 

The  constancy  of  this  angle  may  be  ascribed  to  the 
fact  that  in  the  present  experiment,  the  terrain  was  the 
same  for  all  measurements  (apart  from  very  small  amounts  of 
scattered  rainfall,  of  the  order  of  10  points  on  some  days) 
whereas  in  HAGH ' s  work,  the  terrain  types  investigated  were 
quite  different  samples  which  were  only  apparently  similar. 

Two  possible  doubts  arise  concerning  the  reality  of 
the  knee  . 

(a)  It  may  be  due  to  failure  of  the  transmission  path;  as 
f  falls,  so  that  A  approaches  0°,  PR  vanishes  as  virt¬ 
ually  all  rays  penetrate  the  ionosphere,  Thi6  might  appear 
to  be  a  knee  if  the  calculations  erroneously  give  values  of 
Ag  that  are  too  large.  But  the  discussion  in  Section  8,6 
shows  that  A  would  not  be  in  error  by  so  much;  a  sudden 

decrease  in  P„  when  A  is  over  10°  can  only  be  due  to  G  ory, 
f<  s 


(b)  The  terrain  may  not  be  uniform  to  long  ranges.  Per¬ 
haps  at  some  range  there  is  an  abrupt  change  in  the  terrain, 
and  y  becomes  very  small  at  long  ranges  (corresponding  to 
low  angles  of  elevation) .  But  reference  to  a  map  of  the 
region  from  Charleville  to  Alice  Springs  and  beyond  (LEWIS 
and  CAMPBELL,  1951)  shows  that  whereas  there  are  large  pat¬ 
ches  of  sandy  or  rocky  desert  interspersed  with  "tall  grass 
savana"  and  "tropical  desert  shrub",  there  is  no  range  at 
which  the  terrain  in  the  illuminated  region  changes  abrupt¬ 
ly;  at  every  range,  the  azimuthal  beamwidth  of  the  antenna 
is  wide  enough  to  include  a  fair  mixture  of  these  types. 

As  far  as  large  features  are  concerned,  the  only  one  of  note 
that  corresponds  in  range  to  the  knee  angle  is  the  Macdonnel 
Ranges  west  of  Alice  Springs,  which  would  be  more  likely  to 
cause  an  increase  in  y  than  a  decrease. 

For  sea,  the  knee  angle  is  less  constant  than  for  land 
and  in  half  the  cases  plotted  the  knee  (if  it  existed)  was 
at  an  angle  below  those  at  which  data  was  available.  When 
the  knee  was  observed,  it  occurred  at  angles  ranging  from 
14°  to  9°.  The  presence  of  u  knee  at  14°  shows  that  in  some 
cases  at  least  the  knee  is  a  real  effect,  although  at  angles 
below  10°  it  could  possibly  be  due  to  errors  as  mentioned  in 
9.2  (a)  above.  But  the  problem  is  not  to  explain  away  a 
knee  observed  at  low  angles,  but  to  explain  its  apparent  ab¬ 
sence  in  some  cases.  Probably  the  knee  is  always  present, 
even  if  it  exists  below  the  angles  investigated.  The  knee 
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angle  can  probably  vary  from  about  13°  down  to  almost  0°. 
This  variation  could  be  due  to  the  state  of  the  sea  surface 
which,  unlike  the  land  to  the  west  of  Brisbane,  is  subject 
to  quite  large  variations.  The  number  of  cases  considered 
is  not  large  enough  to  show  any  correlation  between  knee 
angle  and  sea  roughness.  Dr  D.  F.  Martyn  (private  conver¬ 
sation)  is  convinced  that  the  sea  state  does  have  an  effect 
The  effect  known  to  exist  at  u.h.f.  is  discussed  in  Section 
9.4. 

The  presence  of  the  knee  (especially  when  the  knee 
angle  is  known)  is  an  aid  in  distinguishing  between  differ¬ 
ent  types  of  echoes  such  as  IB  and  IF.  In  the  west,  owing 

6 

to  the  knee,  backscatter  is  negligible  below  10°.  E  at  a 

9 

height  of  110  km  therefore  gives  IE  echoes  with  a  limiting 

s 

range  of  about  1000  km.  In  fact,  western  lE^  echoes  obser¬ 
ved  do  show  a  sharp  cut-off  at  about  this  range.  Incident¬ 
ally,  this  observation  gives  support  to  the  conclusion  that 
the  knee  is  not  due  to  a  change  in  ground  type  at  Alice 

Springs,  since  for  IE  echoes,  the  knee  angle  corresponds 

s 

to  a  quite  different  value  of  D, 

For  echoes  from  the  sea,  the  variation  of  the  knee  an 
gle  means  that  the  maximum  range  of  lE ^  cannot  be  as  easily 
predicted,  but  observations  show  that  1E^  echoes  do  not  us¬ 
ually  extend  beyond  1500  km,  which  corresponds  to  an  elev¬ 
ation  angle  of  about  6°.  In  both  east  and  west,  range 
spreading  may  increase  the  maximum  range  of  the  echoes.  E 

6 
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echoes  are  discussed  further  in  Chapter  10. 

9.3  The  shape  of  the  v - A  curve 

above  the  knee 

The  points  shown  on  Figures  50  to  61  were  collected 
on  Fig.  62,  the  ones  above  the  broken  line  being  for  sea, 
and  those  below  for  land.  Points  which  definitely  fell  be¬ 
low  the  knee  were  omitted.  Points  on  Fig.  59  were  not  used 
(morning  data) ,  nor  on  Figures  60  and  61,  when  ripples  were 
known  to  be  disturbing  the  ionosphere.  Also  omitted  were 
some  points  which  were  due  to  P^  =  40  dB,  the  maximum  value 
capable  of  being  read  from  a  swept-ga.in  frame,  and  which 
could  possibly  have  been  due  to  even  stronger  signals. 

The  curves  superimposed  on  the  points  of  the  figure 

o  o 

represent  the  relation  y  <r  1/cos  ip  (=  l/sin~  A  )  .  For 

both  eastern  and  western  echoes,  this  fitted  better  than 

3 

1/cos  ip  or  1/cos  tp  ,  although  on  individual  days  one  of 
the  latter  sometimes  fitted  better.  In  any  case,  the  re¬ 
sults  are  consistent  with  the  relationship  given  by  KATZ 

and  SPETNER  (1958)  (See  Section  4.3).  The  curves  for  y  •' 

2  3 

1,  1/cos  (p  ,  l/cos^tp  and  1/cos  ip  are  illustrated  in  Fig, 

58  as  curves  (a) ,  (b) ,  (c)  and  (d)  respectively,  and  the 

o 

best  fit  in  this  case  is  for  l/cos^ip 

In  Fig,  €2,  the  scatter  of  the  points  about  their  re¬ 
spective  means  is  not  as  bad  as  may  appear,  as  the  y  scale 
is  merely  relative.  Owing  to  day-to-day  changes  in  trans¬ 
mitter  output,  sea  roughness,  receiver  gain  and  CRO  bright- 
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Fig.  62,  Curves  for  y  «  l/cos^0  superimposed  on  data 
selected  from  Figures  50  to  61.  Points  above  broken  line 
are  for  sea,  and  those  below  are  for  land. 
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ness,  the  points  for  various  days  probably  have  slight  zero 
errors  differing  from  one  day  to  the  next,  and  therefore  can 
legitimately  be  moved  up  or  down  the  y  axis  relative  to  one 
another  for  better  agreement.  This  has  been  done  in  Fig, 63, 
where  the  points  shown  in  Fig.  62  have  been  replotted  to 
give  a  better  fit  to  the  two  curves.  It  should  be  emphas¬ 
ized  that  this  procedure  cannot  be  used  to  find  y  (A)  more 

accurately,  but  is  merely  used  to  illustrate  what  the  scatt- 

2 

er  of  points  would  be  y  ec  1/cos  c f  ,  The  scatter  of 
points  is  3  dB  on  either  side  of  the  curve,  and  this  is 
about  the  accuracy  to  be  expected  if  the  probable  errors  are 
2  dB  in  reading  PR  and  1  dB  in  calculating  . 

9.4  Discussion  on  the  knee  angle 

2 

The  relation  y  1/cos  cp  which  seems  to  fit  the  re¬ 
sults  between  25°  and  10°  is  not  true  at  all  angles.  KATZ 
and  SPETNER  (1958) ,  who  derived  such  a  relation  theoretic¬ 
ally  for  low  angles,  derived  other  relations  for  higher  ang¬ 
les.  Further,  at  A  =  0°,  1/sin2 A  becomes  infinite,  so 
there  must  be  yet  another  relation  for  y  (A)  as  A  approach¬ 
es  zero;  the  knee  effect  is  evidence  of  this.  The  variable 
behaviour  of  the  knee  in  the  present  work  does  not  permit 
any  conclusions  to  be  drawn  concerning  y  (A)  below  the  knee 
but  in  some  cases  y  decreases  at  the  rate  of  at  least  10  dB 
per  degree. 

KERR  (1951,  p.  508)  reported  a  knee  effect  for  sea 
return  at  u.h.f.  He  found  a  relationship  such  as  CS  ^  A 


5 


Pig.  63.  Data  from  Pig.  62  normalized  to  fit  the  curves. 
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2 

below  the  knee,  which  is  the  same  as  yccA  if  A  is  small. 

It  may  be  possible  to  describe  y  (A)  both  above  and 
below  the  knee  by  a  single  function.  For  example,  let 

1/Y  oc  A/A2  +  B  sin2  A  . 

2  2 
This  reduces  to  Y  1/sin  A  when  A  is  large,  and  y  ex  A 

when  A  is  small.  At  the  knee  angle,  neither  term  is  predom¬ 
inant  ,  so  le t 

A/ A2  =  B  sin2  A  for  A  =  A  , 
that  is,  A  =  B  A^^  approximately  (  A^  in  radians)  . 

Then  l/y^A^/A2  +  sin2  A 

This  function  is  plotted  in  Fig,  63a,  for  Ak^=  10  ^radians^, 
that  is,  A^  =  10°12'.  The  slopes  above  and  below  the 
knee  agree  with  observations,  but  the  observed  sharpness  of 
the  knee  is  not  reproduced.  There  appears  to  be  a  discon¬ 
tinuity  associated  with  the  knee.  It  seems  that  no  single 
function  will  suffice  to  describe  y  (A)  at  all  angles. 

HAGN  (1962)  suggested  two  possible  causes  for  the  knee 
effect  —  absorption  of  the  signal  by  the  ground,  or  scatter 
at  an  angle  other  than  back.  He  suggested  that  for  vertical 
polarization,  absorption  could  be  due  to  a  Brewster  angle 
effect  . 

In  the  present  work,  the  echo  is  probably  due  to  the 
vertical  component  of  the  outgoing  wave  at  the  point  of 
scatter  (See  Section  5.4),  so  it  is  relevant  to  consider  the 
Brewster  angle.  On  land,  the  Brewster  angle  occurs  at  about 
A  =  17°  for  dry  land  or  rocky  soil,  and  on  sea  the  pseudo- 
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Brei.ster  angle  is  about  0.5°.  At  the  Brewster  angle,  all 
the  energy  is  absorbed,  and  the  reflection  coefficient 
is  zero  (or  nearly  zero  in  the  case  of  the  pseudo-Brewster 
angle) .  But  this  applies  for  smooth  surfaces,  whereas  it 
is  considered  that  the  backscatter  observed  here  is  due  to 
rough  surfaces  with  specular  type  scattorers,  so  that  echoes 
are  received  even  when  the  value  of  A  is  nominally  equal  to 
the  Brewster  angle.  This  explains  why  the  knee  does  not 
occur  at  the  Brewster  angle,  and  suggests  that  the  knee  is 
not  primarily  a  Brewster  angle  effect,  at  least  in  the 
present  investigation. 

The  second  suggestion  by  HAGN  (1962)  which  is  pre¬ 
ferred  by  him,  is  that  the  energy  at  low  angles  is  practic¬ 
ally  all  reflected  in  the  forward  direction,  and  so  is  not 
observed  in  the  back  direction.  He  suggests  that  there  is 
some  critical  angle  at  which  the  surface  looks  very  smooth. 
This  approach  is  certainly  better  than  the  Brewster  angle 
treatment.  The  behaviour  of  v  must  depend  on  the  surface 
roughness  rather  than  on  the  dielectric  nature  of  the  sur¬ 
face.  In  fact,  at  u.h.f.  there  is  some  evidence  that  the 
knee  angle  depends  on  the  roughness  of  the  sea.  Goldstein 
(KERR,  1951,  pp .  505  ff) ,  working  at  10  cm  wavelength,  found 

that  as  the  sea  becomes  more  rough,  the  knee  angle  decreases 
from  4°  to  1°.  KATZITi  (1957)  suggested  that  if  the  knee 
•angle  A^  =  \/4H,  where  H  is  the  height  of  the  scattering 


irregularities  above  a  horizontal  plane,  then  for  A  less 
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than  there  will  always  be  destructive  interference  be¬ 

tween  the  direct  field  and  that  scattered  by  the  sea  sur¬ 
face.  This  is  true  even  for  vertically  polarized  waves 
(SHEARMAN,  1956b) ,  which  at  low  angles  can  exhibit  phase 
reversal  on  reflection.  A  sudden  reversal  of  phase  could 
explain  the  discontinuity  at  the  knee.  In  the  present  work 
however,  if  wave  heights  are  typically  about  3m,  which  is 
the  rms  wave  height  expected  for  a  wind  of  30  knots  (KATZ 
and  SPETNER,  I960) ,  the  corresponding  value  of  A^  would  be 
about  90°,  which  does  not  agree  with  experiment.  To  make 

less  than  10°,  we  would  need  to  assume  a  wave  height  of 
at  least  30  m,  which  is  very  improbable.  KATZIN's  theory 
therefore  fails  to  predict  the  knee  angle  at  16  Mc/s,  and 
according  to  KATZ  and  SPETNER  (1'HiO)  it  also  fails  at  u.h,f 
For  centimetric  radar,  it  has  been  found  that  A,  is 
proportional  to  \  (BECKMANN  and  SPIZZICHINO,  1963,  pp.  414 
and  417;  KATZIN,  1957) .  This  means  that  whereas  at  10  cm, 
is  between  1°  and  4°,  at  18,5  m  (16  Mc/s) ,  A ^  would 
bo  at  least  190°,  which  does  not  agree  with  experiment. 

The  proportionality  therefore  does  not  hold  over  the  range 
of  wavelengths  considered,  but  it  is  nevertheless  true  that 
A.  is  larger  at  h.f.  than  at  u.h.f.  for  similar  sea  states 

K. 

It  seems  that  no  existing  theory  will  suffice  to  des¬ 
cribe  the  knee,  although  effects  noted  at  u.h.f.  seem  to  be 
qualitatively  true  even  at  16  Mc/s.  All  that  can  be  said 
with  confidence  i9  that  t  lie  knee  angle  does  vary  according 


to  the  roughness  of  the  sea,  and  it  seems  to  take  higher 
values  at  h.f.  than  at  u.h.f. 


9.5  Dependence  of  y  on  sea 
r  oug  hne  s  s 

Although  there  was  insufficient  data  to  establish  any 
correlation  between  Ak  and  the  roughness  of  the  sea,  it 
was  possible  to  consider  the  values  of  y  above  the  knee, 
say  at  20°,  and  compare  these  with  sea  states. 

The  Commonwealth  Bureau  of  Meteorology  at  Brisbane 
was  able  to  supply  data  on  windspeeds  at  sea  level  at  Nor¬ 
folk  Island  and  New  Caledonia.  These  islands  are  at  a  range 
of  about  1500  km  from  Brisbane,  and  in  the  present  work  sea 
backscatter  was  obtained  chiefly  from  the  area  between  these 
places.  The  windspeeds  W  were  averaged  to  give  a  measure  of 

the  average  sea  roughness  at  the  scattering  area.  At  u.h.f. 

2 

it  has  been  found  that  y  is  proportional  to  W  or  W  (BECK¬ 
MANN  and  SPIZZICHINO,  1963,  p.  413) ,  so  it  is  to  be  expected 
that  y  will  be  some  increasing  function  of  W  even  at  16 

•  •ic  s  . 

Before  proceeding  to  a  comparison  of  y  and  W,  we 
should  note  the  difference  between  "sea"  and  "swell".  Sea 
refers  to  waves  which  are  still  in  the  process  of  develop¬ 
ment,  with  sharp  crests.  Swell  consists  of  waves  propagat¬ 
ing  some  distance  from  the  region  where  they  were  generated, 
and  these  have  rounded  crests.  Sea  and  swell  normally  ex¬ 
ist  simultaneously  and  may  move  in  different  directions. 
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Table  5  gives  wave  heights  and  percentage  occurrences 
appropriate  to  an  area  1500  km  east  of  Brisbane,  compiled 
from  U.S.A.F.  (1960,  p.12-44) .  The  corresponding  windspeeds 
have  been  estimated  from  the  percentage  occurrences  and  from 
the  fact  that  a  wind  of  30  knots  generates  waves  about  10 
feet  high  (KATZ  and  SPETNER,  1960) . 


TABLE  5 
SEA  STATES 


State 

height 

(ft) 

occurrence 

(%) 

windspeed 
( knots) 

Low  seas 

1-3 

40 

0-10 

Medium  seas 

3-8 

50 

10-20 

High  seas 

over  8 

10 

over  20 

Low  swell 

1-6 

40 

Medium  swell 

6-12 

40 

High  swell 

over  12 

20 

As  it  is  considered  that  the  backscatter  comes  from 


the  sharp  wave  crests  of  the  sea,  and  not  from  the  swell, 
the  backscatter  coefficient  should  correlate  with  the  sea 
and  hence  with  windspeed.  On  the  other  hand,  the  knee  angle 
probably  depends  on  the  total  wave  height  H  (Section  9.4) , 
which  depends  mainly  on  the  swell;  therefore  the  knee  angle 
would  not  be  expected  to  correlate  with  the  windspeed. 

Table  6  lists  wind  data  at  1800  hours,  Brisbane  time, 
along  with  backscatter  results  from  Figures  50  to  58. 
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TABLE  6 


WINDSPEEDS  AND  BACKSCATTER  RESULTS 


Date 


Wind speed 

_ (knots) _ 

Norfolk  New  Cal-  Average 
Island  edonia 


Backscatter  Results 

y  at  20°  A  k 

(  dB)  (  degrees) 


20/3/63 

10 

E 

10 

E 

10 

(  10) 

171 . 2 

13 

21/3/63 

15 

N 

20 

E 

17 . 5 

(10) 

168.9 

13 

22/3/63 

5 

NE 

5 

E 

5 

(  4) 

170.4 

14 

23/3/63 

15 

NE 

15 

E 

15 

(12.5) 

170.4 

8 

24/3/63 

10 

NE 

15 

E 

12 . 5 

(11) 

168 .7 

<  9 

28/3/63 

20 

E 

20 

E 

20 

(  20) 

17  2 

<10 

5/4/63 

5 

W 

10 

SE 

7.5 

(  5) 

168 .7  ’ 

<  9 

6/4/63 

10 

SW 

10 

SE 

10 

(  5) 

167 

<10 

7/4/63 

10 

SE 

10 

S 

10 

(  3.5) 

167.4 

<  7 

The  average  values  in  the  table  were  calculated  init¬ 
ially  without  regard  to  the  direction  of  the  winds.  The 
bracketed  figures  in  the  table  are  the  "effective"  wind- 
speeds  in  an  east- west  direction.  These  effective  speeds 
were  obtained  by  weighting  the  absolute  windspeeds  accord¬ 
ing  to  direction,  and  then  averaging.  The  weights  were  bas¬ 
ed  on  information  concerning  the  dependence  of  y  on  the 
direction  of  the  wind,  given  by  MACDONALD  (1956,  quoted  in 
BECKMANN  and  SPIZZICHINO,  1963,  p.  413).  If  u,  d  and  c 
denote  upwind,  downwind  and  crosswind  values,  then  for  ver- 


tical  polarization  MACDONALD  found  that  y  y  y  ">  Y  c  > 
and  for  horizontal  polarization  y^  >  y^  y  y^,  the  differ¬ 
ences  involved  in  the  inequalities  being  0  to  5  dB.  It  is 
assumed  that  in  the  present  project  the  polarization  is  ef¬ 
fectively  vertical  and  that  y^  =  y^  and  y^  =  0. 

Then  if  y  is  approximately  proportional  to  W,  W  can  be 
weighted  thus: 

W  =  W  or  +  W  or  O.W  .  For  directions  oblique  to 
u  d  c 

the  east-west  direction,  projected  values  are  used. 

Fig,  64  is  a  plot  of  relative  backscatter  coefficient 
against  windspeed  regardless  of  wind  direction.  Allowing 
for  errors  of  -  3  dB  in  y  (Section  9.3) ,  we  find  that  while 
a  smooth  curve  can  be  drawn  through  the  error  bars,  its 
slope  could  be  positive,  zero  or  negative. 

If,  however,  the  effective  east-west  windspeeds  are 
used,  as  in  Fig,  65,  there  may  be  some  correlation.  Even 
if  the  error  bars  are  reduced  to  -  2  dB,  it  is  still  poss¬ 
ible  to  draw  a  smooth  curve  through  them.  This  curve  has  a 
definite  positive  gradient,  averaging  0.2  dB  per  knot.  The 
better  result  obtained  for  the  effective  windspeeds  justif¬ 
ies  weighting  them  according  to  wind  direction,  and  not  on¬ 
ly  supports  the  hypothesis  that  the  scatter  arises  from  the 
sea  wave  crests,  but  also  suggests  that  scatter  is  most  ef¬ 
fective  on  the  sharper  down-wind  sides  of  the  crests. 

If,  then,  the  backscatter  coefficient  depends  on  the 
sea  rather  than  on  the  swell,  v  should  approach  zero  as  W 
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W  (knots) 


Figures  64  and  65.  Variation  of  backscatter  coefficient 
with  windspeed  (weighted  according  to  direction  in  Fig.  65). 
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approaches  zero.  On  Fig.  65,  y  in  dB  should  approach  -oo  , 
and  the  curve  as  drawn  is  capable  of  such  an  interpretation. 
Probably  for  effective  windspeeds  below  3  knots  the  back- 
scatter  coefficient  changes  rapidly  with  sea  roughness,  but 
above  6  knots  is  not  much  influenced  by  the  t  .ate  of  the 
sea.  A  sea  sufficiently  calm  to  suppress  all  backscatter 
would  be  very  rare. 

Experiments  at  .,h,f.  by  DAVIES  and  MACFARLANE  (1946, 
reported  by  BECKMANN  and  SPIZZICHINO,  1963,  p.  412)  are  in 
general  agreement  with  this  hypothesis.  They  found  that  y 
increased  rapidly  with  wave  height  until  it  reached  a  satur¬ 
ation  value,  when  it  remained  constant.  At  16  Me / s ,  it 
seems  that  saturation  is  virtually  complete  for  W  =  10  knots  , 
that  is,  for  wave  heights  of  about  3  feet. 

9 . 6  Summary 

There  are  three  important  results,  namely,  that  the 
sea  scatter  is  stronger  than  land  scatter  by  about  10  dB; 
that  at  some  angle  (13°  for  land  and  lower  for  sea)  there 
is  a  knee  effect  below  which  backscatter  is  much  reduced; 
and  that  above  the  knee  the  curves  for  both  sea  and  land  are 
consistent  with  a  theory  of  quasi-specular  scatter  from  up¬ 
right  objects  such  as  trees  and  wave  crests. 

The  knee  is  probably  due  to  destructive  interference 
between  the  direct  ray  and  the  ground-re  flee t ea  ray  at  the 
point  of  scatter,  due  to  phase  reversal  of  t  lie  reflected 
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ray  at  low  angles.  The  knee  angle  may  therefore  depend  on 
the  height  of  the  scatterers  above  the  reflecting  plane,  and 
in  the  case  of  sea  scatter  this  height  would  depend  on  the 
swell  rather  than  on  the  sea. 

In  contrast,  the  backscatter  coefficient  at  a  given 
angle  above  the  knee  should  depend  on  the  wave  crests,  that 
is,  on  the  sea  rather  than  the  swell.  The  results  are  not 
inconsistent  with  a  correlation  between  the  backscatter  co¬ 
efficient  and  the  windspeed. 
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10,  OCCURRENCE  OF  SPORADIC-E 

10.1  Comparison  of  occurrences 
over  land  and  sea 

Backscatter  sounding  has  been  used  to  investigate  the 
geographical  distribution  and  movement  of  clouds,  and 
has  yielded  good  results  concerning  the  sizes  and  shapes  of 
E^  clouds  and  their  direction  and  speed  of  travel.  One  re¬ 
sult  of  this  investigation,  however,  is  hard  to  explain: 
the  apparent  predominance  of  E^  over  the  sea  compared  to  its 
occurrence  over  the  land.  The  evidence  for  this  apparent 
behaviour  has  been  published  by  EGAN  and  PETERSON  (1961, 
1962) .  Centroids  of  lE^  backscatter  echo  patches  were  plot¬ 
ted  on  maps,  and  in  a  given  period  the  occurrence  over  the 
sea  was  significantly  greater  than  over  the  land.  This  was 
true  whether  the  data  was  taken  at  Stanford,  U.S.A.,  or  at 
Camden,  N.S.W.;  the  IE  echoes  were  mainly  in  the  west  at 
Stanford  and  in  the  east  at  Camden,  but  in  each  case  over 
the  sea.  The  transmitted  frequencies  were  12,  18  and  30 
Me  /  s  . 

The  results  could  be  due  to  a  difference  in  the  occur¬ 
rence  of  E  or  to  a  difference  in  the  backscatter  coeffic- 
s 

lent.  EGAN  and  PETERSON  rejected  the  second  explanation  by 
referring  to  the  results  of  NIELSON  et  al .  (1960),  which 
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indicated  that  the  coefficients  for  sea  and  land  were  not 
markedly  different,  being  in  fact  equal  at  17°,  a  typical 
elevation  angle  for  H^-propagated  rays. 

The  present  work,  which  shows  that  the  backscatter 
coefficient  for  land  is  10  dB  less  than  for  sea,  does  not 
support  NIELSON  et  al .  It  now  seems  that  the  apparent 
difference  in  E^  occurrence  may  be  ascribed  entirely  to  the 
difference  between  the  backscatter  coefficients  for  sea  and 
land.  lE^  echoes  from  the  land  ene  suppressed  by  two  fac¬ 
tors  -  the  lower  backscatter  coefficient  above  the  knee  in 
the  y  (A)  curve,  and  the  existence  of  the  knee  at  the  com¬ 
paratively  high  elevation  of  13°;  many  E^  clouds  have  a  low 
critical  frequency,  and  will  propagate  16  Me / s  waves  only 
below  13°. 

The  main  features  of  this  conclusion  were  published 
(STEELE,  1964)  but  it  should  be  noted  that  in  the  figure 
the  curves  for  Brisbane  were  obtained  by  a  method  less  acc¬ 
urate  than  that  reported  in  Chapter  8,  and  the  sharpness  of 
the  knee  in  each  curve  has  been  lost,  owing  to  excessive 
smoothing.  It  should  also  be  remarked  that  the  Stanford 
curves  with  which  the  Brisbane  ones  were  compared  are  now 
obsolete,  having  been  superseded  by  those  of  HAGN  (1962) , 
who  in  any  case  prefers  no~  to  compare  sea  and  land  scatter 
below  20° . 

The  conclusion  that  there  is  no  difference  between 
the  occurrences  of  E^  over  land  and  sea  agrees  with  the  ob- 
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sorvations  of  HEISLER  and  WHITE HEAD  (1964)  .  They  used  data 
from  vertical  incidence  sounders  in  many  parts  of  the  world, 
including  both  continental-  and  island-based  stations.  They 
found  no  predominance  of  occurrence  over  the  sea. 

Tho  conclusion  arrived  at  here  disagrees  in  part  with 
the  interpretation  of  DOMINICI  (1962,  1963) ,  who  found  that 
observed  at  Rome  by  ground  backscatter  was  most  prominent 
in  the  north-west  and  least  in  the  north-east.  The  echoes 
north-west  of  Rome  are  scattered  at  the  Atlantic  Ocean,  and 
the  north-east  echoes  on  the  land  in  European  Russia,  as 
illustrated  in  RANZI  and  DOMINICI  (1959) .  DOMINICI  was  of 
the  opinion  that  there  was  a  real  difference  in  E^  occurr¬ 
ence,  although  the  observations  were  influenced  to  a  greater 
or  lessor  extent  by  variations  in  the  backscatter  coeffic¬ 
ient  and  variations  in  the  antenna  vertical  pattern  due  to 
the  topography  at  the  antenna  site.  In  the  present  work, 
the  two  last  mentioned  influences  are  not  m  doubt,  and  the 
backscatter  coefficient  is  entirely  sufficient  to  explain 
the  difference. 

The  distribution  of  1H  echoes  around  Brisbane,  scal¬ 
ed  from  range - azimuth  records,  was  plotted  in  n  similar 

way  to  the  Stanford  and  Camden  results  given  by  EGAN  and 
PETERSON.  The  Brisbane  results  are  shown  in  I  igures  6  6  to 
68.  Range  circles  are  at  1O00  kin  and  2000  Km  ground  range. 
The  points  represent  centroids  of  11  echo  patches,  and 
their  oblique  range  plotted  on  the  figures  differs  from  the 
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Pig.  68.  Data  of  figures  66  and  67  combined 
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ground  range  by  less  than  60  km,  provided  the  oblique  range 
is  greater  than  500  km,  assuming  at  a  height  of  110  km. 
The  transition  from  sea  to  land  conditions  is  quite  clear, 
especially  along  the  Queensland  coast.  In  Fig.  66  the  large 
number  of  centroids  of  IE  patches  near  the  New  South  Wales 
coast  may  be  due  to  a  high  backscatter  coefficient  on  the 
Australian  Alps,  but  more  probably  to  a  bias  towards  this 
direction  in  scaling  data  from  films,  as  the  azimuth  mark 
on  the  films  was  at  185°  mag.  az.,  and  many  centroids  which 
were  at  185°  -  10°  were  ascribed  to  185°.  Also,  there  was 
probably  a  genuine  predominance  of  in  the  south,  which 
often  appeared  additionally  on  the  records  as  direct  back¬ 
scatter  from  field-aligned  E^  at  half  the  range  of  the  E 
echo.  For  the  data  of  Fig.  67,  the  azimuth  mark  was  at 
117°.  Fig.  68  combines  the  results  of  Figures  66  and  67. 

The  main  points  of  difference  between  the  results  for 
sea  and  land  are  that  the  land  echoes  are  less  frequent  and 
extend  to  shorter  maximum  ranges  than  the  sea  echoes.  Both 
these  effects  are  consistent  with  what  is  now  known  about 
the  backscatter  coefficient, 

10,2  Comparison  of  IF  echo  strengths 
on  range - azimuth  records 

The  difference  between  sea  and  land  scatter  which  is 
apparent  in  IE  records  is  just  as  pronounced  for  IF  echoes. 

Fig.  69  is  a  typical  range- - azimuth  record  with  range  mark 

every  250  km.  Superimposed  on  this  record  is  a  rang*' - 
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azimuth  map  of  Australia  drawn  to  the  same  scale.  This  map 
is  shown  in  greater  detail  in  Fig.  70,  where  the  azimuth 
scale  is  relative  to  geographic  north  at  Brisbane.  A  short 
line  drawn  on  the  bottom  of  the  record  indicates  the  posit¬ 
ion  of  the  azimuth  mark  (185°  mag.  az.).  The  record  shows 
the  oblique  range,  but  the  map  is  drawn  in  terms  of  ground 
range;  this  leads  to  a  discrepancy  of  about  200  km,  depend¬ 
ing  on  the  range.  A  rough  correction  has  been  made  by  mov¬ 
ing  the  map  200  km  up  the  page. 

It  may  be  seen  that  there  is  a  thinning  out  of  F  reg¬ 
ion-propagated  ground  scatter  coinciding  with  scatter  from 
the  Australian  continent.  The  regular  fading  characteristic 
of  sea  echoes  and  the  random  fading  of  land  echoes  are  also 
evident.  Even  though  the  antenna  pattern  varies  slightly 
from  east  to  west,  there  is  such  a  strong  correlation  of 
backscatter  strength  with  the  sea  that  antenna  differences 
are  negligible. 

This  relative  weakness  of  IF  echoes  from  the  land  was 

observed  on  range - azimuth  records  at  Rome  by  RANZI  and 

D0MINICI  (1959)  ,  though  with  less  azimuthal  resolution,  and 
with  some  misgivings  about  the  vertical  pattern  of  the  an¬ 
tenna  used  (DCNilNICI,  1962,  1963). 

These  results  strengthen  the  case  against  any  predom¬ 
inance  of  over  the  sea,  as  the  same  reasoning  which  led 

to  that  theory  would  lead  to  a  predominance  of  F  region  ion¬ 
ization  ovet  the  sea,  which  is  known  to  be  untrue. 


180  270  360  90 


AZIMUTH  ANGLE  (degrees, geographic) 


Fig.  70.  Range — azimuth  map  used  in  Fig.  69 
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10.3  Occurrence  of  E  over 
_ s _ 

mountain  s 

EGAN  and  PETERSON  (1961,  1962)  also  drew  attention  to 

the  apparent  predominance  of  E  over  mountainous  regions. 

G 

They  admitted  that  this  could  well  be  due  to  a  higher  back- 
scatter  coefficient.  Mountain  ranges  do  not  in  general  give 
rise  to  a  stronger  echo,  unless  they  lie  broadside  on  to  the 
direction  of  propagation  (See  Section  2.2).  In  the  present 
investigation,  no  such  echoes  are  possible,  as  practically 
all  mountain  ranges  are  aligned  radially  with  respect  to 
Brisbane.  The  effect  of  the  Great  Dividing  Range  could  not 
be  investigated,  as  it  is  close  to  the  coastlines  and  it  is 
not  possible  to  eliminate  a  certain  amount  of  sea  scatter 
from  the  results;  by  its  strength,  sea  scatter  would  tend 
to  obscure  any  effect  due  to  the  mountains.  This  was  also 
true  of  New  Zealand,  where  the  high  ranges  of  South  Island 
might  constitute  a  coherent  echoing  surface;  beam  and  pulse 
width  were  not  small  enough  to  resolve  the  island. 

In  general,  mountains  may  be  expected  to  enhance  the 
backscatter  coefficient  in  two  ways:  by  an  average  tilt 

which  might  effectively  decrease  the  knee  angle,  and  by  an 
increase  in  the  number  of  scatterers  per  unit  area.  The 
first  way  seems  likely  only  at  the  near  side  of  a  range  of 
mountains,  preferably  broadside-on.  The  second  seems  un¬ 
likely,  as  the  main  scatterers  are  believed  to  be  objects 
such  as  trees,  which  are  not  necessarily  more  numerous  on 
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mountains . 

Insofar  as  it  is  possible  to  come  to  any  conclusion 

about  this,  it  may  be  noted  that  in  Figures  66  to  68  the 

relative  occurrence  of  It  echoes  seems  to  be  associated 

s 

with  the  coastlines  and  not  with  the  mountains.  If  the 
mountains  do  increase  the  backscatter  coefficient,  the  eff¬ 
ect  is  small . 

10,4  The  "south-east  patch" 

In  early  work  on  field-aligned  irregularities  ob¬ 
served  by  backscatter,  attention  was  drawn  to  an  apparent 
predominance  of  ionization  south-east  from  Brisbane,  manif¬ 
esting  itself  as  a  patch  echo  on  range - azimuth  records. 

Figures  66  to  68  were  prepared  partly  to  investigate  the 
reality  of  this  patch,  and  the  result  was  negative.  All 
patches  which  were  not  obviously  IF  echoes  or  direct  back¬ 
scatter  front  field-aligned  H  were  plotted,  and  no  predom¬ 
inance  of  patches  appeared  in  the  south-east.  There  was, 
however,  a  larger  number  in  the  south-east  than  in  the 
south-west,  and  as  has  been  shown,  this  can  be  ascribed  to 
the  relative  backscatter  coefficients  of  che  ground.  It  is 
now  considered  that  the  "south-east  patch"  is  really  the 
normal  occurrence,  and  that  it  only  achieves  prominence 

when  attention  is  concentrated  on  southerly  directions,  as 
is  done  when  scaling  echoes  from  field-aligned  irregularit¬ 
ies.  When  t he  original  records  which  suggested  the  patch 
were  examined,  there  was  nothing  to  suggest  that  this  was 
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anything  other  than  a  normal  11  echo. 

10 . 5  Summary 

Although  it  has  been  thought  that  occurs  predomin¬ 
antly  over  the  sea,  the  present  work  shows  that  this  belief 
was  largely  unfounded.  Backseatter  records  suggest  an  app¬ 
arent  predominance  which  may  now  be  accounted  for  as  due  to 
the  difference  in  the  backseatter  coefficient.  A  similar 
apparent  predominance  is  observed  in  the  case  of  IF  echoes. 

There  did  not  appear  to  be  any  predominance  of  E  ov¬ 
er  mountainous  regions,  but  resolution  was  insufficient  for 
any  certainty  about  this.  From  what  is  now  believed  about 
the  behaviour  of  the  backseatter  coefficient,  it  may  be 
predicted  that  the  knee  in  the  backseatter  curve  would  be 
lower,  arid  therefore  echoes  might  be  possible  from  longer 
ranges. 
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11.  TRANS EQUATORIAL  PROPAGATION 

11.1  Introduction 

Long-range  one-hop  propagation  to  the  north  of  Bris¬ 
bane  has  been  described  by  THOMAS  (1961,  1962a,  1962b,  1964) 
and  by  THOMAS  and  McINNES  (1962,  1964) .  Backscatter  paths 
between  the  directions  290°  Mag.  and  70°  Mag.  are  consider¬ 
ed.  The  typical  transequatorial  echo  appears  in  the  late 
afternoon  towards  the  north-ea6t,  first  as  a  weak  echo  at 
an  azimuth  of  60°  to  70°,  and  at  a  range  of  9000  to  10,000 
km.  It  gradually  extends  towards  the  north  at  closer  ran¬ 
ges  ,  appearing  as  the  right  arm  of  a  L’  on  the  range - azi¬ 

muth  plot.  It  may  remain  in  position  for  3  or  4  hours,  ly¬ 
ing  roughly  along  the  line  of  magnetic  latitude  20°  north, 
as  shown  in  Fig.  71, 

The  predominance  of  the  transequatorial  echo  in  the 
north-east  is  stressed.  Peak  occurrence  is  at  10°  to  30°, 
and  very  little  of  the  echo  is  detected  further  west  than 
about  320°  (TUC..AS,  1962a). 

11.2  Interpretation 

The  reason  for  this  east-west  asymmetry  is  partly  giv¬ 
en  by  THOMAS  (1962b) .  Analysis  of  ionospheric  height  con¬ 
tours  shows  that  an  upward  bulge  favourable  to  long  range 
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270  300  330  360  30  60  90 

AZIMUTH  ANGLE  (degrees,  magnetic) 

Fig.  71.  Range — azimuth  map  showing  the  relation  of 
transequatorial  echo  occurrence  at  20°  north  magnetic 
latitude  to  the  continent  of'  Asia. 
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propagation  occurs  progressively  between  70°  and  310°,  but 
not  at  290°,  According  to  THOMAS  et  al .  (1962) ,  the  bulge 

finally  disappeared  at  long  ranges  at  about  300°.  But  this 
does  not  entirely  explain  the  decrease  in  occurrence  of 
transequatorial  echoes  at  320°. 

It  might  be  suggested  that  the  layer  tilts  which  per¬ 
mit  transequatorial  propagation  have  an  east-west  asymmetry 
that  militates  against  propagation  to  the  north-west  at  ab¬ 
out  310°,  even  though  the  bulge  is  present.  This  suggestion 
seems  unlikely,  as  in  the  late  afternoon  and  evening  very 
strong  interference  appears  on  the  records  at  310°,  This 
has  sometimes  been  identified  as  the  Overseas  Telecommunic¬ 
ations  Comrission  service,  and  the  direction  of  arrival  in¬ 
dicates  that  the  signals  come  probably  from  London,  or  at 
least  from  Hongkong  or  Manila.  This  shows  that  long  dist¬ 
ance  propagation  from  the  north-west  is  particularly  good 
by  at  least  one  mode,  which  must  cross  the  equator  in  the 
tilted  region  of  the  ionosphere  to  the  north-west  of  Bris¬ 
bane.  A  tilted  region  which  allows  transequatorial  propag¬ 
ation  to  the  north-east  and  north,  and  later  allows  good 
propagation  of  some  kind  to  the  north-west,  would  surely  be 
expected  to  allow  t r an soq u a t o r ia ]  propagation  to  the  north- 
we  s  t  . 

It  might  also  be  suggested  that  I)  region  absorption 
is  high  near  the  scattering  region  in  the  north-west,  3ut 
geographical  symmetry  suggests  that  this  absorption  would 


be  just  as  high  in  the  vicinity  of  the  sounder  at  the  time 
of  onset  of  transequatorial  propagation  to  the  north-east, 
so  that  the  arms  of  the  U  should  nevertheless  be  equal. 

It  is  now  suggested  that  the  rarity  of  transequator¬ 
ial  echoes  at  about  310°  is  duo,  wholly  or  in  part,  to  the 
low  backscatter  coefficient  of  the  land  to  the  north-west. 
Examination  of  Fig.  71  shows  that  at  the  range  at  which 
transequatorial  echoes  terminate  (7000  km,  320°  Mag.)  lies 
the  coast  of  China,  and  further  west  the  line  of  20°  north 
magnetic  latitude  remains  on  the  land,  from  Shanghai,  ac¬ 
ross  China  and  Tibet  to  Delhi  at  a  range  of  10,000  km  from 
Brisbane.  Even  if  the  sea  covered  this  region,  we  would 
expect  transequatorial  echoes  to  be  weak,  by  comparison 
with  equidistant  echoes  to  the  north-east.  The  10  dB  diff¬ 
erence  between  land  and  sea  backscatter  is  therefore  suff¬ 
icient  to  suppress  practically  all  the  transequatorial  echo 
in  the  north-west.  Even  more  significant,  however,  is  the 
fact  that  the  more  important  of  the  transequatorial  paths 
occur  at  low  elevation  angles  (THOMAS  and  McINNES,  1962, 
1964)  ,  especially  below  8°.  This  means  that  the  most  pow¬ 
erful  transequatorial  signals  will  arrive  at  the  ground  at 
an  angle  below  the  knee  angle  for  land,  and  will  not  be  re¬ 
turned  to  the  sounder. 

In  further  support  of  this  explanation,  it  is  observ¬ 
ed  that  in  the  north-west  not  only  are  transequatorial  ech¬ 
oes  scarce,  but  also  IF,  2l;  and  3F  (when  the}  occur  to  the 


17  4 


north)  (THOMAS,  1962) .  If  the  lack  of  echoes  in  the  north¬ 
west  were  due  to  the  shape  of  the  equatorial  bulge  alone, 
this  would  not  explain  the  shortage  of  IF  echoes,  which  occ¬ 
ur  at  ranges  short  of  the  bulge.  The  IF  echoes  are  markedly 
scarce  west  of  about  320°  Mag.,  which  coincides  with  the 
Queensland  coastline.  It  seems  that  the  geographical  asym¬ 
metry  in  the  occurrence  of  transequatorial  echoes  not  only 
accompanies  an  ?. symmetry  in  IF  occurrence,  but  has  the  same 
explanat ion . 

11.3  Summary 

The  comparative  lack  of  one-hop  transequatorial  echoes 
to  the  north-west  of  Brisbane  may  be  ascribed,  wholly  or  in 
part,  to  the  fact  that  the  north-west  echoes  would  be  scatt¬ 
ered  on  the  land  between  Shanghai  and  Mew  Delhi,  and  there¬ 
fore  would  be  much  weaker  than  those  returned  from  the  sea 
to  the  north  and  north-east. 


12.  THE  PROSPECTS  OF  OBTAINING  DISTANT 


IONOSPHERIC  DATA  FROM 
BACKSCATTER  SOUNDINGS 

12,1  Historical  background 

The  usefulness  of  backscatter  sounding  to  find  the 
maximum  usable  frequency  for  a  given  communication  path  (or 
conversely,  the  skip  distance  obtaining  on  a  given  frequen¬ 
cy)  was  foreseen  even  before  the  backscatter ing  range  was 
definitely  identified  with  the  ground  (SMITH,  1945;  BENNER, 
1949).  Later,  the  backscatter  P.P.I,  display  gave  promise 
of  a  simple  "communication  zone  indicator"  ( de  BETTENCOURT, 
1952)  to  give  instantaneous  information  on  possible  commun¬ 
ication  zones  and  modes  of  transmission. 

PETERSON  (1951)  and  WILKINS  and  SHEARMAN  (1957)  poin¬ 
ted  out  the  need  for  great  care  in  deducing  D  from  p' 

s  inj  n 

In  particular,  the  accuracy  of  the  determination  depends  on 

the  correct  identification  of  modes  (IF  or  IE  )  and  the  eff- 

s 

ect  of  off-beam  echoes  on  the  observed  p'  .  With  the  an- 

r  min 

tennas  and  techniques  of  the  present  project  these  problems 
are  minimized,  but  they  would  be  rather  serious  if  inexpen¬ 
sive  equipment  were  used.  In  SHEARMAN  (1961,  discussion) 
it  was  mentioned  by  P.A.C.  Morris  that  backscatter  sounding 
was  successfully’  used  to  find  the  fade-out  time  of  fixed- 
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frequency  transmissions  between  Singapore  and  Hongkong,  by 
pulsing  the  spare  sideband  of  a  normal  multiplex  telegraph¬ 
ic  circuit.  However,  backscatter  sounding  has  not  been 
widely  adopted  as  an  operational  tool  for  commercial  serv¬ 
ices. 

Even  more  ambitious  is  the  proposal  to  obtain  iono¬ 
spheric  par  a me  ters  such  as  f  ,  h  and  y  from  backscatter 

o  o  m 

sounding.  This  would  be  very  desirable,  as  routine  vertic¬ 
al  incidence  sounding  is  possible  over  less  than  a  quarter 
of  the  earth's  surface.  Oblique  sounding  could  possibly  be 
used  to  extend  our  measurements  of  the  ionosphere  into  isol¬ 
ated  regions, 

SHEARMAN  (1956b)  pointed  out  that  measurements  of  p'min 
do  not  give  unique  values  of  layer  height  and  critical 
frequency.  To  supplement  the  p'  observations,  therefore, 

EGAN  (1960)  and  SHEARMAN  (1961)  recorded  echo  pulse  shapes 
(A-scan)  and  compared  them  with  calculated  pulse  shapes. 
SHEARMAN  calculated  the  echo  pulse  shapes  using  values  of 
hQ ,  P  and  y ^  given  in  ionospheric  prediction  maps,  and 
correcting  these  values  so  that  they  agreed  with  actual  ver¬ 
tical  incidence  measurements.  Although  he  neglected  ionos¬ 
pheric  focusing,  D  layer  absorption,  magneto-ionic  splitting 
and  variations  in  the  backscatter  coof f icient ,  he  found  good 
agreement  with  observations. 

EGAN  (1960)  neglected  splitting  and  variations  in  the 

backscatter  coefficient,  and  assumed  that  h  and  y  were  the 

’  o  Jm 


same  at  the  control  point  as  measured  above  the  backscatter 
sounder.  He  then  calculated  the  echo  shapes  and  ranges  for 
various  values  of  f  ,  and  selected  the  one  which  agreed  best 
in  range  and  shape  with  that  observed  by  backscatter.  The 
value  of  f  obtained  in  this  way  agreed  well  with  that  de¬ 
duced  by  interpolating  vertical  incidence  data. 

12.2  Comments  on  previous  work 

(a)  For  reasons  given  in  earlier  chapters  of  the  present 
work,  the  calculation  of  the  complete  echo  shape  can  only  be 
done  accurately  by  taking  into  account  all  known  factors, 
including  antenna,  ionospheric,  geomagnetic  and  ground  eff¬ 
ects.  While  this  would  be  theoretically  possible,  the  am¬ 
ount  of  data  required  would  be  prohibitive  except  in  specia¬ 
lly  selected  cases.  The  measure  of  success  obtained  by  EGAN 
(1960)  and  SnEAR.'.iAN  (1961)  is  thought  to  be  largely  fortuit¬ 
ous  , 

( b)  Even  if  the  calculation  could  be  done  accurately,  it 

would  still  not  be  possible  to  deduce  unique  values  of  fQ, 

h  and  y  .  An  accurate  knowledge  of  p'  could  allow  the 

o  in  min 

deduction  of  one  of  these  parameters  provided  the  other  two 
were  known.  Figures  35  to  40  could  be  used  for  this.  The 
shape  of  the  echo  pulse  might  allow  the  deduction  of  two 
parameters  if  the  third  were  known,  but  the  shape  is  not 
very  sensitive  to  changes  in  these  parameters,  so  the  resul¬ 


ts  would  be  rather  inaccurate. 


12.3  Future  possibilities 

While  the  calculation  of  echo  pulse  shapes  holds  litt¬ 
le  promise  for  the  deduction  of  ionospheric  parameters,  the 
measurement  of  angles  of  arrival  of  backscatter  echoes  at 
the  receiver  should  make  possible  the  deduction  of  two  par¬ 
ameters  if  the  third  is  known.  Calculated  values  of  A  , 

s 

or  of  A  corresponding  to  the  p’^  ray,  could  be  marked  al¬ 
ong  the  curves  of  Figures  35  to  40.  Then,  for  example,  giv¬ 
en  h  ,  we  could  use  the  observed  values  of  p '  .  and  A  to 

o  min  s 

read  off  unique  values  of  f  (or  f„)  and  y  . 

o  E  m 

A  could  be  measured  by  strobing  the  echo  pulse,  and 
measuring  the  angle  of  arrival  of  the  skip  ray;  the  compon¬ 
ent  of  the  echo  pulse  due  to  the  skip  ray  is  concentrated 
around  the  highly-focused  region  just  beyond  P'm^n»  but  oth¬ 
er  angles  will  also  contribute  to  this  section  of  the  echo. 
The  angle  of  arrival  of  the  P*m^n  ray  might  be  more  useful; 
it  could  be  found  by  selecting  the  leading  edge  of  the  echo 
pulse  and  measuring  its  angle  of  arrival.  Such  a  sample 
would  be  easier  to  isolate  by  strobing  and  would  contain  a 
narrower  spread  of  angles.  Unfortunately  the  power  contain¬ 
ed  in  this  part  of  the  echo  is  relatively  weak. 


The  accuracy  of  this  method  may  be  estimated  by  con¬ 
sidering  typical  data.  Suppose  p'  =  1500  -  20  km  and  A 

mm 

for  the  p’  ray  =  1S°  -  1°.  Then  if  we  are  given  h  = 

min  6  o 

200  km,  the  values  of  f  and  v  may  be  read  from  Fig,  72. 

’  o  -  m 

The  shaded  area  on  the  figure  shows  the  limits  of  probable 
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Pig.  72.  Part  of  Pig.  37  with  additional  curved 
parametric  in  A  appropriate  to  the  minimum  path-length 
ray,  illustrating  the  deduction  of  f  and  y  . 
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error  in  the  data,  and  hence  the  limits  of  error  of  the  re¬ 
sult,  We  have,  f  =  6.84  -  0.4  Mc/s  and  y  =45-15  km. 

o  m 

If  were  known,  a  similar  method  would  permit  the 

deduction  of  f  and  h  .  The  error  to  be  expected  in  h 
o  o  c 


would  be 


u  ^ 

6  h  -  — - 7——  *  for  constant  f  and  p'  .  . 

o  dy  / d  h  o  min 

Jm'  o 

From  Section  8,3  ( b)  and  Fig.  41  it  is  seen  that  dy^/ 
dh^  typically  equals  -0,75,  so  6  hQ  =  +  20  km. 

More  generally,  if  b  A  is  the  error  in  A  ,  in  degrees, 
6fQ  =  io.4  6  A  Mc/s, 

b  ym  =  -  15  6A  km  (given  h^), 

S h  =  +  20  Q A  km  (given  y  ) . 
o  m 

Much  depends  on  the  accuracy  with  which  A  can  be 

measured . 

For  sporadic-E,  which  may  be  assumed  to  exist  in  a 
thin  layer,  we  can  take  y^  -  0.  If  the  rays  penetrate  the 
layer  at  high  angles,  P'm^n  and  Ag  could  give  unique  val¬ 
ues  of  hQ  and  fQ,  at  the  point  of  penetration.  If  the  rays 
do  not  penetrate  the  cloud,  as  when  the  cloud  is  at  a  long 
range  from  the  sounder  and  the  rays  meet  its  nearest  edge 

at  a  small  grating  angle,  p'  will  then  be  a  measure  of 

min 

the  range  of  the  cloud,  and  not  of  its  critical  frequency. 

Fig.  73  is  a  set  of  curves  parametric  in  h  ,  with  y 

o’  J  m 

=  0,  for  use  with  Sporadic-E  data.  For  example,  suppose 


300  400  500  600  700  800  900  1000  1100  1200 

Pig.  73.  Curves  to  find  f  K  and  h  E  from  1G  Me/s  data. 
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p'  =  750  -  20  km  and  A  =  15°  1  1°.  Then  f  =  5.06  - 
min  s  o 

0.2  f'ic/s  and  hQ  =  107  -  8  1cm. 

More  generally ,  the  errors  are 
0 f  =  -0.2  6  A  Me  /  s  , 

Oh  =  -  8  £»A  km . 
o 

If  f  is  not  16  Me / s  ,  the  curves  of  Fig.  73  may  still 
be  used  if  the  values  of  P  obtained  are  multiplied  by  f/16. 

Swept- frequency  sounding  is  necessary  if  the  ionos¬ 
pheric  parameters  at  particular  ranges  arc  to  be  measured 
continuously,  since  with  fixed  frequency  sounding  the  con¬ 
trol  point  range  varies  considerably  as  the  ionospheric  con¬ 
figuration  varies.  Angle  of  arrival  measurements  would  be 
very  difficult  with  swept-f requency  backscatter ,  owing  to 
the  problem  of  keeping  the  angle  of  arrival  equipment  in 
tune  with  the  sounder.  A  crude  estimate  of  angle  of  arrival 
could,  however,  be  obtained  by  noting  that  for  16  Me / s  back- 
scatter  from  the  land  west  of  Brisbane  the  knee  angle  is  re¬ 
markably  constant  at  about  13°.  The  maximum  value  of  p' 

mm 

for  the  western  echo  (at  its  appearance  or  disappearance) 

would  therefore  correspond  to  an  angle  A  approximately 

8 

equal  to  10° . 

12.4  Summar y 

The  early  promise  that  backscatter  sounding  would  be¬ 
come  a  widely  used  tool  in  obtaining  maximum  usable  frequen¬ 
cies  for  commercial  circuits  has  not  been  realized,  owing  to 
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difficulties  in  interpreting  the  data.  Similarly,  the  de¬ 
duction  of  the  ionospheric  layer  parameters  f ^ ,  h^  and  y^ 
by  comparing  calculated  and  observed  shapes  of  echo  pulses 
is  highly  complex  and  can  achieve  little.  By  measuring  the 
angle  of  arrival  of  selected  parts  of  the  echo  pulse,  it 
should  be  possible  to  deduce  two  parameters  if  the  third 
were  known;  this  could  be  most  useful  in  observations  of 
sporadic-E . 
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APPENDIX  I 

LIST  OF  PRINCIPAL  SYMBOLS 

a  2  2. 

A  =  x  y  n  /  r 

J  n  '  o 

A  =  O 

o  o 

C  =  cos  i 

o 

c  =  free-space  velocity  of  electromagnetic  waves. 

D  =  ground  range  to  scatter  source,  =  .  * 

D  =  skip  distance, 

s 

IE  =  echo  returned  by  one-hop  E  mode, 
s  s 

F  =  focusing  factor. 

IF  =  echo  returned  by  one-hop  F  mode, 

f  =  transmitted  frequency  (always  16  Mc/s  here) . 

f  =  "no  field"  critical  frequency. 

critical  frequency  of  the  E  ray. 

f^  -  gyrof requency . 

f  =  critical  frequency  of  the  ordinary  ray. 

G  --  power  gain  of  aerials  for  transmission  or  reception, 

e  =  n  /  2  -  it) 

H  -  height  of  ocean  waves. 

h  true  height  of  base  of  parabolic  layer, 

h  '  =  virtual  height  of  base  of  parabolic  layer. 


i 

o 


L 


n 


P 

P' 


P' 

P* 

r* 

‘  O 
b  t 
W 
X 


man 


m 


A 

A, 


angle  of  incidence  of  rays  on  base  of  layer. 

loss  of  power  due  to  deviative  absorption, 
sin  i 

o 

calculated  relative  peak  echo  power  if y  is  constant, 
relative  power  received, 
phase  path-length. 

equivalent  free-space  path-length  to  scatter 
source,  sometimes  referred  to  as  oblique  or 
group  range. 

minimum  value  of  p*  . 

radius  of  the  earth  (about  6370  km) . 

R  +  h 

o 

duration  of  transmitted  pulse, 
windspeed . 

2  o 

f  Vf~ 
f  /  f 

o 

f„/f 

semi-thickness  of  a  parabolic  layer. 

backseat  ter  coefficient,  -  n  /  cos  O 

o 

angle  of  elevation  of  any  ray  at  the  ground , 
knee  angle  in  y  (A)  curve. 

angle  of  elevation  of  the  skip  ray. 


0 


angle  between  ray  and  geomagnetic  field. 
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\  =  wavelength  (At  16  Me / s  ,  this  is  18,75  m)  . 

V  =  collision  frequency. 

p  =  y/  2  ,  also  reflection  coefficient  of  F  region. 

0  ^  =  radar  cross-section  por  unit  area  of  the  ground. 

(f  =  Tt  /  2  -  A 


*  Subscripts  1  and  2  denote  ionospheric  and  sub- 

ionospheric  contributions  respectively. 
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APPENDIX  II 


FORMULAS  USED  TO  CALCULATE  P0  AND  A 

C  s 


( a )  Appleton  and  Be^non  Equations 

As  given  in  SHEARMAN  (1956b) ,  these  are: 


x  cos  1 


2x  y  tanh  ; 
m  ( 


2 

x  y 


m  sin  i 


R  +  h 


n - r~  sin  i  p, 

R  +  h  o  1 


sin  (  H/2  -  A  -  i  ) 


2  R  (  n/2  -  A  -  iQ  ) 


Using  other  symbols  given  in  Appendix  I,  we  can  simp¬ 
lify  the  form  of  these  equations.  Also,  the  further  formula 

„  R  -  h 

R  =  _ o _ 

sini  sin  (  71  /  2  +  A  ) 

o 

enables  i^  to  be  expressed  in  terms  of  A  . 


Thus ,  sin  l 


cos  A  =  n 


Also,  since  g  -  n/2  -  i  , 

’  "  o 


tan  g  =  cot  iQ  -  -  , 

-1  ,C 

g  tan  ( - ) 
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Hence,  the  equations  may  be  written 


pl' 

=  2 

X 

ym 

m 

tanh 

(  x 
(1  ■ 

H 

Q 

R 

n 

r 

o 

pi‘ 

1 

=  2 

p 

sin 

(  R  - 

A  ) 

P2 

n 

Cl 

Q 

=  2 

R 

(  g 

-  A  ) 

( b)  To  find  the  penetration  angle 

By  the  secant  law  of  ionospheric  reflection,  at  the 
penetration  angle,  x  =  sec  i  . 


T hen  cos  i 

o 


so  t  ha  t  n 


-  1  /  x  =  J  1  -  n2 

“J  1  -  1/x2  . 


> 


r  n 

But.  cos  A  =  — ^ —  ,  hence  A  for  penetration  can  be 
found;  this  is  the  angle  of  elevation  at  the  ground  for 
which  a  ray  penetrates  the  ionosphere. 


( c )  Focusing  F actor 

(i)  First  order  theory:  Re-writing  SHEARMAN  (1986b), 
and  simplifying, 


dD 
d  A 


cot  A 


D  dD 
R  d  A 


where 


R  sin  A 

!•  C 

o 


r  (1 
A 

o 

o 

x"y 

J  m 

A  )  > 

+  1  f 

2  R 

r 

o 

- 

V 

( 1 

-  A) 

o  o 

-  x4' 

2  R 

n 
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(ii)  Second  order  theory,  used  at  skip  distance: 

.-12 


10.59  .  10 


cos  A 


\ 


1/3 


sin  A 


sin<R>  <  II 


>■  A 


D 


D 


( d )  Echo  power 

SHEARMAN  (1956b)  has  the  equivalent  of 
D 

PR  cc  f  G2  l-2  Aq  R  sin  .|  dD 

D-  6D 

For  a  single  hop,  R  sin  is  approximately  equal  to  D,  to 
within  2%  and  this  leads  to  an  error  of  only  0.1  dB  in  . 
Making  this  approximation,  and  replacing  A^  by  y  sin  A  ,  we 
have 


D 

PRocf  G  “  F 2  y  sin  A  D  dD 

D-  6D 


in  the  calculations,  y  was  assumed  constant,  giving  the  re- 
sul  t 

D 

Pr  oc  f  G2  F2  sin  A  n  dD 

D-&D 
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